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Introduction 


u  Fracture  of  advanced  materials,  such  as  intermetallic  alloys,  ceramics,  and  composite 
materials  nearly  always  begins  at  interfaces  of  one  sort  or  another.  Fibrous  fractures,  for 
example,  initiate  at  the  interfacews  between  second  phase  particles  (or  reinforcement  fibers, 
whiskers,  etc)  and  the  ductile  matrix;  the  conditions  for  this  to  occur  are  sensitve  to  both 
matrix  constitutive  properties  as  well  as  to  the  chemistry  and  structure  of  the  interfaces 
between  the  two  phases.  Interfacial  cleavage  involves  the  initiation  of  micrto-cracks  and 
their  propagation  along  interfaces  and  transgranular  cleavage  often  initiates  at  interfacial 
sites,  e.g.  twin  boundaries  in  intermetallic  compounds,  especially  those  with  non-cubic 
crystal  structures.  Here  again,  the  conditions  required  for  this  to  occur  are  sensitive  to 
the  constitutive  properties  of  the  matrix  materials  as  well  as  to  the  interfacial  properties 
themselves.  Intermetallic  alloys,  in  particular,  are  known  *o  possess  complicated  thermal 
and  strain  rate  dependencies  that  can  lead  to  either  briwe  or  ductile  behavior  which  in 
any  event  are  very  much  dependent  upon  the  loading  conditions.  _  > 


Our  research  over  the  past  three  years  has  been  concerned  with  providing  a  theo¬ 
retical  and  comutational  framework  for  analyzing  the  deformation  of  complex  crystalline 
materials,  like  intermetallic  alloys,  and  in  particular  of  the  failure  of  interfaces  within 
these  materials  and  between  these  materials  and  other  phases.  In  addition,  we  have  been 
developing  experimental  methods  to  study  the  failure  of  interfaces  between  metals  and 
metals  and  ceramics.  Our  experimental  studies  are  to  be  directly  correlated  with  our 
theoretical/computational  work.  Specifically,  our  aims  are  to 


1)  conduct  computational,  and  complimentary'  experimental  studies  of  interface  crack 
mechanics,  including  crack  growth,  at  interfaces  separating  crystalline  phases  described  by 
physically  based  slip  theories; 


2)  extend  our  physically  based  slip  theories  to  specifically  account  for  the  dislocation 
micromechanisms  that  are  peculiar  to  more  complex  crystalline  materials  such  as  ordered 
intermetallic  compounds; 

3)  conduct  analytical  and  computational  modeling  of  the  deformation  and  failure  o: 
ordered  intermetallic  alloys,  and  especially  of  the  deformation  and  failure  near  interfaces; 
and 


4)  further  develop  our  vacuum  diffusion  bonding  methods  to  prepare  metal/ ceramic, 
alloy/alloy,  and  intermetalhc/metal  (ceramic)  interfaces  with  controlled  defects. 

Our  ongoing  experimental  work  has  been  concerned  with  the  Nb/A^Oz  interface.  We 
have  used  vacuum  diffusion  bonding  methods  to  prepare  pre-cracked  interfaces  in  this 
system  and  have  been  conducting  experiments  to  measure  the  fractuare  toughness.  We 
have  extended  these  studies  to  the  interfaces  of  N^Al  and  Nb  and  TiAl/Nb  for  the 
purposes  of  measuring  interface  toughness  and  also  to  correlate  the  interface  toughness 
with  ductility  and  toughness  of  composites.  We  have  also  formulated  a  flow  theory  for 
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intermetallic  alloy  materials  that  accounts  for  deviations  from  Schmid’s  rule.  We  have 
previously  shown  that  such  deviations  have  important  implications  for  the  stability  of 
uniform  flow. 

During  the  next  year,  we  plan  to: 

1)  Continue  our  work  of  the  past  year  on  the  analyses  of  the  mechanics  of  cracks  on 
interfaces.  Our  intention  is  to  perform  computational  analyses  of  the  very  near  crack  tip 
fields  under  a  complete  range  of  mixed  mode  loadings  and  of  the  full  fields  for  cracked 
specimen  geometries  of  the  type  we  are  studying  experimentally.  We  intend  to  explore 
the  effects  of  geometric  constraints  on  the  development  of  the  crack  tip  fields,  and  in 
particular,  the  influence  of  interface  geometry  on  the  buildup  of  hydrostatic  stresses.  The 
role  of  hydrostatic  stress  in  determining  the  failure  modes  of  the  interfaces,  and  the  resulting 
ductility  of  the  microstructure  in  specimens  containing  the  interface,  will  also  be  assessed. 

2)  Extend  our  computational  studies  of  localized  deformation  to  ordered  intermetallic 
alloy  systems,  in  particular  Ni3Al  and  later  TiAl  alloys.  As  noted  above,  we  have  for¬ 
mulated  a  constitutive  theory  for  such  materials  and  have  analyzed  the  behavior  in  single 
crystal  Ni$Al.  We  are  now  implementing  the  theory  in  our  FEM  codes  to  perform  analyses 
of  deformation  processes  ;n  these  materials.; 

3)  Residual  stresses,  due  to  incompatible  thermal  strains,  are  known  to  play  an  im¬ 
portant  role  in  the  development  and  propagation  of  flaws  and  cracks  along  interfaces.  We 
have  formulated  a  finite  strain  theory  for  thermal  strains,  and  have  implemented  this  the¬ 
ory7  in  our  finite-element  models  for  composite  behavior,  and  interface  fracture.  We  intend 
to  carry  out  simulations  of  crack  growth  along  interfaces  of  bimaterial  couples  such  as 
Nb/Al203  for  which  the  thermal  strain  mismatch  is  minimal  and  intermetallic  compounds 
such  as  TiAl,Nb$Al ,  and  Ni3Al/Nb  for  which  thermal  mismatch  is  significant;  and 

4)  We  will  continue  our  experiments  on  the  fracture  of  metal/ceramic  and  metal/metal 
interfaces  prepared  by  diffusion  bonding  methods.  We  are  currently  conducting  experi¬ 
mental  studies  of  crack  growth  along  Nb/Ali  O3  interfaces  which  are  prepared  by  vacuum 
diffusion  bonding.  We  will  be  exploring  the  use  of  the  bend  specimen  to  subject  interface 
crack  tips  to  a  wide  range  of  mixed  mode  loadings  in  addition  to  the  compact  tension 
specimen  we  developed  in  the  initial  phase  of  our  project.  We  will  extend  these  studies  to 
interfaces  involving  intermetallic  compound/metal  and  intermetallilc  compound/ceramic 
couples. 

Summary  of  Research 

In  the  past  three  years,  we  have  completed  a  series  of  analyses  on  1)  the  structure 
of  interface  crack  tip  fields;  2)  the  full  field  analyses  of  the  elastic-plastic  fields  associated 
with  interface  cracks  and  flaws;  3)  the  mic^omechanical  and  macromechanical  process  of 
crack  and  flaw  propagation  along  interfaces;  and  4)  deformation  at  the  tips  of  cracks  lying 
on  interfaces  with  structure.  We  have  also  completed  the  formulation  of  a  flow  theory'  for 
ordered  intermetallic  alloys.  Our  computations  have  been  carried  out  for  stationary  cracks 
as  well  as  for  cracks  propagating  along  interfaces  where  the  phenomenology7  of  propagation 
is  governed  by7  adhesive  laws  ascribed  to  the  interface.  Nine  papers  have  been  submitted 
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for  publication.  Four  of  these  have  appeared  in  print  in  1990;  one  additional  has  been 
accepted  for  publication. 

R.J.  Asaro: 

1)  Cracktip  Fields  in  Ductile  Crystals  (with  J.R.  Rice  and  D.E.  Hawk),  International 
Journal  of  Fracture  42,  301  (1990). 

R.J.  Asaro  and  C.F.  Shih: 

2)  Elastic- Plastic  Analysis  of  Cracks  on  Bimaterial  Interfaces:  Part  I  -  Small  Scale 
Yielding,  J.  Appl  Mechanics  55,  299  (1988). 

3)  Elastic- Plastic  Analysis  of  Cracks  on  Bimaterial  Interfaces:  Part  II  -  Structure  of 
Small  Scale  Yielding  Fields ,  J.  Appl.  Mechanics  III,  763  (1990). 

4)  Elastic- Plastic  Analysis  of  a  Collinear  Array  of  Cracks  on  a  Bimaterial  Interface,  J. 
Material  Science  and  Engineering  A  107,  145  (1989). 

R.J.  Asaro,  N.P.  O’Dowd  and  C.F.  Shih: 

5)  Cracks  on  Bimaterial  Interfaces:  Plasticity  Aspects ,  in  proceedings  of  the  Acta/Scripta 
Metallurgica  Conference  on  Bonding,  Structure  and  Mechanical  Properties  of 
Metal/Ceramic  Interfaces,  1989. 

6  Elastic- Plastic  Analysis  of  Cracks  on  Bimaterial  Interfaces:  Part  III  -  Large  Scale 
Yielding,  submitted  to  J.  Appl  Mechanics  (1989). 

7)  Elastic- Plastic  Analysis  of  Cracks  on  Bimaterial  Interfaces:  Part  IV  -  Interfaces  with 
Structure,  submitted  to  J.  Appl  Mechanics  (1990). 

R.J.  Asaro,  N.P.  O’Dowd,  C.F.  Shih  and  A.G.  Varias: 

8)  Failure  of  Bimaterial  Interfaces,  submitted  for  publication,  Mat.  Sci.  and  Eng.  A 
126,  65  (1990). 

9)  Crack  Growth  on  Bimaterial  Interfaces,  in  proceedings  of  the  Acta/Scripta  Metallug- 
ica  Conference  on  Bonding,  Structure,  and  Mechanical  Properties  of  Metal/Ceramic 
Interfaces  (1990) 

A  tenth  manuscript  authored  by  R.J.  Asaro  and  M.  Dao,  entitled  The  Mechanics  of 
Plastic  Flow  in  Ni^Al  is  under  preparation;  the  work  associated  with  it  was  completed 
during  July/August  1990  at  the  University  of  California,  San  Diego. 

Our  initial  studies  under  this  project  were  concerned  with  the  structure  of  plane  strain 
interface  crack  tip  fields.  Specifically  in  Parts  I  and  II  of  our  four  part  series  of  articles  (i.e., 
refs.  2  and  3  listed  above),  the  structure  of  the  caxck  tip  field  that  develop  under  conditions 
pertaining  to  small  scale  yielding  was  determined.  In  reference  4,  the  stress  and  strain  field 
that  evolve  under  conditions  producing  contained  yieldling  and  fully  plastic  conditions 
wre  examined  through  a  study  of  collinear  interface  crack  geometries,  these  (primarily) 
numerical  studies  led  to  an  understanding  of  the  asymptotic  structure  of  interface  crack 
tip  fields  along  with  a  correspondence  of  these  fields  to  the  well  characterized  small  strain 
(HRR)  fields  in  homogeneous  media.  In  particular,  it  was  found  that  the  small  scale 
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yielding  fields  of  the  interface  crack  axe  members  of  a  family  paremeterized  by  a  plastic 
phase  angle ,  £,  defined  as 


where  the  parameters  and  X  axe  defined  vis-a-vis  the  linear  elastic  interface  crack  tip  field 
as  follows. 

In  terms  of  polar  coordinates,  r,  6,  centered  at  the  crack  tip,  the  crack  tip  fields  can 
be  organized  in  the  form 


aij  —  &o 


a<70eor 


!/(»+!) 


h,j(8,r  ;  C,n) 


where  n  is  a  power  law  strainhardening  index  prescribed  in  connection  with  J2  deformation, 
or  J2  flow  theory.  htJ  is  a  bounded  function  which,  however,  varies  with  r(=  r/^j/a^). 
For  — 7r/6  <  f  <  7r/6  however,  htJ  varies  slowly  with  r  so  that  the  near  tip  field  conforms 
closely  with  that  of  an  HER  field.  K  is  the  complex  stress  intensity  factor  defined,  with 
respect  to  the  elastic  crack  tip  asymptotic  field,  as 


atj  = 


\/27r  r 


Re{K 


r,e}  <7 


(0;  e)  +  Im{Kr,e}  e) 


£  is  a  dimensionless  bimaterial  constant  and  L  is  a  crack  dimension.  Our  results  showed 
that  for  cases  involving  two  elastic-plastic  media  that  the  fields,  in  both  materials,  axe 
parts  of  a  single  asymptotic  field;  the  intensity  of  the  stresses  axe  determined  by  the  yield 
strength  of  the  weaker  material  whereas  the  angular  distribution  of  stress  about  the  crack 
tip  is  set  by  the  strainhaxdening  characteristics  of  the  more  weakly  hardening  material. 

In  reference  4,  analyses  involving  contained  yielding  and  fully  plastic  states  were  carried 
out.  An  interesting  discovery  was  that  the  fields  scaled  with  the  value  of  the  J  integral 
when  J  was  normalized  by  yield  strength  and  crack  length.  In  particular,  a  value  J, 
prescribed  as  J/(cr0c0X)  was  found  to  scale  the  fields  that  developed  between  colinear 
arrays  under  conditions  that  produced  deformation  states  ranging  from  small-scale  yielding 
-  to  contained  yielding  -  to  fully  plastic  behavior. 

Finite  deformation  analyses  were  carried  out  in  references  5,  6,  7,  8  and  in  Part 
IV.  Computational  studies  were  performed  using  crystal  plasticity  theory  as  well  as  phe¬ 
nomenological  J2  flow  theory.  Our  finite  strain  calculations  have  included  studies  of  the 
deformation  of  cracks  on  interfaces  with  structure.  An  example  of  this  is  Fig.  1  which 
shows  the  deformed  meshes  for  cracks  lying  on  an  interface  between  two  materials  de¬ 
scribed  by  a  rate  dependent,  finits  strain  J2  flow  theory.  The  materials  on  either  side  of 
the  interface  differed  by  a  factor  of  two  in  yield  strength;  their  properties  were  taken  to 
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match  those  of  the  dissimilar  steels  that  we  studied  experimentally.  The  interface  region 
had  properties  intermediate  between  these  two  phases. 

Figure  1  indicates  the  strong  influence  of  mixed  mode  loading  in  that  our  calculations 
were  carried  out  for  three  values  £  =  —  7t/6,  0,  7t/6.  The  crack  tip  openings,  and  associated 
extents  of  plastic  deformation,  are  greatest  for  the  mixed  mode  sta6tes.  In  fact,  for  a 
given  value  of  the  applied  far  field  J,  the  crack  tip  opening  was  largest  for  the  mixed  mode 
state  corresponding  to  £  =  —  tt/6.  Figure  2  shows  contours  of  hydrostatic  stress  which 
also  illustrate  the  effects  of  mixed  mode  loading.  A  particularly  remarkable  feature  is  the 
shift  in  hydrostatic  stress  with  mode  from  the  softer  to  the  harder  materials  as  the  mode 
shifts  from  that  corresponding  to  £  =  0.  These  particular  calculations  were  part  of  a  set 
of  model  studies  of  our  experiments  concerned  with  dissimilar  steel  interfaces.  One  such 
iinterface  is  shown  in  Fig.  3.  Note  that  this  interface  is  rather  broad  and  is  characterized 
by  an  interdiffusion  zone  of  about  60  /xm  in  thickness.  Our  observations  showed  that,  as  the 
initial  crack  began  to  propagate,  it  quickly  travelled  to  the  boundary  of  the  transition  zone 
and  the  softer  steel.  Given  that  the  fracture  mode  was  fibrous,  this  is,  in  fact,  expected 
since  our  computations  showed  that  the  hydrostatic  stresses  are  larger  in  the  transition 
zone  toward  the  softer  steel.  Fibrous  fractures,  which  are  driven  by  void  initiation  and 
growth,  are  accelerated  by  hydrostatic  tensile  states  of  stress. 

Recent  Work 

Our  more  recent  work  has  been  concerned  with  the  formulation  of  a  flow  theory' 
for  ordered  intermetallic  alloys.  NijAl  was  proposed  as  an  initial  example.  One  of  the 
particularly  interesting  features  of  the  dislocation  mechanics  in  such  crystalline  materials  is 
that  they  display  strong  departures  from  Schmid’s  rule  of  a  critiacal  resolved  shear  stress. 
Deviations  from  the  Schmid  rule  of  this  type  represent  non-normality  in  the  flow  rule. 
Asaro  and  Rice  have  shown  that  there  are  indeed  important  implications  of  non-normality 
of  ths  type  regarding  stable,  uniform,  plastic  flow. 

As  part  of  this  project,  Dao  and  Asaro  have  formulated  a  flow  rule  for  both  rate 
independent  and  rate  dependent  materials.  For  rate  independent  materials,  the  flow  rule 
takes  the  form  (for  each  slip  system) 


7m,— Q;f 

where  r  is  the  shear  stress  phrased  on  slip  system  coordinates  as  shown  in  Fig.  4.  This 
figure  also  indicates  some  typical  micromechanical  processes,  involving  dislocations,  that 
give  rise  to  the  non-  schmid  effects.  Q  is  the  dyad 

Q  =  sm  +  a  :  H  :  L-1 

where  a  is  the  matrix  of  non-Schmid  factors,  i.e. 
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For  rate  dependent  materials,  the  flow  law  instead  takes  the  form 

7  =  function(Q  :  <r;  material  state). 

Dao  and  Asaro  have  analyzed  the  behavior  of  Ni^Al  within  this  framework  and  have 
calculated,  based  on  experimental  data,  values  for  their  non-Schmid  coefficients. 

We  are  currently  in  the  process  of  implementing  our  theory  in  one  of  our  FEM  codes. 
Our  initial  objectives  will  be  to  study  the  development  of  flow  patterns  in  single  and 
polycrystalline  microstructures  of  these  materials  and  to  then  analyze  the  behavior  of 
interfaces  between  Ni$Al  and  Nb. 


a,z  \ 
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List  of  Participants  and  Status 

S.  Schoenfeld,  Graduate  Student 
M.  Dao,  Graduate  Student 

Both  of  these  students  are  in  their  second  year. 

Other  Sponsored  Research 

’’Fundamental  Studies  of  Localized  Plastic  Deformation”,  NSF,  71,900. 

Awards 

R.J.  Asaro  has  been  awarded  the  TMS  Champion  H.  Mattewson  Gold  Metal  for  1991 
for  his  research  on  micromechanics. 
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Robert  J.  Asaro,  a  professor  in  the  Department  of  Applied 
Mechanics  and  Engineering  Sciences  (AMES)  at  the  University  of 
California,  San  Diego  (UCSD) ,  has  been  awarded  the  TMS  Champion 
H.  Mathewson  Gold  Medal  for  1991. 

The  Mathewson  Gold  Medal  is  presented  to  an  author,  or 
authors,  of  a  paper  published  in  "Metallurgical  Transactions," 
which  is  considered  to  be  the  most  notable  contribution  to 
metallurgical  science  during  the  period  under  review.  The  award 
was  established  in  1950  to  honor  Dr.  Champion  H.  Mathewson,  a 
noted  scientist  and  teacher. 

Formal  presentation  of  the  award  will  be  made  during  the 
annual  dinner  of  the  Minerals,  Metals  and  Materials  Society,  to 
be  held  Feb.  19,  1991,  at  the  Radisson  Suite  Hotel  in  New 
Orleans,  La. 

The  Minerals,  Metals  and  Materials  Society  is  a  member 
society  of  the  American  Institute  of  Mining,  Metallurgical  and 
Petroleum  Engineers,  Inc. 
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Vbstract.  Results  on  I  lie  .isvmph'lic  analysis  of  track  lip  lields  in  elastic-pla  -lie  single  crystals  arc  presented  and 
Mime  preliminary  results  of  limtc  dement  solutions  lor  cracked  solid'  of  this  type  are  summari/ed  In  the  eases 
'i  ndicd .  ins  oh  me  plane  strain  tensile  and  anti-plane  shear  cracks  in  ideally  plastic  I  c  c  and  h  c  c  cry  stals.  analyzed 
within  conventional  small  displacement  gradient  assumptions,  the  asymptotic  analyses  re\eal  striking  discontinuous 
fields  at  the  crack  lip 

for  the  stationary  crack  the  sttvss  state  is  lotind  to  he  locally  uniform  in  each  of  a  fain 1 1>  of  angular  sectors  al 
tile  crack  tip.  hut  to  jump  discontinuous!)  at  sector  boundaries,  which  arc  also  the  surfaces  of  shear  discontinuities 
in  the  displacement  held,  l  or  the  quasi-statically  growing  crack  the  stress  state  is  lull)  continuous  from  one 
near-tip  angular  sector  to  the  nest,  hut  now  some  of  the  sectors  involve  elastic  unloading  from,  and  reloading  to. 
a  Yielded  state,  and  shear  discontinuities  of  the  velocity  held  develop  at  sector  boundaries.  In  all  anti-plane  case 
studied,  inclusion  of  inertial  terms  for  (dynamically  l  growing  cracks  restores  a  discontinuous  stress  field  al  the  tip 
which  mines  through  tile  material  as  an  elastic-plastic  shock  wave  for  high  symmetry  crack  orientations  relative 
lo  the  crystal.  the  discontinuity  surfaces  are  sometimes  coincident  with  the  active  crystal  slip  planes,  hut  as  often 
lie  perpendicular  to  the  family  of  active  slip  planes  so  that  the  discontinuities  correspond  tv'  a  kinking  mode  of 
shear  . 

I  he  finite  clcmcntVtudics  so  far  attempted,  simulating  the  ideal!)  plastic  material  model  in  a  small  displacement 
gradient  tv  pe  program.  appear  to  he  consistent  w  ith  the  as)  mptotic  analyses,  Small  scale  y  lelding  solutions  confirm 
the  expected  discontinuities,  within  limits  of  mesh  resolution,  of  displacement  for  a  stationary  crack  and  ol  velocity 
lor  quasi-static  growth. Yurther.  the  discontinuities  apparently  extend  well  into  the  near-tip  plastic  /one  \  finite 
element  formulation  suitable  for  arbitral)  deformation  has  been  used  to  solve  for  the  plane  strain  tension  of  a 
I  ay lor ■-hardening  crvstal  panel  containing,  a  center  crack  with  an  initially  rounded  tip.  I  his  shows  elfects  due  to 
lattice  rotation,  which  distinguishes  the  regular  versus  kinking  shear  modes  of  crack  up  relaxation,  and  holds 
promise  for  exploring  the  mechanics  of  crack  opening  al  the  tip. 


I.  Introduction 

This  paper  summarizes  recent  analytical  and  numerical  investigations  into  the  nature  of  the 
near-crack-tip  stress  and  deformation  fields  in  ductile  single  crystals.  Ductile  crystals  deform 
plastically  by  the  motion  of  dislocations  on  a  limited  set  of  slip  systems.  A  continuum 
representation  of  this  plastic  deformation  consistent  with  the  Schmid  rule,  which  stales  that 
How  on  it  sy  stem  is  activated  w  hen  the  shear  stress  resolved  on  that  system  reaches  a  critical 
value,  is  used  in  the  analyses  to  be  presented.  This  formulation  leads  to  a  yield  surface  in 
stress  space  consisting  of  planar  facets  joined  at  vertices  and  to  an  "associated”  plastic 
straining  relation. 

General  methods  of  constructing  asymptotic  near-tip  fields  for  such  crystals,  with  either 
stationary  or  quasi-statically  growing  cracks,  have  been  obtained  in  the  ideally  plastic  case 
for  both  anti-plane  strain  (mode  III:  |l|)  and  tensile  plane  strain  (mode  I:  [2] )  cracks.  The 
results,  as  illustrated  for  common  crack  orientations  in  fee  and  bee  crystals,  lead  to 
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striking  predictions  of  discontinuities  at  the  crack  tip.  Full  scale  elastic-ideally  plastic 
solutions  to  the  near-tip  stress  and  deformation  fields  have  been  given  for  stationary  cracks 
in  the  mode  III  study.  These  show  that  all  How  is  confined  to  planar  plastic  zones  emanating 
from  the  crack  tip.  across  which  both  displacement  and  stress  are  discontinuous.  Asymptotic 
analy  sis  of  dynamic  crack  growth,  i.e.,  including  inertia,  has  ’'een  developed  as  well  for  that 
mode  [3]  and  reveals  that  an  elastic-plastic  shock  discontinuity  moves  along  with  the  tip.  In 
addition,  asymptotic  fields  of  the  HRR  type  have  been  developed  for  stationary  cracks  in 
crystals  showing  Taylor  hardening,  with  a  power-law  stress-strain  relation  at  large  strain,  in 
mode  III  [4]  and  mode  I  [5]. 

The  mode  I  asymptotic  analysis  [2],  based  on  ideal  plasticity  and  a  "small  displacement 
gradient"  formulation,  shows  that  for  material  at  yield  the  stress  state  is  constant  within 
angular  sectors  whose  boundaries  are  certain  crystallographic  directions  on  which  discon¬ 
tinuities  in  either  displacement  (stationary  crack)  or  velocity  (quasi-statically  growing  crack) 
are  possible.  A  direct  comparison  of  the  mode  I  analyses  is  made  here  with  the  numerical 
results  of  Hawk  and  Asaro  [6].  When  comparing  different  types  of  crystals,  or  crack 
orientations  within  a  given  crystal,  the  structure  of  the  dislocations  necessary  to  produce  the 
same  continuum  field  discontinuity  is  different.  Furthermore,  certain  dislocations  structures 
may  induce  rotation  of  the  lattice  relative  to  the  material,  thus  changing  the  resolved  shear 
stresses  on  slip  systems  and  causing  a  geometric  hardening  or  softening  of  the  crystal.  This 
is  important  particularly  when  large  deformations  are  taken  into  account. 

Recently,  full  scale  elastic-plastic  solutions  for  mode  1  cracks  in  ductile  crystals  have  been 
obtained  numerically  by  the  detailed  finite  element  analyses  of  Hawk  and  Asaro  [6].  Some 
of  their  results  are  briefly  summarized  here  for  several  different  loading  cases.  These  analyses 
model  the  constitutive  behavior  of  the  crystal  with  a  visco-plastic  formulation  in  the  nearly 
rate-independent  limit.  The  numerical  analysis  of  Asaro  [6]  models  from  small-scale  to 
general  yielding  a  center  cracked  panel  with  a  blunted  crack  tip  under  uni-axial  tension.  The 
slip  systems  of  the  crystal  are  idealized  by  a  planar  double-slip  model.  The  effect  of  Taylor 
type  hardening  of  the  slip  systems  and  large  displacement  gradients  (e.g..  lattice  rotation) 
are  included.  In  the  small-scale-yielding  analysis  of  Hawk  [6],  both  a  stationary  and  quasi- 
statically  propagating  perfectly  sharp  tensile  crack  are  simulated.  The  crystal  is  modeled  as 
elastic-ideally  plastic  using  a  complete  description  of  the  slip  systems  in  a  fee  crystal.  Small 
displacement  gradients  are  assumed. 


2.  Constitutive  law 

The  crystals  considered  can  undergo  both  elastic  and  plastic  deformation.  The  plastic 
deformation  is  consistent  with  a  continuum  description  of  single  crystals  [7-10].  The  total 
strain-rate  is  taken  as  the  sum  of  the  elastic  strain-rate  and  plastic  strain-rate 

i:„  =  k;;  4-  ( 1 ) 

Under  plane  strain  conditions,  which  arc  possible  for  the  high-symmetry  crack  orientations 
considered  here.  rx,,  and  <x:i  are  zero.  The  plastic  deformation  of  the  crystal  occurs 

by  the  motion  of  dislocations  along  certain  preferential  crystallographic  planes.  The  move¬ 
ment  of  these  dislocations  causes  a  permanent  dilationless  straining  of  the  crystal.  A  slip 
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system  denoted  by  y.  is  defined  by  two  unit  vectors  giving  the  slip  direction  s'”  and  the  normal 
to  the  slip  plane  n'11.  The  parameter  ym  describes  the  amount  of  shear  strain  on  each  slip 
system.  These  preferential  slip  systems  defined  by  s'”  and  nm  vary  according  to  the  crystal 
structure  (e.g.,  face  centered  cubic,  f  c  e.  or  body  centered  cubic,  bee).  It  is  possible  to  express 
the  plastic  strain  rate  as  a  sum  over  all  the  \  slip  systems  of  the  crystal  as 


Here  /fj'  is  termed  the  Schmid  factor  and  is  determined  from  the  slip  direction  s'”  and  the 
normal  to  the  slip  plane  nm  as 

+  .v;”//,'”)  <3> 

The  resolved  shear  stress,  r1”,  on  a  system  is  expressible  in  terms  of  these  Schmid  factors  as 


The  plastic  deformation  of  a  crystal  is  said  to  obey  the  Schmid  rule  if  the  instantaneous 
shear  response  y”'  of  any  given  system  i  depends  on  the  current  stress  state  only  through  the 
combination  in  (4).  In  a  time-independent  plasticity  formulation,  a  necessary  but  not 
sufficient  condition  for  slip  to  occur  is  that  the  resolved  shear  stress  rm  reaches  a  critical  value 
g'”.  This  criteria  results  in  a  yield  surface  in  stress  space  w  hich  consists  of  planar  facets  that 
join  at  vertices  at  which  two  or  more  systems  are  simultaneously  active  Equation  (2)  then 
assures  that  the  flow  rule  is  of  an  "associated"  type.  Rice  [II]  showed  that  for  r/g/</-plastic 
incompressible  solids  of  this  class  in  plane  flow,  the  yield  surface  can  be  represented  as  a 
curve  in  the  reduced  stress  space  of(<rM  -  <x::)/2  and  ar_.  Such  a  yield  surface  is  schematic¬ 
ally  shown  in  Fig.  I  and  for  crystals  as  described  here,  it  is  polygonal.  In  *fn  c/W/V-plastic 


I  iv.  I.  Schematic  of  yield  surface  lor  a  ductile  crystal  undergoing  plane  straining. 
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material  the  same  result  applies  for  stress  states  in  large  sustained  plastic  strain.  However, 
in  such  a  case  limited  deformation  may  occur  on  slip  systems  corresponding  to  segments  of 
the  yield  surface  sensitive  to  the  value  of  (an  +  it,, )  but  which  cannot  produce  sustained 
plane  flow.  We  will  term  the  sides  of  the  polygon  as  "flats"  and  the  point  of  junction  of  two 
sides  as  a  "vertex”.  The  flow  law  of  (2)  being  associative,  the  direction  of  plastic  straining, 
with  components  (fjj  —  f:{j )/ 2  and  f:('; ,  is  normal  to  the  yield  surface  on  a  "flat"  and  within 
the  forward  fan  of  normals  at  a  "vertex". 

It  will  also  be  convenient  to  represent  the  slip  systems  of  the  crystal  as  traces  in  the  .v, .  v, 
plane.  For  example,  in  Fig.  2a.  a  unit  cube  representing  the  lattice  of  an  fee  crystal  is  shown. 
The  slip  planes  of  a  fee  crystal  are  the  1111}  planes  indicated  by  the  various  shaded  planes 
in  the  cube.  The  slip  directions  are  <101)  directions  which  are  the  diagonals  of  the  cube 
faces.  If  the  crack  lies  in  the  plane  (0  1  0)  and  the  crack  tip  lies  along  the  diagonal  of  the  cube 
face  as  indicated,  then  the  traces  of  the  slip  systems  (i.e.  intersection  of  the  slip  planes  with 
the  .v, .  .\\  plane)  are  as  indicated  in  Fig.  2b.  It  is  then  possible  to  represent  the  slip  direction 
s  and  slip  plane  normal  n  as  the  directions  S  and  N  respectively  in  the  .v, .  ,v2  plane  as  shown 
in  the  figure.  Only  those  traces  shown  correspond  to  systems  which  can  produce  sustained 


Fig.  2.  (a)  Face  centered  cubic  slip  systems,  and  orientation  of  the  crack  plane  and  crack  tip  for  case  illustrated 
below,  (bl  Crack  on  ((I  I  0)  plane  growing  in  [1  0  I]  direction.  Traces  are  shown  of  those  slip  plane  families  which 
can  accommodate  sustained  plane  straining. 
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plane  flow.  The  solid  line  traees  correspond  to  individual  J I  I  1 1  planes,  and  the  dashed  lines 
to  two  such  planes  which  intersect  along  a  <  1  1  ()>  slip  direction  common  for  each  and  which, 
by  equal  coincident  slip,  produce  plane  flow. 

The  finite  element  analyses  use  a  visco-plastic  formulation  where  slip  can  occur  on  any 
system  so  long  as  resolved  shear  stress  on  that  system  is  non-zero.  A  simple  power-law 
relation  [12-14]  can  be  used  to  describe  the  rate  of  slip  on  each  system  as 

=  iiV  sgn  (r1  )  ^  —  J  (» 

As  the  exponent  m  -»  0  the  visco-plastic  model  approaches  the  rate-independent  formu¬ 
lation.  Use  of  a  visco-plastic  law  in  a  finite  element  formulation,  pioneered  in  the  work  of 
Pierce  et  al.  [14],  eliminates  certain  problems  of  uniqueness  and  considerations  of  when  a 
point  is  at  yield.  The  reflects  the  current  level  of  strain  hardening  in  the  crystal.  In  [6]. 
results  are  presented  for  both  ideal  plasticity  i.e.,  g'"  is  a  constant.  r„,  and  cases  where  gm 
is  a  function  of  the  sum  of  the  magnitudes  of  the  slips  (T  =  E,  |‘/”|).  which  coincides  with 
Taylor  hardening. 


3.  Asymptotic  analysis 

The  asymptotic  analyses  (e.g..  lim,  <xj/((r.  0.  /))  of  Rice  [2]  are  summarized  here  for  the  case 
of  ideal  plasticity  in  a  "small  displacement  gradient"  formulation  of  the  mode  1  problem. 
Based  on  the  equations  of  equilibrium  with  a  bounded  crack  tip  stress  state  and  with  the 
further  condition  that  the  stress  state  is  at  yield  relative  to  a  slip  plane  trace  in  the  direction 
S,,  the  following  requirement  is  found: 


(A',e,  )(S,c,  )(a'u  +(?::)  =  0 


(6) 


Greek  indices  range  from  1  to  2.  As  indicated.  Fig.  2b.  the  e7  are  the  components  of  the  radial 
unit  vector  in  the  .v, .  as  plane;  (e, .  e: )  =  (cos  0,  sin  0).  The  terms  S ,  and  N,  are  the 
components  of  the  traces  of  the  slip  direction  and  slip  plane  normal  as  previously  defined. 
The  <t'(,  denote 


lim 

r  .0 


<  <r,iM.  0) 

CO 


(7) 


Equation  (6)  implies  rr J ,  +  ad  =  0  for  all  0  except  for  four  special  values  when  e  is  aligned 
with  either  N  orS.  Based  on  the  form  of  the  stresses  consistent  with  equilibrium  for  material 
at  yield,  this  statement  further  implies  that  within  sectors  bounded  by  these  four  special 
values  of  0,  all  a'7ll  are  zero;  i.e.  the  stress  state  is  independent  of  0.  For  a  tensile  crack  it  is 
shown  that,  in  order  to  meet  boundary  conditions,  either  (i)  the  stresses  in  certain  angular 
sectors  around  the  crack  tip  arc  not  at  yield  or  (ii)  that  the  stresses  change  discontinuously 
at  the  special  values  of  0  as  mentioned  above. 
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3.1.  Stationary  crack 

For  the  stationary  crack,  as  r  -»  0,  the  yield  condition  can  be  met  in  all  angular  sectors 
around  the  crack  tip,  thus  requiring  discontinuous  jumps  in  the  stress  state.  The  stress 
discontinuity  must  be  from  one  point  on  the  yield  surface  to  another,  and  considering  the 
nature  of  restrictions  on  such  discontinuities,  this  path  must  be  a  straight  line  in  the  stress 
space  and  the  line  must  lie  everywhere  along  the  yield  locus  (and  hence  correspond  to  a  flat 
segment  in  Fig.  1 )  since  otherwise  elastic  unloading  would  occur.  The  solution  which  satisfies 
the  jump  conditions  for  the  stress  discontinuity  is  one  where  the  stress  state  changes  from 
vertex,  to  vertex  on  the  yield  surface.  The  deformation  fields,  consistent  with  the  above  stress 
state,  must  have  a  shear  type  of  discontinuity  in  the  displacement  field  along  the  same 
direction,  e.  for  which  the  stress  is  discontinuous.  These  correspond  to  concentrated  shear 
parallel  to  the  slip  plane  traces  when  e  is  aligned  with  S.  and  to  kink-like  shear  perpendicular 
to  the  traces  when  e  is  aligned  with  N.  The  specifics  of  constructing  such  solutions  may  be 
found  in  [2],  Specific  examples  of  these  solutions  will  be  given  in  the  discussion  of  the  finite 
element  solutions  of  [6], 

3.2.  Growing  crack 

For  the  quasi-statically  growing  crack,  Drugan  and  Rice  f  1 5]  demonstrated  for  this  type  of 
material  that  discontinuities  in  stresses  and  displacements  cannot  exist.  However,  it  is 
possible  to  have  velocity  discontinuities.  It  is  therefore  necessary  to  have  bordering  on  the 
plastic  sectors  angular  sectors  in  which  elastic  unloading  (and  perhaps  reloading)  occur. 
Construction  of  the  solution  for  a  particular  crystal  orientation,  again  for  which  the  details 
are  given  in  [2],  shows  a  complicated  pattern  of  plastic  sectors  and  elastic  unloading  and 
reloading  sectors  is  necessary  to  model  the  crack  tip  fields  for  quasi-static  growth.  A  specific 
example  will  be  presented  for  comparison  with  the  finite  element  solutions  of  Hawk  [6]  in 
the  next  section  which  shows  that  one  of  the  physically  active  system  for  the  stationary  crack 
becomes  inactive  in  the  growing  case. 

3.3.  Note  on  role  of  S  and  I\ 

In  the  asymptotic  analyses  presented,  based  on  the  "small  displacement  gradient"  formula¬ 
tion,  the  role  of  the  slip  direction  S  and  slip  plane  normal  N  are  interchangeable  since  they 
come  in  through  the  symmetric  form  of  (3).  In  certain  cases  [2]  fee  crystals  and  bcc  crystals 
differ  only  by  an  interchange  of  S  and  N.  Therefore  solutions  obtained  for  such  fee  crystals 
are  also  valid  for  bcc  crystals  as  well  (aside  from  a  slight  scaling).  However,  the  structure 
of  the  actual  dislocations  necessary  to  produce  the  shear  deformation  is  quite  different  when 
S  and  N  are  reversed.  For  example,  shown  in  Fig.  3a  and  3b  is  shear  parallel  to  the  slip  planes 
and  shear  perpendicular  to  the  slip  planes  respectively.  The  motion  indicated  in  Fig.  3a  can 
be  generated  by  the  emission  of  dislocations  from  the  crack  tip  along  the  slip  direction  as 
shown  in  Fig.  3c.  Notice  however  that  in  Figure  3b  the  slip  planes  must  form  a  kink  in  order 
to  accommodate  the  deformation.  Experiments  in  Fe-3%  Si  crystals  [16]  show  results 
consistent  with  kink-like  shear.  This  kink  requires  dislocation  dipoles  as  in  Fig.  3d.  These 
dipole  loops,  illustrated  as  pairs  of  dislocations,  cannot  be  swept  out  from  the  crack  tip  and 
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Fig.  3.  (a)  Shear  in  band  lying  parallel  to  the  active  slip  systems,  (b)  Shear  by  a  kinking  mode  in  a  band  lying 
perpendicular  to  the  active  system,  (c)  Dislocations  can  be  generated  at  the  crack  tip  and  swept  out  along  slip 
planes,  or  can  be  generated  front  internal  sources,  to  produce  the  slip-plane-parallel  shear  band  of  (a),  (d) 
Dislocations  dipole  loops  must  be  nucleated  from  internal  sources  and  expand  as  illustrated  to  produce  the  kinking 
shear  band  of  (b). 


must  therefore  depend  on  the  availability  of  internal  sources  to  produce  them.  Requiring  the 
availability  of  internal  sources,  as  opposed  to  allowing  dislocations  to  be  swept  out  from  the 
crack  tip  along  slip  planes,  may  influence  whether  the  fracture  of  a  particular  crystal  is  brittle 
or  ductile.  It  is  seen  also  that  the  lattice  between  the  pair  of  dislocations  in  Fig.  3d  is 
apparently  rotated.  This  can  cause  either  a  geometric  hardening  or  softening  of  the  material 
because  the  resolved  shear  stress  tm  on  each  slip  system  varies  as  the  Schmid  factor,  p'f . 
changes  due  to  rotation  of  S  and  N  directions.  The  effect  of  such  lattice  rotation  will  be  seen 
later  in  the  large  deformation  results  of  Asaro  [6]. 

3.4.  Some  further  analyses  of  mode  III  cracks 

The  analyses  discussed  so  far  have  been  for  tensile  cracks,  which  is  the  physically  more 
interesting  case.  However  they  are  approximate  in  that  their  validity  is  only  in  the  limit  as 
r  ->  0.  The  simpler  nature  of  the  equations  for  mode  111  (anti-plane  strain)  type  loading 
allows  for  more  complete  solutions  [1.3]  which  we  briefly  summarize  here.  Three  such 
solutions  are  shown  in  Fig.  4  exclusively  for  the  fee  crystal  with  the  orientation  of  the  crack 
the  same  as  in  Fig.  2a  with  the  crack  along  the  (0  1  0)  plane  with  crack  tip  along  the 
[10T]  direction.  Under  mode  III  loading  conditions  the  only  non-zero  displacement  is 
m,(.Y|  .  -Y:)  which  gives  rise  to  only  two  non-zero  stress  components  <7,,  and  cr,:.  In  this  case 
only  the  solid-line  slip  plane  traces  in  Fig.  2b  can  provide  anti-plane  straining.  The  yield 
surface  can  thus  be  represented  in  <t,,.  <t,:  stress  space  and  for  this  orientation  is  diamond 
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Fig.  4.  Summary  of  analyses  for  anti-plane  shear  of  a  cracked  ideally  plastic  crystal:  (a)  Stationary  crack;  (b) 
Quasi-static  crack  growth;  (c)  Dynamic  crack  growth  asymptotic  field  (inertial  effects  included). 

shaped.  The  assembled  stress  sectors  for  these  solutions  will  be  labeled  with  A.  B.  etc.  with 
corresponding  points  or  trajectories  labeled  on  the  yield  surface.  As  is  often  the  case,  similar 
features  are  seen  to  carry  over  from  the  mode  111  solution  to  the  mode  1  solution. 

The  stationary  case  is  shown  in  Fig.  4a  with  the  corresponding  yield  surface  and  stress 
trajectories  as  indicated.  Only  two  sectors  A.  B  in  the  upper  half  plane  (a  corresponding 
sector  to  B  labeled  B'  is  shown  in  the  lower  half  plane)  of  the  constant  stress  type,  as 
previously  discussed  for  plane  strain  in  the  asymptotic  limit  as  r  -+  0.  are  necessary  to 
construct  the  full  solution.  Both  of  the  corresponding  stress  points  are  on  the  yield  surface 
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as  indicated.  These  sectors  are  separated  by  a  line  of  discontinuity  in  both  displacement  and 
stress.  The  stress  discontinuity  can  be  thought  of  as  a  rapid  transition  (in  a  vanishing  small 
sector)  from  point  A  to  point  B  as  indicated  by  the  arrow  on  the  yield  surface.  The  complete 
solution  of  Rice  and  Nikolic  [I]  shows  that  these  sectors  are  actually  elastic  with  all  plasticity 
collapsed  into  a  shear  zone  along  the  line  of  the  displacement  discontinuity.  The  size  of  this 
plastic  zone  is  found  exactly  in  their  analysis. 

The  quasi-static  crack  growth  case  is  shown  in  Fig.  4b.  Three  sectors  A.  B'  and  C"  in  the 
upper  half  plane  (w  ith  corresponding  sectors  B  and  C  in  the  lower  half  plane)  are  necessary 
to  assemble  the  solution.  Sector  A  is  of  the  constant  stress  type  with  the  state  of  stress  at  the 
yield  surface  as  indicated  on  the  yield  surface  plot.  As  necessary  from  the  asymptotic  solution 
this  sector  is  bounded  by  an  elastic  sector  C .  Sector  A  is  separated  from  C'  by  a  velocity 
discontinuity  (expected  to  extend  over  a  finite  size  region  based  on  the  approximate  analysis 
of  Rice  and  Nikolic  [1])  along  which  all  the  plastic  flow  takes  place  in  the  region  forward  of 
the  crack  tip.  As  shown  by  the  lightly  shaded  region  in  Fig.  4b.  as  the  crack  grows  this 
velocity  discontinuity  also  leaves  behind  a  wake  of  plastic  deformation.  Sector  C'.  as 
indicated  by  the  corresponding  path  within  the  yield  surface,  unloads  elastically  and  reloads 
to  yield  by  the  time  it  reaches  the  other  plastic  sector  B'  creating  a  thin  wedge  of  plastic 
deformation  along  the  crack  face,  shown  in  darker  shading.  The  exact  details  of  the 
construction  of  the  angles  for  these  sectors  is  presented  in  [1], 

Finally,  the  case  of  dynamic  crack  growth  where  the  effects  of  inertia  are  taken  into 
account  is  show  n  in  Fig.  4c  [3].  In  this  case  the  entire  field  consists  of  constant  stress  sectors 
separated  by  an  elastic-plastic  shock  (possible  only  in  single-crystal-like  materials  with  flat 
segments  along  their  yield  locus).  Both  velocity  and  stress  are  discontinuous  across  the 
shock.  Remarkably,  the  strain  accumulated  in  crossing  the  shock  is  finite  at  the  crack  tip. 
It  is  of  the  order  of  the  strain  at  first  yield  divided  by  the  elastic  Mach  number  associated 
with  the  crack  speed.  The  angle  1 1/,  in  Fig.  4c.  is  proportional  to  the  Mach  number  at  low 
speed  and  approaches  {n  -  ()„)  2  at  the  sonic  speed.  There  must  be  a  complicated  transition 
between  the  near  tip  field  with  the  inertia  included.  Fig.  4c.  and  that  without  Fig,  4b.  This 
is  an  example  of  nonuniform  asymptotic  limits  as  one  considers  r  ->  0  and  %/(w:!G  -»  0 
where  a  is  the  density.  G  is  the  shear  modulus  and  v  is  the  crack  speed.  E.g..  the  quasi-static 
growth  case  may  be  considered  as  letting  q  -*  0.  or  v  -*  O' .  before  letting  r  —>  0. 


4.  Comparison  of  asymptotic  and  numerical  analyses  for  mode  I 

4.1.  Small-scale-yielding  results 

Under  small-scale-yielding  conditions,  a  stationary  and  quasi-statically  propagating  plane 
strain  mode  I  crack  in  a  fee  crystal  has  been  investigated  by  the  finite  element  analysis  of 
Hawk  [6],  The  material  is  elastic-ideally  plastic,  modeled  by  a  visco-plastic  formulation  as 
in  (5).  with  g1'1  =  t„  and  m  =  0.005.  Small-scale-yielding  conditions  exist  when  the  size  of  the 
plastic  zone  is  much  smaller  than  the  region  over  which  the  elastic  singular  stress  field 
dominates.  The  dominance  of  the  elastic  singularity  allows  the  crack  tip  region  to  be  modeled 
as  an  infinite  solid  with  a  semi-infinite  crack  where  the  stresses  approach  those  of  the  elastic 
singularity  as  r  -»  x .  The  finite  element  mesh  of  both  analyses  is  shown  in  Fig.  5.  Sinee  the 
mode  I  problem  is  sy  mmetric,  only  the  upper  half  plane  is  modeled.  The  finite  element  mesh 
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Fig.  5.  Finite  element  mesh  for  "small  displacement  gradient"  analysis  of  small-scale  yielding  in  an  ideally  plastic 
crystal. 

is  shown  in  three  sections  A,  B  and  C  for  clarity.  Section  A  fits  into  the  rectangular  space 
of  section  B  and  section  B  in  turn  fits  into  the  space  in  section  C.  The  ratio  of  the  size  of  the 
smallest  rectangular  element  of  region  A  to  the  outermost  radius  of  region  C  is  approximately 
10  4.  Each  quadrilateral  element  is  actually  a  so  called  cross  triangle  element  [17]  made  up 
of  four  constant-strain  triangles  (formed  from  the  diagonals  of  the  element)  which  as  a  group 
behave  well  under  incompressible  conditions.  The  crack  opens  to  the  left  with  nodes  ahead 
of  the  crack  tip  constrained  from  vertical  movement  by  the  symmetry  boundary  condition. 
Tractions  corresponding  to  the  elastic  K.  field  singularity  are  applied  to  the  outer  boundary 
of  section  C.  The  results  are  primarily  from  section  A  and  the  region  immediately  surround¬ 
ing  it.  The  details  of  the  finite  element  method  are  presented  in  [6]. 

The  analysis  of  the  stationary  crack  in  an  fee  crystal  was  performed  with  the  orientation 
of  the  crack  the  same  as  shown  in  Fig.  2a.  The  crack  lies  in  the  (0  I  0)  plane  with  the 
crack  tip  along  the  [1  01]  direction.  Elasticity  was  modeled  as  isotropic.  The  yield  surface 
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Fig.  6.  Stationary  crack  in  ideally  plastic  crystal  of  Fig.  2:  (a)  Yield  surface:  (b)  Asymptotic  structure  of  crack  tip 
field;  (c)  Finite  element  results:  grey  zone  with  T  ^  0.01  (r„  O')  essentially  denotes  plastic  zone;  black  zone  with 
T  S  10  (r„  G)  shows  zone  of  more  concentrated  plastic  strain. 

(appropriate  to  sustained  flow)  for  this  orientation  is  shown  in  Fig.  6a.  The  asymptotic 
solution  obtained  by  Rice  [2]  for  the  stationary  crack  case  is  shown  in  Fig.  6b.  In  the 
stationary  case,  the  entire  crack  tip  region  may  deform  plastically.  For  the  upper  half  plane, 
the  solution  consists  of  four  constant  stress  sectors  labeled  A.  B.  C,  and  D  separated  by  stress 
and  displacement  discontinuities  as  indicated  in  Fig.  6b.  The  stress  state  of  each  sector  is 
similarly  labeled  as  points  on  the  yield  surface  A.  B.  C  and  D  in  Fig.  6a.  The  discontinuity 
between  A  and  B  at  54.7°  and  the  one  between  C  and  D  at  125.3°  correspond  to  the  slip 
directions  for  the  solid  traces  in  Fig.  2b.  The  discontinuity  between  sectors  B  and  C 
corresponds  to  the  normal  of  the  slip  plane  indicated  by  the  dashed  traces  in  Fig.  2b.  Each 
discontinuity  in  stress  represents  a  rapid  transition  from  point  to  point  on  the  yield  surface 
as  shown  by  the  arrows. 

Let  the  sum  over  the  slip  systems  of  the  magnitude  of  slip.  Z,  |-/”|  be  denoted  as  T.  This 
quantity  is  an  overall  measure  of  the  amount  of  plastic  straining  and  is  shown  in  Fig.  6c 
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normalized  by  the  elastic  strain  at  yield,  t „/6\  The  grey  zone,  formed  by  shading  individual 
triangular  elements  w  hose  value  of  r/(r0  G )  is  between  0.01  and  10,  is  representative  of  the 
overall  plastic  zone  shape.  The  black  bands  are  those  elements  whose  value  of  T  is  greater 
than  10.  Therefore,  the  most  intense  straining  is  along  the  predicted  discontinuities  of  the 
asymptotic  solution.  It  is  interesting  to  note  that  all  the  sectors  are  predicted  to  be  stressed 
to  yield  as  r  -*•  0,  but  only  the  interfaces  between  them  arc  proven  to  be  deforming  plastically 
in  the  asymptotics.  The  finite  element  solution  indicates  that  sector  D  bounding  the  crack 
face  is  relatively  free  of  plastic  deformation  while  the  others  deform.  Some  plastic  strain 


Fix.  7-  (a)  Deformed  mesh  for  stationary  crack  (factor  of  50  amplification  of  displacements);  note  shear  zones 
coincident  with  directions  shown  in  Fig.  6(b).  <b)  Line  lengths  proportional  to  displacements  accumulated  in  a 
load  increment:  suggests  that  displacement  discontinuities  predicted  by  asymptotic  analysis  extend  well  into  plastic 


zone. 
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ahead  of  the  crack  occurs  on  slip  systems  other  than  those  of  Fig.  2b.  i.e..  on  systems  which 
cannot  undergo  sustained  plane  flow.  The  deformation  of  the  finite  element  mesh  around  the 
crack  tip  is  shown  in  big.  7a  with  the  displacements  magnified  by  a  factor  of  50.  The  most 
intense  deformation  is  in  those  elements  which  lie  along  the  discontinuity  between  sectors  A 
and  B.  The  increment  of  displacement,  scaled  up  for  visibility,  over  a  load  increment  is  shown 
in  Fig.  7b.  Four  distinct  sectors  exist  with  the  motion  in  each  sector  fairly  uniform.  This 
indicates  the  majority  of  deformation  is  caused  by  the  movement  of  nodes  near  the  predicted 
discontinuities. 

A  quasi-statically  propagating  crack  in  a  fee  crystal  is  simulated  using  a  node  release 
technique  the  details  of  w  hich  are  given  in  [6],  The  crack  propagates  to  the  right  through  *■) 
elements  under  constant  load  from  an  initial  stationary  position  to  the  center  of  section  A 
in  Fig.  5.  The  same  orientation  of  the  crack  with  respect  to  the  crystal  is  used  in  the 
propagating  crack  case  as  in  the  stationary  crack  case  above.  The  yield  surface  for  this 
orientation  is  repeated  in  Fig.  8a  for  reference.  As  discussed  in  the  previous  section  on 
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Fin.  .V,  Quasi-statically  growing  crack  in  ideally  plastic  crystal  of  fig.  2:  (a)  Yield  surface  and  stress  trajectory,  (hi 
Asymptotic  structure  of  near  tip  field  ( c )  finite  element  results  showing  plastic  activity,  as  measured  by  increment 
of  T.  during  one  finite  element  step  of  crack  growth.  Cirev  /one  corresponds  essentially  with  plastically  active  /one 
during  growth:  black  /one  has  more  concentrated  plastic  straining. 
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asymptotic  solutions,  neither  stress  nor  displacement  discontinuities  may  exist  and  if  a  sector 
is  plastic  it  must  border  an  elastic  sector.  The  assembly  of  sectors  for  this  orientation  [2]  is 
shown  in  Fig.  Hb.  For  the  upper  half  plane,  the  solution  again  consists  of  four  sectors  labeled 
A.  B.  C'.  and  D.  Sector  A  is  a  constant  stress  plastic  sector  and  it  is  separated  from  an  elastic 
sector  B  by  a  discontinuity  in  velocity  at  54.7°.  As  indicated  by  the  stress  trajectory  in  Fig.  Ha. 
sector  B  invokes  an  elastic  unloading  and  reloading  back  to  yield  by  the  time  a  second 
velocity  discontinuity  is  reached  at  125.3°.  Sector  C  is  another  elastic  unloading  and  reload¬ 
ing  sector.  Finally  ,  sector  L)  is  plastic. 

The  crack  leaves  behind  a  wake  of  plastic  deformation  which  makes  the  quantity 
T  not  as  illustrative  as  it  had  been  in  the  stationary  crack  case.  The  increment  in  T  (r„  (i ) 

(  =  AT  (r„  (/))  from  just  before  the  last  node  release  to  just  after  it.  is  shown  in  Fig.  He.  which 
eliminates  much  of  the  accumulated  plastic  wake  and  shows  only  the  regions  which  are 
plastically  active  during  growth.  Elements  which  sustained  AT  (r„  G )  between  0.001  to  1  are 
shaded  gray  and  those  with  increments  greater  than  1  are  in  shaded  black.  There  are  very 
striking  similarities  between  the  predicted  asymptotic  solution.  Fig.  Hb.  and  the  finite 
element  solution.  Fig.  He.  Sector  A  is  plastic  and  is  bounded  by  an  intense  band  of  plastic 
deformation  along  the  predicted  velocity  discontinuity  .  Sector  B  is  relatively  free  of  plastic 
deformation  as  predicted.  At  the  second  velocity  discontinuity  at  125.3°  where  the  pre¬ 
dicted  stress  trajectory  is  tangent  to  the  yield  surface  an  intense  band  of  plastic  defor¬ 
mation  is  encountered.  Subsequently,  an  elastic  sector  roughly  similar  in  angular  extent 
to  sector  (  is  seen  in  the  finite  element  solution.  Finally,  elements  along  the  crack 
laces  are  loading  plastically.  The  deformed  mesh  after  the  crack  growth  has  taken 
place  is  shown  m  Fig.  4a.  It  is  interesting  to  note  that  while  the  crack  was  stationary  (as  load 
was  increased)  the  solution  was  the  same  as  in  Fig.  7a.  as  shown  by  the  distinct  kink  in  the 
crack  profile  at  the  crack's  initial  position.  However,  once  crack  growth  began  and  a 
displacement  discontinuity  was  no  longer  allowed  the  crack  profile  was  smooth.  The  incre¬ 
ment  in  nodal  displacements,  again  scaled  up  for  visibility,  from  just  before  the  final  crack 
growth  increment  to  just  after  is  shown  in  Fig.  4b.  Two  rather  distinct  bands  of  discontinuity 
are  seen  close  to  the  crack  tip  corresponding  to  the  lines  of  the  velocity  discontinuities.  The 
smooth  change  of  direction  of  the  increments  directly  above  the  crack  clearly  indicates  that 
the  discontinuity  seen  at  40°  in  the  stationary  case  is  no  longer  present  in  the  growing  crack 
case. 

4.2.  (  enter  crack  panel  results 

A  square,  center  cracked  panel  under  plane  strain  conditions  has  been  modeled  from 
small-scale-yielding  to  general  yielding  in  the  finite  element  analyses  of  Asaro  [6].  Three 
sections  of  the  finite  element  mesh  are  shown  in  Fig.  10.  Section  A  is  imbedded  in  the  dense 
part  of  section  B  which  is  in  turn  imbedded  in  the  densest  part  of  section  C\  Section  C  is  then 
extended  outwards  in  a  similar  pattern  to  slightly  more  than  20  times  the  si/e  shown,  with 
the  outer  boundaries  of  successive  "'rings”  of  elements  forming  two  adjacent  sides  of  a 
succession  of  squares.  Symmetry  reduces  the  analysis  tojusi  iiuu  of  one  quarter  of  the  speci¬ 
men.  represented  by  section  C  as  extended  to  a  large  square  with  sides  of  dimension  b.  The 
total  length  of  the  crack  is  2a  with  the  ratio  of  a  b  equal  to  0.01.  The  ratio  of  the  si/e  of 
elements  in  section  A  to  the  si/e  of  the  elements  in  the  outermost  mesh  is  approximately 
2^10  ' .  The  crack  tip  is  initially  rounded  as  we  can  see  from  the  mesh  :n  section  A.  Stresses 


Crack  tip  fields  in  ductile  crystals 


315 


0.060 


-0.030  0.000  0.030 

xy(K/t0)2 


ha  v.  (.il  Deformed  mesh  (factor  of  100  amplification)  after  several  steps  of  crack  growth  at  constant  far-lieki 
stress  intensity  A, .  (h)  Line  lengths  proportional  to  displacement  increments  during  crack  growth  over  one  finite 
element  step  ol  crack  growth  Suggests  that  velocity  discontinuities  predicted  tn  asymptotic  analysis  extend  well 
into  the  plastic  /one. 


corresponding  to  uniform  tension  a'  perpendicular  to  the  crack  are  applied  to  the  outer 
boundary. 

Finite  rotations  of  both  material  and  the  crystal  lattice  are  taken  into  account  in  this 
work.  The  computational  procedures  used  in  the  finite  element  calculations  follow  the  initial 
work  of  [14]  and  [IS].  The  elements  used  are  the  same  cross  triangle  type  [17]  discussed  in 
conjunction  with  the  small  scale-yielding  analysis.  The  plasticity  is  modeled  by  a  visco-plastic 
formulation  with  the  slip  rate  on  each  system  as  given  by  (5)  with  m  =  0.005.  A  Taylor  type 
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tin.  10.  finite  element  mesh  lor  analysis  of  tensile  load  of  a  panel  containing  a  center  crack  with  initially  rounded 
tip.  Used  lor  crystal  models  with  double  slip,  as  in  insets  of  Figs.  11(a)  and  (b).  Strain  hardening  and  arbitrary 
displacement  gradients  are  included  in  the  analysis. 


of  hardening  of  slip  systems  is  included  with  g1'1  a  function  of  T  =  I,  |yu’|, 

X'7'  =  [I  4-  0.8  tanh  (11.1  D]  (8) 

where  gnjG  -  0.0026. 

The  crystal  is  idealized  by  a  planar  double  slip  model  for  two  orientations  of  the  crack  as 
shown  by  the  insets  of  Fig.  1  la  and  I  lb.  The  first  orientation  corresponds  to  that  in  Fie.  2b 
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Fig.  II.  Contours  of  equivalent  shear  strain  T  near  the  tip;  outermost  is  T  =  0.005,  next  is  T  =  0.01  and 
innermost  is  F  =  0.02:  (a)  For  35.3°  angle  between  slip  planes  and  tensile  direction;  (b)  For  54.7°  angle. 


for  a  fee  crystal,  i.e.,  a  (0  1  0)  crack  growing  in  the  [101]  direction,  when  we  neglect  the  slip 
plane  traces  shown  as  the  dashed  lines.  It  also  corresponds  to  a  bcc  crystal  with  the  crack 
on  the  ( 1  0  1 )  plane  growing  in  the  [0  1  0]  direction  when  we  similarly  neglect  certain  slip  plane 
traces  parallel  to  the  crack.  The  angle  between  the  .v,  axis  and  the  traces  of  the  slip  direction 
S,  denoted  by  (ft.  is  equal  to  35.3°  in  this  case.  The  second  orientation  has  (ft  —  54.7°.  It 
corresponds,  with  neglect  of  certain  systems  as  above,  to  a  crack  lying  on  the  (10  1)  plane 
growing  in  the  direction  [0  1  0]  in  a  fee  crystal,  or  to  a  crack  on  the  (0  1  0)  plane  growing  in 
the  [1  0  1]  direction  in  a  bcc  crystal.  See  [2]  for  fuller  discussion  of  yield  surfaces  in  these 
cases.  These  orientations  with  only  the  slip  systems  indicated  correspond  to  a  diamond 
shaped  yield  surface,  e.g..  as  in  Figs.  6a  and  8a  but  without  the  horizontal  cut-ofis.  An 
identical  diamond-shaped  surface  applies  for  the  two  cases  shown  in  Fig.  11a  and  b.  The  fee 
crack  orientation  coinciding  with  Fig.  1  la  is  sometimes  observed  in  fatigue  studies  on  ductile 
Cu  and  A1  crystals,  whereas  the  bcc  crack  orientation  noted  to  coincide  with  Fig.  1  lb  is  a 
common  cleavage  crack  orientation,  e.g..  in  Fe-3%  Si  crystals  [16].  The  term  g„  in  (8)  is  the 
critical  resolved  shear  stress  in  the  .v, ,  ,v:  plane  of  deformation  for  the  double  slip  model.  It 
corresponds  to  g„  =  (2/>/3)rll  for  the  fee  interrelations  and  to  g„  =  t„  for  bcc  [2]. 
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Contours  of  constant  plastic  strain  (as  measured  by  the  accumulated  sum  of  the  slips. 
T  =  I,  |  v111  |)are  shown  in  Figs.  1  la  and  1  lb  for  the  orientations  (f>  =  35.3°  and  <j>  —  54.7° 
respectively.  The  outermost  contour  is  for  T  =  0.005.  the  next  is  for  T  —  0.01  and  the 
innermost  for  f  =  0.02.  Away  from  the  immediate  near-tip  region,  the  plastic  deformation 
is  confined  to  two  discrete  bands  of  deformation  emanating  from  the  crack  tip.  For  a 
diamond  shaped  yield  surface,  the  vertical  (0  =  90°)  discontinuity  present  in  the  small-scale 
yielding  case  (associated  with  the  slip  system  represented  by  the  dashed  line  in  Fig.  2b)  does 
not  exist.  The  deformation  in  Fig.  11a  corresponds  to  shear  bands  with  the  line  of  dis¬ 
continuity  parallel  to  the  slip  direction  S.  However,  the  deformation  in  Fig.  lib  has  the 
discontinuities  aligned  with  the  slip  plane  normal  N  which  was  discussed  earlier  as  giving  rise 
to  kink  bands  and  lattice  rotation;  see  Fig.  3b  and  3d.  Not  until  we  are  within  two  to  three 
initial  root  radii  of  the  initial  blunted  crack  tip  does  the  straining  pattern  deviate  significantly 
from  that  predicted  from  the  small-scale  yielding  perfectly  sharp  crack  tip  results,  based  on 
the  "small  displacement  gradient"  formulation.  The  solutions  for  the  two  cases  in  Fig.  1 1  are 
identical  within  that  formulation,  and  the  slight  broadening  of  the  strain  contours  in 
Fig.  1  lb,  compared  to  those  in  Fig.  1  la.  reflects  geometrical  hardening  due  to  the  lattice 
rotation  in  the  Fig.  1  lb  case. 

The  maximum  extent  of  the  plastic  zone  r  varies  nearly  linearly  with  a(a'  lg„)2  x  A'f  for 
values  of  rnia  0.018,  such  that  rp  ^  0. 14  a(&'  jg„f.  (The  corresponding  numerical  factor 
is  0. 1 5  in  the  ideally  plastic  fee  solution  of  Hawk  ([6]  Fig.  6c).  and  both  results  are  close 
to  the  factor  0.16  based  on  the  approximate  ([19]  Eqn.  (37b))  model  for  yield  on  a  pair  of 
inclined  shear  planes  at  a  crack  tip,  when  the  angle  of  those  planes  is  taken  as  55°  with  the 
cracking  direction.)  For  larger  plastic  zone  sizes  in  the  Asaro  solution,  departures  from  such 
linearity  occur  indicating  a  departure  from  self-similar  growth  of  the  plastic  zone.  />..  a 
violation  of  small-scale  yielding,  presumably  explainable  over  some  range  of  a '  prior  to  full 
plasticity  in  terms  of  the  nonsingular,  crack-parallel  stress  ( =  -  &'  )  of  the  elastic  crack  tip 
field  (e.g..  [20]). 

Lattice  rotations  from  the  initial  orientations,  caused  by  plastic  deformation  at  the  crack 
tip.  are  shown  in  Figs.  12a  and  12b.  These  show  the  deformed  orientations  of  the  traces  of 
one  of  the  slip  plane  families  in  the  double-slip  model.  For  the  case  with  4>  =  35. 3C.  Fig.  12a 
indicates  that  although  quite  large  rotations  occur  at  the  crack  tip  itself,  very  little  reorien¬ 
tation  of  the  lattice  occurred  in  the  shear  band  regions.  This  is  compatible  with  flow  parallel 
to  slip  plane  traces  as  in  Fig.  3a.  In  contrast  to  this.  Fig.  12b  indicates,  for  the  case  where 
<(>  =  54.7°.  rather  discrete  bands  of  lattice  reorientation  occur  that  more-or-less  coincide 
with  the  bands  of  concentrated  straining  of  Fig.  1  lb.  By  sighting  along  the  slip  plane  traces 
in  this  latter  case,  one  may  observe  a  zone  of  kinking  shear,  as  in  Fig.  3b.  extending  out  from 
the  crack  tip.  This  is  further  illustrated  in  Fig.  12c.  also  for  the  (/>  =  54.7°  case,  in  which 
contours  of  constant  lattice  rotation  are  shown.  The  outermost  contour  corresponds  to  a  2° 
rotation,  and  it  protrudes  out  from  the  crack  tip  in  the  directions  of  the  bands  of  con¬ 
centrated  straining  in  Fig.  1  lb;  no  such  outward  protrusion  of  regions  of  significant  lattice 
rotation  occurs  in  the  <f>  —  35.3°  case.  By  comparing  Figs.  1 2a  and  b.  it  is  seen  that  although 
the  overall  plastic  zone  sizes  differ  little  in  the  two  cases,  the  patterns  of  large  deformation 
in  crack  tip  opening  are  quite  different  and  seem  to  be  dominated  by  the  geometrical 
hardening  and  softening  from  lattice  rotations.  The  results  also  show  the  development  of 
features  resembling  "high-angle  boundaries",  or  at  least  narrow  transition  zones  between 
regions  of  significantly  different  lattice  rotation,  in  the  near  tip  region. 
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Hr.  12.  Rotated  orientations  of  slip  plane  traces  (from  one  of  the  two  slip  systems  in  the  double-slip  model)  near 
the  crack  tip:  (a)  For  35.3°  angle  between  planes  and  tensile  direction;  shows  little  evidence  of  lattice  rotation 
except  near  blunting  tip  (b)  For  54.7°  angle;  shows  significant  lattice  rotation  as  in  the  kinking  shear  zone  of 
Fig.  3(b).  (c)  Also  for  54.7°  angle:  show  s  contours  of  constant  lattice  rotation  near  the  tip.  of  w  hich  the  outermost 
corresponds  to  2°  clockwise. 
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5.  Conclusion 

The  analytical  and  numerical  investigations  into  the  near-crack-tip  fields  of  ductile  crystals 
have  been  shown  to  be  consistent.  The  results  show  that  for  ideally  plastic  crystals  obeying 
the  Schmid  rule,  the  near-crack-tip  deformation  fields  are  characterized  by  discontinuities  in 
either  displacement  (stationary  crack)  or  velocity  (quasi-statically  growing  crack).  These 
discontinuities  are  either  parallel  or  normal  to  crystallographic  planes  on  which  slip  can  take 
place  (i.e.  by  the  motion  of  dislocations)  for  the  fee  and  bcc  orientations  considered.  The 
state  of  stress  is  constant  within  angular  sectors  bounded  by  these  discontinuities  for  material 
at  yield. 

The  asymptotic  methods  developed  [1,2]  provide  a  general  way  of  constructing  such  fields 
around  the  crack  once  the  yield  surface  has  been  determined  for  the  particular  crystal  and 
orientation  of  the  crack.  In  the  case  of  a  mode  HI  crack  [1.  3]  a  more  complete  analysis  is 
possible  of  stationary  cracks  (exact  full  analysis)  and  quasi-static  crack  growth  (asymptotic 
plus  approximate  analysis)  and  even  dynamic  crack  growth  (asymptotic).  The  general 
features  seen  in  the  mode  III  cases  carry  over  to  the  physically  more  interesting  mode  I  plane 
strain  asymptotic  analyses  [2], 

The  numerical  investigation  of  Hawk  [6]  by  detailed  finite  element  analyses  provided 
confirmation  that  the  features  seen  in  the  mode  I  asymptotic  analyses  (stationary  and 
quasi-static  growth)  are  valid  over  a  finite  size  region.  The  finite  element  analysis  of  Asaro 
[6]  of  a  center  crack  with  initially  rounded  tip  in  a  tension  panel  from  small-scale  yielding 
to  general  yielding,  including  hardening  and  full  account  of  arbitrary  displacement 
gradients,  shows  the  effects  of  lattice  rotation  and  the  limits  of  small-scale  yielding.  However, 
the  overall  feature  of  concentrated  deformation  along  the  predicted  discontinuities  is 
retained  throughout. 

Finally,  the  analyses  summarized  here  may  provide  some  insight  into  why  certain 
crystals  undergo  ductile  fracture  while  others  are  brittle.  One  possible  factor  is  the  struc¬ 
ture  of  the  dislocations  necessary  to  produce  the  predicted  slip  patterns  which  in  one 
case  involves  a  shear  band  parallel  to  the  active  slip  system  and  in  another  involves 
a  band  perpendicular  to  the  active  system  that  deforms  by  kinking  shear.  In  the  parallel 
shear  case  it  is  possible  for  these  dislocations  to  be  generated  at  the  crack  tip  and  swept 
out  along  the  slip  planes  while  in  the  kinking  case  it  is  necessary  to  have  present,  and  to 
activate,  internal  sources  to  generate  the  necessary  dislocation  dipoles.  The  finite  strain 
pattern  associated  with  large  ductile  opening  at  the  crack  tip  is  also  significantly  different  in 
the  two  cases. 
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In  Part  I  we  found  that  although  the  near  tip  fields  of  cracks  on  bimaterial  interfaces 
do  not  have  a  separable  form  of  the  HRR  type,  they  appear  to  be  nearly  separable 
in  an  annular  zone  within  the  plastic  zone.  Furthermore,  the  fields  bear  strong 
similarities  to  mixed  mode  HRR  fields  for  homogeneous  medium.  Based  on  our 
numerical  results,  we  have  been  able  to  identify  a  clear  mathematical  structure.  We 
found  that  the  small-scale  yielding  crack  tip  fields  are  members  of  a  family  param¬ 
eterized  by  a  near  tip  phase  angle  £,  and  that  the  fields  nearly  scale  with  the  value 
of  the  J-integral.  In  Part  II,  the  original  derivation  of  the  mathematical  structure 
of  the  small-scale  yielding  fields  is  elaborated  upon.  The  issue  of  crack  face  contact 
is  addressed  and  the  phenomenology  is  described  in  terms  of  the  phase  parameter 
£.  Crack  tip  plastic  deformation  results  in  an  open  crack  for  a  range  of  £  which  is 
nearly  symmetric  about  the  state  corresponding  to  pure  remote  tension.  Plane-strain 
plastic  zones  and  crack  tip  fields  for  the  complete  range  of  £  are  presented.  Over 
distances  comparable  to  the  size  of  the  dominant  plastic  zone,  the  stress  levels  that 
can  be  achieved  are  limited  by  the  yield  stress  of  the  weaker  (lower  yield  strength) 
material.  On  the  other  hand,  the  stresses  well  within  the  plastic  zone  are  governed 
by  the  strain-hardening  behavior  of  the  more  plastically  compliant  ( tower  strain¬ 
hardening)  material.  We  observe  that  the  extent  of  the  annular  zone  where  the  fields 
are  nearly  separable  (i.e.,  of  the  HRR  form)  is  dependent  on  the  remote  load 
combinations  and  the  material  combination.  When  the  tractions  on  the  interface 
are  predominantly  tensile,  there  are  no  indications  of  crack  face  contact  over  any 
length  scale  of  physical  relevance.  Instead,  the  crack  tip  opens  smoothly  and  crack 
tip  fields  as  well  as  the  crack  opening  displacement  are  scaled  by  the  I -integral.  The 
paper  concludes  with  a  discussion  on  the  range  of  load  combinations  which  could 
be  applied  to  two  fracture  test  specimen  geometries  to  obtain  valid  fracture  toughness 
data. 


1  Introduction 

In  Part  I  of  this  article  (Shih  and  Asaro,  1988),  numerical 
solutions  were  presented  for  the  elastic-plastic  fields  of  a  crack 
on  the  interface  between  a  nonlinear  power-law  hardening 
material  and  a  rigid  substrate.  Specifically,  the  problem  of  an 
infinite  crack  embedded  in  an  infinite  bimaterial  body  (see  Fig. 
1)  subject  to  combinations  of  remote  tension  and  shear  was 
analyzed  under  loading  conditions  that  caused  small-scale 
yielding.  The  calculations  were  performed  for  a  material  de¬ 
scribed  by  a  small  strain,  isotropic  J2  deformation  theory.  This 
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facilitated  making  the  connection  between  the  interface  crack 
solutions,  the  existing  framework  for  nonlinear  fracture  me¬ 
chanics,  and  specific  solutions  for  crack  tip  fields  in  homo¬ 
geneous  media  (e.g.,  Hutchinson,  1983).  In  one  sense  the  work 
complimented  analyses  such  as  that  of  Knowles  and  Sternberg 
(1983)  who  studied  the  behavior  of  an  interface  crack  between 
two  neo-Hookean  sheets.  In  these  analytic  asymptotic  solu¬ 
tions,  as  well  as  the  numerical  solutions  of  Part  I,  it  was  shown 
that  pathological  features  of  the  linear  elastic  solutions,  such 
as  the  near  tip  oscillations  in  stresses  and  displacements  are 
strongly  mitigated  by  nonlinear  kinematics  or  material  behav¬ 
ior.  The  linear  elastic  asymptotic  solutions,  which  provided 
the  point  of  departure  for  our  numerical  analyses,  were  dis¬ 
cussed  in  Part  I. 

The  numerical  analyses  described  in  Part  I  also  provided 
clarification  of  the  structure  of  the  small-scale  yielding  fields. 
In  particular,  it  was  found  that  the  interface  crack  small-scale 
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Fig.  1  Crack  on  blmaterial  Interface 


yielding  fields  are  members  of  a  family  parameterized  by  a 
near-tip  phase  angle ,  £,  defined  as  5  =  <t>  +  e  In  (QQ/c^i). 
Q  is  a  complex  stress  intensity  factor  defined  by  the  linear 
elastic  asymptotic  solution,  <j>  is  its  phase  angle,  L  is  a  char¬ 
acteristic  crack  dimension,  o0  is  the  yield  stress  of  the  weaker 
material,  and  e  is  a  elastic  bimaterial  constant.  We  found  that 
the  dependence  of  the  crack  tip  fields  on  radial  distance  is 
nearly  of  the  HRR  form  for  asymptotic  fields  in  homogeneous 
media  (Hutchinson,  1968;  Rice  and  Rosengren,  1968)  and  that 
the  small-scale  yielding  fields  nearly  scale  with  the  value  of  the 
7-integral  (Rice,  1968).  This  structure  is  reviewed  in  Section  2 
of  this  paper. 

In  the  present  paper  we  describe  the  small-scale  yielding 
fields  in  more  detail  and  also  demonstrate  how  the  pheno¬ 
menology  of  crack  face  contact  can  be  described  in  terms  of 
£.  In  the  absence  of  plastic  deformation,  the  range  of  £  for 
which  the  crack  is  open  is  skewed  about  the  state  corresponding 
to  pure  remote  tension.  In  contrast,  crack  tip  plastic  defor¬ 
mation  results  in  an  open  crack  range,  which  is  more  nearly 
symmetrical  about  the  state  corresponding  to  pure  remote  ten¬ 
sion.  We  make  contact  with  a  dimensional  analysis  by  Rice 
(1988)  and  a  subsequent  study  by  Zywicz  and  Parks  (1988), 
which  contain  results  which  are  consistent  with  our  general 
framework. 

The  plan  of  the  paper  is  as  follows.  In  the  next  section  the 
structure  of  the  small-scale  yielding  fields  is  described.  Nu¬ 
merical  methods  and  the  finite  element  boundary  problems  are 
described  in  Section  3.  Results  of  plane-strain  crack  tip  fields 
for  two  material  combinations  involving  an  elastic-plastic  solid 
bonded  onto  an  elastic  substrate  and,  secondly,  to  an  elastic- 
plastic  material  with  different  material  properties,  are  pre¬ 
sented  in  Section  4.  The  range  of  validity  of  linear  elasticity 
and  small-scale  yielding  solutions  as  they  pertain  to  two  frac¬ 
ture  test  specimens  are  discussed  in  Section  5. 

2  Structure  of  Fields  Under  Small-Scale  Yielding 

2.1  Linear  Elasticity  Solutions.  Solutions  to  specific  prob¬ 
lems  of  cracks  lying  along  bimaterial  interfaces  of  isotropic 
media  have  been  given  by  Cherepanov  (1962),  England  (1965), 
Erdogan  (i96j),  and  Rice  and  Gih  (1965).  More  recently,  Ting 
(1986)  has  presented  a  framework  for  determining  the  degree 
of  singularity  and  the  nature  of  the  asymptotic  fields  for  the 
general  interfacial  crack  between  two  elastic  anisotropic  ma¬ 
terials.  Park  and  Earmme  (1986),  Hutchinson,  Mear,  and  Rice 
(1987),  and  Suo  and  Hutchinson  (1988)  have  obtained  solutions 
for  several  elastic  interfacial  crack  problems,  and  Rice  (1988) 


has  reexamined  elastic  fracture  mechanics  concepts  for  inter¬ 
face  cracks. 

A  crack  of  total  length  L  lying  on  the  interface  between  two 
semi-infinite  slabs  of  isotropic  elastic  solids  with  differing  ma¬ 
terial  moduli  is  shown  in  Fig.  1 .  The  shear  moduli  and  Poisson’s 
ratios  are  n2,  and  v2,  and  r  and  0  are  polar  coordinates 
centered  at  the  crack  tip.  There  are  three  independent  non- 
dimensional  material  parameters  but,  as  has  been  shown  by 
Dundurs  (1969),  the  solution  for  this  class  of  problems  depends 
on  only  two  parameters.  In  plane  strain  they  are 


0i  =  U*,(l  -  f*)  -  M2(l  -  •'iMmiU  -  +  M<1  -  »i)l  (1) 

02=  [MlO-2»'2)-/*2(l-2l'l)]/[/il(l-J'2)  +  /i2(l-<'l)]-  (2) 


These  parameters  vanish  for  identical  materials  across  the  in¬ 
terface  and  change  sign  when  the  materials  are  interchanged. 

At  small  distances  from  the  crack  tip,  the  in-plane  stresses 
have  the  singular  form 


(3) 


where  i  =  \  -  1 ,  Q  is  the  complex  stress  intensity  factor, 
I  (r/Ly<  I  =  1,  and  5,/0;e)  is  the  universal  complex  dimen¬ 
sionless  angular  function  which  depends  on  the  bimaterial  con¬ 
stant  f  defined  by  (Williams,  1959)  as 


e  = 


(4) 


In  (3),  L  is  a  characteristic  dimension  of  the  crack  geometry. 
As  an  example,  L  is  identified  as  the  crack  length  of  the  ge¬ 
ometry  depicted  in  Fig.  1.  It  is  convenient  to  write  Q  as 

Q  =  e,  +  /Q2=IQIe'*  (5) 

where  IQ  I  and  <f>  are  the  amplitude  and  phase  of  Q.  Thus, 
tractions  on  the  bond  line  are 


t-  (<Jyy  +  iaxy)e=0- 


+  f  In (r/L)) 


(6) 


and  displacement  jumps  across  crack  faces  take  the  form 
Au  =  A  uy  +  iAux=  (Uy  +  iuJe^,,- (uy+iux)e= 

=  —  ^ - /QIa/f  pii-e  +  6  +  t  ln(r/L>)  (7) 

y/(l  +4f2)  V2x 

where  0  =  tan  ~'2e,  and  for  plane  strain  A  =  [(1  -  yj/ni 
+  (1  -  «’2)//i2]/(4cosh2xe). 

Different  normalizations  of  the  singular  crack  tip  fields  which 
result  in  stress  intensity  factors  that  differ  by  phase  angles  and 
scaling  constants  involving  e  have  appeared  in  the  literature 
on  interface  cracks.  Rice  (1988)  introduced  the  stress  intensity 
factor  K  which  is  related  to  Q  by 

K  =  QL  ie,  IKI  =  IQI,  KK  =  QQ  and  v>  =  4>-e  In  L  (8) 


where  ip  is  the  phase  of  K.  It  is  clear  from  (3)  that  K  or  QL~ " 
uniquely  characterizes  the  crack  tip  fields.  In  other  words,  K 
fully  accounts  for  the  effects  of  both  load  and  geometry  on 
the  crack  tip  field  whereas  Q  and  a  characteristic  length  L  are 
required  for  the  same  purpose.  The  merits  of  various  defini¬ 
tions  of  stress  intensity  factors  and  possible  approaches  for 
recording  and  using  fracture  toughness  data  have  been  dis¬ 
cussed  by  Rice  (1988). 

For  the  geometry  depicted  in  Fig.  1  and  stressed  by  remote 

traction  T(T=  ITIe**  =  <r"  +  /a"),  Q  =  (l  +  i2t)T\J nL/2  (Rice 
and  Sih,  1965).  The  connection  between  the  phase  of  Q  and 
T  is  <(>  =  \p  +  tan”  '2t.  Thus,  changes  in  crack  length  at  fixed 
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0(  IT  I  is  free  to  vary)  do  not  change  the  phase  of  Q.  For  later 
use,  we  note  that  <p  =  —  «lnl  +  0  +  tan  ' 1 2 e  for  the  geometry  of 
Fig.  1. 

A  relation  which  will  be  of  use  later  in  the  analysis  is  the 
energy  release  rate  for  the  crack  advancing  along  the  interface 
(Malyshev  and  Salganik,  1965 ;  Willis,  1971) 

8  =  AQQ  =  AKK  (9) 

where  A  is  the  parameter  defined  for  (7). 

The  linear  elasticity  solution  for  the  displacement  jumps  (7) 
predicts  that  overlapping  of  crack  faces  always  occurs,  though 
the  zone  of  overlapping  crack  faces  is  very  small  compared  to 
L  for  load  states  in  the  range  -45  deg  s  0  <  45  deg,  i.e., 
Qi  >  I  Qi  I  ■  To  redress  this  physically  objectionable  behavior, 
investigators  have  proposed  various  models  and  approaches. 
Comninou  (1977a,  b)  and  Comninou  and  Schmueser  (1978) 
reformulated  the  linear  elasticity  boundary  value  problem  to 
allow  a  zone  of  contact  to  develop  at  the  crack  tip.  Solutions 
for  a  range  of  remote  load  combinations  were  obtained. 
Achenbach  et  al.  (1979)  introduced  a  Dugdale-Barenblatt  strip 
yield  zone  at  the  crack  tip  which  eliminated  crack  face  over¬ 
lapping  as  well  as  stress  singularities  altogether.  Ortiz  and 
Blume  (1988)  have  obtained  an  inner  solution,  based  on  a  zone 
of  decohesion  (and  sliding)  at  the  crack  tip,  which  does  not 
lead  to  interpenetrating  crack  faces  in  the  region  of  dominance 
of  this  inner  solution.  Along  similar  lines  of  inquiry,  Needle- 
man  (1987)  has  implemented  a  decohesion  model  in  a  finite 
element  study  of  (rigid)  inclusion  debonding  which  takes  ac¬ 
count  of  finite  geometry  changes. 

We  take  the  view  that  while  the  elasticity  solutions  are  invalid 
on  the  scale  of  the  contact  zone,  they  still  provide  an  accurate 
description  of  the  fields  in  an  annular  region  surrounding  the 
contact  zone,  as  long  as  the  size  of  the  contact  zone  and  plastic 
zone  (if  any)  is  small  compared  to  the  zone  of  dominance  of 
the  elastic  fields  given  by  (3).  In  this  sense,  crack  face  contact 
on  a  size  scale  which  is  much  smaller  than  the  crack  length  (or 
the  relevant  crack  dimension)  can  be  treated  as  a  small-scale 
nonlinearity.  Along  this  line  of  argument,  Atkinson  (1982)  had 
used  (3)  as  the  outer  solution  and  matched  asymptotic  expan¬ 
sions  to  derive  an  analytic  contact  zone  solution  for  the 
Comninou  (1977a)  model.  We  will  return  to  the  issue  of  crack 
face  contact  following  the  discussion  on  the  structure  of  the 
small-scale  yielding  fields. 

2.2  Small-Scale  Yielding  Formulation.  In  the  small-scale 
yielding  formulation,  the  actual  crack  problem  is  replaced  by 
a  semi-infinite  crack  in  an  infinite  media  with  the  asymptotic 
boundary  condition  that  at  large  r  the  field  approaches  that 
given  by  (3).  The  original  derivation  of  the  results  to  be  dis¬ 
cussed  here  is  contained  in  an  earlier  publication  (Shih  and 
Asaro,  1988).  It  is  elaborated  upon  here  and  contact  is  made 
with  subsequent  studies  (and  alternative  approaches)  which 
lead  to  results  which  are  in  agreement  with  our  main  conclu¬ 
sions. 

Let  a0i  and  a02  be  the  yield  strengths  of  materials  1  and  2. 
It  is  convenient  to  designate  the  yield  strength  of  the  weaker 
solid  by  a0,  i.e.,  <r0  =  min  (<r0i>  <zo2)-  Within  the  boundary  layer 
formulation,  the  stresses  depend  on  Q L~'’  and  ct0,  and  on 
dimensionless  material  parameters,  e.g.,  aai/a0,  etc.  By 

dimensional  analysis, 

6'  phaseHz)1: 

dimensionless  material  parameters^  .  (10) 

Here,  f,  is  a  dimensionless  function  of  the  dimensionless  ar¬ 
guments  and  it  may  be  noted  that  QQ  has  dimensions  of  square 
of  stress  times  length.  Henceforth  we  will  assume  that  the 
implicit  dependence  of  f0  on  material  parameters  is  under¬ 
stood. 


Directing  our  attention  to  the  dependence  on  o,j  on  Q,  L 
and  distance  from  the  crack  tip  we  write 

°u  =  °o4  ,0,4>  +  e  In  (0  )  (11) 

where  <t>  +  t  In (r/L)  is  the  phase  of  Q ( r/L ) '* .  To  reveal  the 
structure  of fu,  it  is  advantageous  to  express  the  third  argument 
as  the  sum  of  a  constant  phase  angle  and  a  variable  phase 
angle  which  depends  on  distance  r.  To  this  end,  the  third 
argument  is  modified  by  incorporating  a  dependence  on  the 
first  argument,  viz., 

~  -(8) 

+  ln(^|)-  (12) 

We  should  point  out  that  the  phase  can  depend  on  In  r  in  a 
manner  which  is  more  complex  than  that  suggested  by  the  last 
term  in  (12).  In  Part  I  we  introduced  a  near-tip  phase  param¬ 
eter,  £,  defined  by 


1=0  +  6 


£  varies  linearly  with  0,  depends  weakly  on  QQ  and  L,  and 
does  not  involve  r.  Making  use  of  (13),  the  right-hand  side  of 
(12)  can  be  written  as 


£  +  phase 


MY1 

Vqq/ 


where  Z  is  a  dimensionless  combination  of  material  parameters 
yet  to  be  determined.  We  use  the  latter  result  to  restate  (11) 
as 

(i’#-pha”[(i)14  051 

For  the  small-scale  yielding  analysis,  itjs  convenient  to  phrase 
the  remote  load  in  terms  of  IKI  (or  KK/ojj)  and  its  phase 
Use  (8)  in  (13)  to  obtain 

/KK\ 

£  =  y>  +  eln^— J.  (16) 

The  characteristic  crack  dimension,  L ,  does  not  appear  ex¬ 
plicitly  in  the  equation  (16)  for  £,  but  has  been  absorbed  in 
the  phase  of  K.  Now,  the  general  result  in  (15)  has  the  alter¬ 
native  representation: 

phase((S)T0-  (i7) 

Since  £  is  the  phase  angle  of  a  complex  quantity,  fu  has  a 
periodicity  of  2*  with  respect  to  the  argument  £,  i.e., 

fij{ . :  £)=/,{ . ;f  +  rnir)  m  =  2,4,6.  (18) 

Furthermore,  due  to  the  linearity  of  the  equilibrium  and  strain- 
displacement  equations 

M . ■A)=-fJ< . iS  +  mrr)  1,3,5.  (19) 

Thus,  £  serves  as  the  phase  parameter  of  the  fields  in  the  small- 
scale  yielding  formulation  just  as  $  (or  0  -  rlnZ.)  is  the  phase 
angle  of  the  linear  elastic  singular  fields.  Using  a  different 
approach,  Zywicz  and  Parks  (1988)  introduced  a  phase  pa¬ 
rameter  f0  which  differs  inconsequentially  from  (16)  by  a  con¬ 
stant  involving  e. 

The  mathematical  structure  of  the  fields  expressed  by  (15), 
and  (17)— (19)  was  derived  without  regard  to  contact  between 
the  crack  faces,  in  fact,  displacement  jumps  across  the  crack 
faces  must  have  the  form, 
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-JrA 

i  S‘  \  KK  ’ 


(20) 


where  e0  is  a  reference  yield  strain  and  the  dimensionless  func¬ 
tion  gi  has  the  periodic  structure  expressed  by  (18)  and  (19). 
In  this  paper  our  interest  is  restricted  to  the  range  of  £  where 
the  crack  faces  are  not  in  contact.  The  definition  of  an  open 
crack  will  be  introduced  shortly  and  the  validity  of  the  proposed 
fields  will  be  established  in  Section  4  in  terms  of  the  range  of 
£.  For  the  present,  it  is  useful  to  examine  the  mathematical 
structure  of  the  fields  under  the  fictitious  assumption  that  crack 
faces  may  overlap. 

The  effective  stress  a,  (o^  =  is,/s,J/2,  where  s,y  is  the  deviatoric 
stress),  has  the  form 

o^o/P(£t,0,  Phase[(^)1:0'  (21) 


The  elastic-plastic  boundary  in  the  weaker  material  is  the  locus 
of  points  where  ae  equals  a0.  Substituting  these  values  in  (21) 
and  rearranging  leads  immediately  to  the  following  results  for 
the  plastic  zone  (Shih  and  Asaro,  1988): 


rPm  = 


KK 


Rm) . 


(22) 


Here,  R(8;  £)  is  a  dimensionless  angular  function  which  de¬ 
pends  on  £  and  on  dimensionless  material  parameters.  Since 
af  is  quadratic  in  the  stress  components  which  have  periodicity 
expressed  by  (18)  and  (19),  the  angular  function  R  has  a  pe¬ 


riodicity  of  t  with  respect  to  £, 

/?(«;£)  = /?(«;£  +  mu)  m=  1,2,3 .  (23) 

Rice  (1988)  has  derived  the  relation, 

KK 

rp  =—/?(<*>-«  ln(L/r„))  (24) 

al 


by  a  similar  dimensional  analysis.  In  (24),  rp  is  a  characteristic 
dimension  of  the  plastic  zone;  for  example,  rp  can  be  taken  to 
be  the  maximum  extent  of  the  plastic  zone.  Rice’s  implicit 
equation  for  r^can  be  made  explicit  by  replacing  rp  in  the 
argument  by  KK/<^,  whereupon  a  result  similar  to  (22)  is  re¬ 
covered. 

The  function  R(d;  £)  can  be  determined  directly  from  plots 
of  small-scale  yielding  plastic  zones.  This  has  been  carried  out 
for  a  range  of  phase  angles  £  in  Part  I  (Shih  and  Asaro,  1988) 
and  for  several  material  systems  in  this  paper.  Alternatively, 
an  estimate  of  R(8;  £)  can  be  obtained  by  approximating  /,  in 
(21)  by  using  the  linear  elasticity  fields  in  (3).  The  latter  pro¬ 
cedure  has  been  adopted  by  Zywicz  and  Parks  (1988)  to  es¬ 
timate  plastic  zone  sizes  and  shapes. 


2.3  Plausible  Form  of  Asymptotic  Fields.  The  size  of  the 
dominant  plastic  zone  (which  develops  in  the  weaker  materia!) 
is  controlled  by  the  yield  strength  of  the  weaker  material  a0 
(see  (22)).  Over  length  scales  which  are  comparable  to  the 
dominant  plastic  zone,  the  stress  levels  that  can  be  achieved 
in  both  materials  are  set  by  a0.  On  the  other  hand,  we  will 
argue  that  the  form  of  the  asymptotic  crack  tip  fields  is  gov¬ 
erned  by  the  hardening  characteristics  of  the  material  with  the 
least  hardening  capacity. 

We  assume  that  each  material  of  the  material  pair  is  char¬ 
acterized  by  a  Ramberg-Osgood  stress-strain  relation  where 
the  uniaxial  plastic  strain  is  related  to  the  stress  by  e^/e0  = 
a(<j/a0)n  where  f0  =  on/E.  The  plastic  properties  of  material 
I  (top  material)  are  designated  by  a,,  ct0i ,  and  n,  and  those 
for  material  2  are  designated  by  a2,  a01,  and  n2.  The  reference 
strains  are  defined  by  <01  =  <r0|/£|  and  ta2  =  °m/E2.  We 
designate  the  larger  of  n,  and  n2  by  n ,  i.e. ,  n  =  maxfrr,,  n2) 
To  avoid  ambiguity  in  the  discussions  to  follow,  material  1 


will  be  taken  to  be  the  lower  hardening  material  and  its  elastic 
and  plastic  properties  are  designated  by  v ,  E,  a,  aa,  t0,  and  n. 

Let  £■,,  and  Ea  designate  the  tangent  moduli  of  material  1 
and  2,  respectively.  As  r  —  0,  and  assuming  that  the  stresses 
are  singular,  Ea/Ea  —  0.  This  suggests  that,  as  r  —  0,  the 
material  system  behaves  increasingly  like  that  of  a  plastically 
deforming  material  which  is  bonded  to  a  rigid  substrate.  We 
will  now  argue  that  the  form  of  the  asymptotic  fields  is  gov¬ 
erned  by  the  strain-hardening  characteristics  of  material  1 . 

We  assume  that  the  response  of  the  plastically  deforming 
material  is  described  by  a  J2  deformation  theory  of  plasticity. 
Path  independence  of  the  ./-integral  (Rice,  1968)  is  easily  dem¬ 
onstrated  for  the  interface  crack  once  it  is  recognized  that  the 
contribution  to  the  /-integral  from  the  upper  interface,  x2  = 
0+,  is  negated  by  the  contribution  from  the  lower  interface, 
x2  =  0  •  Now,  take  a  circular  contour  centered  at  the  tip  and 
use  path  independence  to  shrink  the  contour  onto  the  crack 
tip.  Since  /  is  nonzero,  as  it  must  be  for  a  deformation  theory 
solid  (in  this  case  the  /-integral  is  precisely  the  energy  release 
rate),  we  conclude  that  the  strain  energy  density  W  ~  0(1  /r) 
as  r  —  0.  This  requires  the  product  of  stress  and  strain  to  have 
a  1/r  singularity.  If  it  is  further  assumed  that  the  fields  have 
a  product  dependence  on  the  first  argument  of/y  (17),  then 
singular  solutions,  if  they  exist,  would  have  the  form 

(KK\  ,/(n+ 

^r)  hn(8,  phase!/'-)  ;  £)  (25) 

where  r  =  /-/(KK/cto),  and  the  implicit  dependence  of  on 
material  parameters  is  understood. 

The  dimensionless  function  htJ  has  a  periodic  dependence 
(with  period  2ir)  on  £,  and  the  arguments  leading  to  (25)  require 
that  h,j  be  bounded.  Though  htJ  is  a  bounded  function,  it  can 
nevertheless  oscillate  rapiu'y  as  r  —  0.  Under  small-scale  yield¬ 
ing,  /  =  9  (Rice,  1968)  and  noting  the  relation  between  9  and 
KK  (9),  the  singular  fields  can  be  arranged  in  the  form 

(/  \  i/(n+ it 

- )  /i,y(0,  phase  |  r'c)  ;  £)  .  (26) 

aa0eor/ 

where  htJ  =  (a/AE)W{n*  and  A  has  been  defined  in  con¬ 
nection  with  (7).  It  may  be  noted  that  (26)  is  consistent  with 
the  product  form  of  the  linear  elastic  fields  (3)  involving  a 
bounded  hi2  which  depends  on  the  phase  |(r/Z,)'e ) . 

It  is  instructive  to  examine  the  dependence  of  h:j  on  r  for 
two  special  cases.  Suppose  that  both  materials  are  elastic  and 
e  ^  0.  At  fixed  8,  the  stresses  have  the  form  a^a  r~  hlh,j(r\ 
£,  n).  Comparing  this  form  with  (3),  we  conclude  that  the 
variation  of  A,y  with  r  is  bounded  but  rapidly  oscillatory  as  r 
—  0.  Now  suppose  that  one  material  is  rigid-perfectly  plastic. 
Then,  hIJ  is  again  bounded  and  the  term  multiplying  it, 
r  1/1 "+ 11  ( n  —  oo),  is  also  bounded. 

As  a  third  example,  we  take  materials  1  and  2  to  be  strain- 
hardening  materials  with  n,  >  n2.  Then  by  (26).  the  stresses 
above  the  bond  line  (8  =  0+)  are  of  the  form  a, y  - 
r~  11 h,j  ( r;  n,).  For  the  moment,  suppose  that  the  stresses 
below  the  bond  line  (8  =  O')  have  the  form  atJ  ~  r~  ,/("2  +  "/i,y 
(r;  n2).  Continuity  of  traction  across  the  bond  line  requires 
that  r  %/r;  n,)  =  r  'n"^')h2j{ry  n2)  =  0(r*).  For  n, 
*  n2,  the  equality  can  be  met  only  if  h2j  is  singular  as  r  —  0, 
which  contradicts  the  original  assumption  that  hl}  is  bounded. 
In  fact,  if  the  stresses  are  assumed  to  have  a  power  dependence 
on  rsuch  that  traction  is  continuous  across  the  bond  line,  i.e., 
<i,,  -  rxh,l  for  both  materials  1  and  2,  then  continuity  of  dis¬ 
placements  across  the  bond  ,;ne  cannot  be  satisfied.  Similarly, 
if  the  displacements  are  continuous  across  the  bond  line  and 
are  of  the  form  u,  ~  r^u,  for  materials  1  and  2,  then  continuity 
of  traction  cannot  be  satisfied.  However,  if  we  assume  that  as 
r  —  0  the  material  system  behaves  like  that  of  a  plastically- 
deforming  material  bonded  to  a  rigid  substrate,  then  (26)  proves 
to  be  an  admissible  form. 
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It  is  convenient  to  write  (26)  in  the  form 
/  J  \  l/*"+  *> 

°ij  =  o0  ( - )  hu(6,  r;  f  ,«)  (27) 

WoCo/-/ 

where  hy  is  a  bounded  function  with  respect  to  r  and 

?=*+e  ln  @  =*+e  ln  fif)  =*+‘ ln  (x^)- 

(28) 

The  form  of  the  fields  (27)  has  been  corroborated  by  the  nu¬ 
merical  solutions  of  Section  4.  These  full  field  solutions  also 
show  that  for  load  states  in  the  range  -  ir/6  <  £  <  ir/6,  hu 
is  a  slowly  varying  function  of  r  over  all  physically-relevant 
length  scales.  Thus,  the  near  tip  fields  nearly  scale  with 
Juln+l\  since  h,j  is  only  weakly  dependent  on  J  through  the 
phase  parameter  £.  Moreover,  these  near  tip  fields  of  an  in¬ 
terface  crack  display  strong  similarities  to  mixed  mode  HRR 
fields. 

When  the  plane  ahead  of  the  crack  is  stressed  by  tension 
and  shear,  the  HRR  singularity  (Hutchinson,  1968;  Rice  and 
Rosengren,  1968)  for  homogeneous  media  has  the  form  (Shih, 
1974) 

(J  \l/(n+l) 

- )  o,j{6\  \F,n)  .  (29) 

a<To<or/ 

Here,  M”  is  the  mixity  of  the  plastic  singular  fields  defined  by 
the  relative  magnitudes  of  and  a#  along  the  radial  line  9 
=  0  as  r  —  0.  The  angular  functions  atJ  for  several  values  of 
M 'p  have  been  tabulated  by  Symington  et  al.  (1988).  Under 
small-scale  yielding,  the  effects  of  load  and  geometry  on  the 
crack  tip  field  of  homogeneous  media  are  characterized  by  J 
and  A^.  The  structure  of  the  near  tip  fields  (27)  suggests  that 
J  and  £  are  the  characterization  parameters  of  the  near  tip 
fields  of  cracks  on  bimaterial  interfaces. 

2.4  Range  of  £  for  an  Open  Crack.  According  to  the  linear 
elasticity  solution,  overlapping  of  crack  faces  will  occur  at 
some  distance  for  all  values  of  <t>  if  «  /  0.  Nevertheless,  for 
a  range  of  </>  the  zone  of  interpenetration  is  confined  to  a 
distance  from  the  crack  tip  that  is  smaller  than  physically 
relevant  length  scales.  For  this  reason,  we  define  a  crack  to  be 
open  if 

Am2>0  for  rcsr<:L  (30) 

where  rc  is  small  compared  to  L.  The  size  of  the  contact  zone 
can  be  estimated  by  using  the  elasticity  solution  (7).  From  (7), 
the  opening  gap  between  the  crack  faces,  A w2  is  positive  over 
distances  r/L  which  satisfy 

-  +  e  \n(r/L)  -  0<  *72.  (31) 


(b) 

Fig.  3  (a)  Small-scale  yielding  formulation  and  (b)  finite  element  mesh 
of  upper  half  of  crack  tip  region  and  crack  tip  conventions 


Table  1 


t 

0(deg) 

r/L 

tl(deg) 

0.025 

2.86 

102 

6.6 

10  5 

9.89 

io-4 

13.2 

10  6 

19.8 

0.050 

5.71 

10  2 

13.2 

10  5 

19.8 

10  4 

26.4 

io-6 

39.6 

0.100 

11.3 

10  2 

26.4 

10  3 

39.6 

10  4 

52.8 

10  6 

79.2 

0.175 

19.3 

10  2 

46.2 

10  5 

69.3 

10- 

92.3 

10  6 

139.0 

For  the  case  e  =  0,  (31)  states  that  the  crack  is  open  as  long 
as  the  Mode  I  stress  intensity  factor  is  positive.  Let  U  =  -  e 
In (rc/L)  and  recall  that  0  =  tan~  '2e.  Then,  by  combining  (30) 
and  (3 1 ),  the  open  crack  range  in  terms  of  <t>  is 

—  ir/2  +  fi+  9<<ji  <  rr/2  +  0.  (32) 

Figure  2  depicts  the  open  crack  region  in  the  Q  plane.  Values 
of  fi  and  0  for  four  values  of  e  and  four  values  of  rc/L  are 
tabulated  in  Table  1.  Since  <t>  =  +  0  for  the  geometry  in 

Fig.  1,  the  range  of  phase  angles  of  the  remote  load  for  an 
open  right  crack  tip  is 

-  Tr/2  +  fl<^<x72.  T  (33) 

The  characteristic  length  in  the  small-scale  yielding  for¬ 
mulation  is  KK/co  (or  QQ/ol).  Therefore,  crack  face  contact 
under  small-scale  yielding  must  be  phrased_in  terms  of  the 
normalized  distance  f  (  =  r/(KK/ol  -  r/(QQ/ffo))-  The  form 
of  the  desired  result  can  be  obtained  with  the  help  of  the 
elasticity  result  (31).  Using  the  plastic  phase  angle  £,  defined 
by  (13),  we  rewrite  (31)  as 

-  x/2<£  +  e  ln  f—  0<  t/2.  (34) 

Under  small-scale  yielding,  the  crack  is  open  crack  if 

At/2&0  for  c,  <r<c2  (35) 
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Fig.  4  Elastic-plastic  material  bonded  onto  elastic  substrste— effective 
stress  contours  and  plastic  zones  of  six  load  states 


where  c,  «  I  and  c2  »  1.  Combining  (34)  and  (35)  and 
approximating  9  by  2t,  the  range  of  £  for  an  open  crack  is 

-  t/2  +  e(2  -  In  c,)  <  £  s  ir/2  +  «(2  -  In  c2).  (36) 

Suppose  In  c,  =  -2  and  In  c2  =  2  (c2  =  1/c,  =  7.5),  then 
(36)  reduces  to  -t/2  +  4e  <  £  <  t/2.  if  In  c,  =  -4  and 
In  c2  =  4  (c2  =  1/c,  =  55),  the  range  is  decreased  by 
4*.  -  t/2  +  6«  s  (  <  t/2  -  2«.  A  more  accurate  estimate 
of  the  open  crack  range  based  on  the  small-scale  yielding  so¬ 
lutions  is  discussed  in  Section  4. 

3  Problem  Description  and  Numerical  Procedures 

3.1  Boundary  Conditions  and  Mesh  Design.  Computations 
were  performed  in  Part  I  for  a  full  boundary  value  problem 
involving  an  interface  crack  lying  on  a  bimaterial  interface  of 
the  center-cracked-panel  geometry,  and  for  small-scale  yielding 
fields  produced  by  loading  the  crack  tip  with  a  remote  K  stress 
field  which  is  referred  to  as  the  boundary  layer  formulation. 


The  calculations  to  be  discussed  in  Section  4  are  all  of  the 
latter  type,  although  contact  will  be  made  with  the  full  bound¬ 
ary  value  problem  presented  in  Part  I.  Loading  is  prescribed 
by  imposing  the  stresses  of  the  asymptotic  linear  elastic  field, 
given  by  (3), 


°ij  -  Re 


(37) 


on  the  remote  boundary,  Rb  as  depicted  in  Fig.  3(a). 

The  finite  element  mesh  is  constructed  with  9-node  biquad¬ 
ratic  Lagrangian  elements.  An  arrangement  of  wedge-shaped 
elements  is  employed  at  the  crack  tip;  this  innermost  ring  of 
elements  has  dimensions  of  10" l5^.  The  upper  half  of  the 
crack  tip  mesh  and  the  spatial  coordinates  are  shown  in  Fig. 
3(h).  Accommodating  the  15  orders  of  magnitude  in  element 
size  while  still  preserving  the  aspect  ratio  of  the  element  is 
accomplished  by  exponential  scaling  of  element  size  in  the 
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Fig.  5  Elastic-plastic  material  bonded  onto  elastic  substrate— hydro¬ 
static  stress  contours  ot  six  load  states 


radial  direction.  Each  decade  of  radial  distance  is  spanned  by 
four  annular  strips  of  elements.  In  this  manner,  the  domain 
between  10' 15  Rh  and  Rb  is  spanned  by  60  strips  of  elements. 
Within  each  strip,  twelve  equally-sized  elements  span  the  in¬ 
terval  -  rr  <  6  <  0  (the  bottom  material),  and  another  twelve 
elements  span  the  interval  0  <  6  <  ir  (the  top  material).  The 
mesh  has  a  total  of  1464  elements  and  5980  nodes.  A  more 
complete  description  of  a  similar  finite  element  mesh  and  the 
element  arrangements  has  been  given  in  Part  I. 

3.2  Material  Laws.  The  deformable  media  are  described 
by  a  J2  deformation  theory  with  a  Ramberg-Osgood  stress- 
strain  behavior.  In  uniaxial  tension  the  material  deforms  ac¬ 
cording  to 

f/e0  =  a/a0  +  a(o/o0)n,  (38) 

where  an  and  «0  are  the  yield  stress  and  strain,  a  is  a  material 
constant  (taken  to  be  0.1),  and  n  is  the  power-law  strain- 
hardening  exponent.  Under  multi-axial  stress  states  the  strain 
is  given  by 


~rs"+ 


—  _  +  —  a 

IE  1  2 


Here  5,,  is  the  stress  deviator,  ae  =  (3sy  s,/2)l/2  is  the  effective 
stress,  and  v  and  E  are  the  isotropic  elastic  constants.  In  (38) 
the  connection  a0  =  Et0  was  used.  Numerical  solutions  for  n 
=  2  through  n  =  10  have  been  obtained;  in  the  present  paper, 
however,  only  results  for  n  =  5  and  selected  results  for  n  = 
10  are  presented. 


3.3  Numerical  Procedure.  The  numerical  procedure  em¬ 
ployed  in  this  study  is  discussed  in  Part  1.  Only  a  summary  of 
the  procedures  will  be  provided  in  this  section.  As  noted  in 
Section  3.2,  the  cracked  body  is  modeled  by  9-node  Lagrangian 
quadrilateral  elements.  To  avoid  numerical  difficulties  which 
arise  from  nearly  incompressible  plastic  deformation,  the  strain- 
displacement  matrix,  B,  is  formed  using  the  fl-bar  method 
(Hughes,  1980).  The  deviatoric  part  of  B,  Bdcv,  is  evaluated  at 
a  3  x  3  array  of  quadrature  points,  while  the  volumetric  part. 
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Fig.  6  Elastic-plastic  material  bonded  onto  elastic  substrate— normal 
displacement  jumps  of  six  load  states.  Crack  face  contact  (if  any)  is 
indicated. 


B™1,  is  evaluated  at  a  reduced  array  of  2  x  2  quadrature  points. 
The  desired  matrix  is  given  by 

B  =  Bd'v  +  Bvo1  (40) 

where  Bvo1  is  associated  with  the  3  x  3  quadrature  points  by 
a  bilinear  interpolation/extrapolation.  The  tangent  stiffness 
matrix  for  an  element  is  given  by  the  sum  of  the  inner  products 
at  each  Gauss  point,  viz., 

v 

k'=£(BDB)wp  (4!) 

p* 1 

where  D  is  the  matrix  of  material  moduli  evaluated  at  the  3 
x  3  quadrature  points  and  wp  is  the  appropriate  weight. 

The  solution  to  the  nonlinear  boundary  value  problem  is 
obtained  by  the  Newton-Raphson  method  which  is  second- 
order  convergent  if  a  close  initial  estimate  of  the  solution  is 
available.  Such  an  estimate  is  generated  by  parameter  tracking 
(Shih  and  Needleman,  1984). 

It  is  advantageous  to  work  with  area/domain  integrals  in 
finite  element  computations.  In  this  regard  it  may  be  noted 
that  with  the  help  of  a  weighting  function,  a  contour  integral 


(path-independent  or  not)  can  be  restated  as  an  area/domain 
integral  (e.g.,  Moran  and  Shih,  1987).  To  this  end  consider 
the  pa,lvindependent  ./-integral  defined  as  in  Eshelby  (1956), 
or  specifically  for  cracks  as  in  Rice  (1968),  by 


Here,  W  is  the  strain  energy  density,  oIJ  and  (du,/dXj)  are  'he 
stresses  and  displacement  gradients,  respectively,  F  is  a  -intour 
beginning  at  the  bottom  crack  face  and  ending  on  the  top  crack 
face,  and  n  is  the  outward  pointing  normal  to  T.  A  weighting 
function  <7,  (jr, ,  x2 )  is  assigned  a  value  of  unity  on  the  inner 
contour  T,  a  value  of  zero  on  the  outer  contour  C„  and  is  a 
smooth  (C”)  function  in  the  domain  bounded  by  T,  C„,  and 
the  crack  faces  C,  and  C  (these  contours  are  shown  in  Fig. 
4  of  Part  I).  The  integral  in  (42)  can  then  be  restated  as  an 
integral  over  a  closed  contour  C,  viz., 

•/=  jf  °,j  ^7  -Whi,  rnfl\dC .  (43) 

Here,  m  is  the  unit  normal  to  C  pointing  away  from  the  en- 
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Fig.  7  Variation  of  normalized  stresses  with  log  of  normalized  distance 
in  elastic-plastic  material  (9  =  7.S  deg)  and  in  elastic  substrate  (0  = 
-  7.5  deg)  for  three  load  states 


closed  area;  on  the  contour,  F,  m  points  opposite  to  n.  The 
derivation  assumes  that  the  crack  faces  are  traction-free.  Ap¬ 
plying  the  divergence  theorem  to  the  integral  in  (43),  and  using 
the  equations  of  equilibrium,  the  equivalent  domain  represen¬ 
tation  of  (42)  is 


J= 


l 


de)'  ^ 
—  dA 

dXj 


(44) 


where  A  is  the  area  enclosed  by  C.The  finite  element  form  of 
(44)  is  given  in  Part  I. 


4  Numerical  Solutions  for  Small-Scale  Yielding 

Full  field  numerical  solutions  under  plane-strain  assump¬ 
tions  have  been  obtained  for  a  range  of  bimatcrial  constants 
and  different  combinations  of  plastic  properties  for  materials 
1  and  2.  In  the  interest  of  length,  only  two  material  combi¬ 
nations  are  discussed  in  this  section. 

4.1  Elastic-Plastic  Material  on  Elastic  Substrate.  Material 
1  is  elastic-plastic  while  material  2  is  linear  elastic.  The  elastic 
properties  are  v,  =  v2  -  0.3  and  E2/E ,  =  2.5;  the  values  of 


the  Dundur’s  parameters  are  then  /3,  =  -0.429  and  02  = 
-0.122,  whereas  e  =  0.039.  Plastic  properties  of  the  top  ma¬ 
terial  are  characterized  by  yield  stress  a0,  strain-hardening  ex¬ 
ponent  n  which  is  equal  to  5,  and  a  which  has  a  value  of  0.1. 

Numerical  solutionsare  presented  for  seven  states  of  remote 
load  as  defined  by  KK/ooT  =  0.5,  and  the  elastic  phase  na- 
rameter  $  =  -  90  deg,  -  60  deg,  -  30  deg,  0  deg,  30  deg,  60 
deg,  90  deg.  The  small-scale  yielding  phase  parameter,  £,  de¬ 
fined  by  (13),  has  these  values:  -  1.60.  -  1.07,  -0.551, 
-  0.027, 0.496, 1 .02,  and  1 .54.  Under  small-scale  yielding  con¬ 
ditions, 

y=g  =  AQQ  =  AKK  (45) 

where  for  plane  strain,  A  =  [(1  -  v,)//*,  +  (1  -  i<2),  >  2J/(4 
cosh2  r e).  The  value  of  J  was  computed  using  (44)  and  a 
“plateau”  weighting  function.  The  J  values,  which  were  ex¬ 
tracted  from  annular  domains  with  mean  radii  that  differed 
by  up  to  fifteen  orders  of  magnitude  (10  l!),  agreed  to  better 
than  two  significant  digits  with  the  exact  result  given  by  (45). 
The  essentially  exact  path-independence  of  the  numerically 
determined  J  and  the  excellent  agreement  with  (45)  is  indicative 
of  the  quality  of  the  numerical  solutions. 
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Fig.  8  Variation  of  stresses  normalized  by  HRR  torm  with  log  of  nor- 
mallzed  distance  In  elastic-plastic  material  (f  =  7.5  deg)  and  in  elastic 
substrate  (8  =  -  7.5  deg)  for  three  load  states 


Plastic  zones  for  six  ofthe  loading  states,  £  =  -1.60,  -1.07, 
-0.551,  -0.027,0.496,  1.02  (KK /<&  =  0.5,0  =  -90  deg, 
-  60  deg,  -  30  deg,  0  deg,  30  deg,  60  deg)  are  shown  in  Fig. 
4.  Distances  are  normalized  by  KK/ct^  in  accordance  with  (22), 
and  the  stress  levels  indicated  in  the  plots  are  normalized  by 
<r0.  The  computed  plastic  zone  for  £  =  -1.60  is  identical  to 
that  determined  for  £  =  1.54  (as  demanded  by  (23))  and  the 
latter  is  therefore  not  shown.  The  material  above  the  bond  line 
is  elastic-plastic  and  thus  the  contour  labeled  by  ae  =  1 .0  is 
the  elastic-plastic  boundary;  the  contour  labeled  by  ae  =  1.5 
is  the  contour  for  ae  =  1.5ct0.  Maximum  extent  of  the  plastic 
zone,  Rp,  varies  between  0.15KK/cto  (for  £  =  0)  and  0.6KK/ 
ctq  (for  £  *  ±  1.5).  Contours  in  the  elastic  substrate  has  no 
physical  significance  other  than  being  indicative  of  the  effective 
stress  levels.  It  is  noteworthy  that  despite  the  significant  dif¬ 
ference  in  the  asymptotic  linear  elastic  fields  of  cracks  on 
bimaterial  interfaces  and  cracks  in  homogeneous  media,  the 
overall  size  and  shape  of  the  plastic  zones  are  similar  to  those 
found  for  the  homogeneous  medium  subject  to  combined  K, 
and  Kfi  loading  (see  Shih,  1974,  for  a  discussion  of  these  latter 
Fields). 

Hydrostatic  stress  contours  for  the  six  load  states  just  dis¬ 
cussed  are  shown  in  Fig.  5.  The  stress  levels  indicated  in  the 


plots  are  normalized  by  <?0.  Positive  values  correspond  to  hy¬ 
drostatic  tension  while  negative  values  correspond  to  hydro¬ 
static  compression.  The  stress  contours  for  £  =  -  1.60  are 
identical  to  those  for  £  =  1.54;  however,  the  sign  of  the  stress 
levels  have  to  be  changed  as  required  by  (19).  It  may  be  noted 
that  the  hydrostatic  stress  in  materials  1  and  2  are  very  de¬ 
pendent  on  the  load  phase  even  in  the  vicinity  of  £  =  0.  For 
£  =  -0.551,  the  top  material  endures  high  hydrostatic  tension. 
The  top  and  bottom  materials  are  subjected  to  high  hydrostatic 
tension  for  the  cases  £  =  -0.027,0.496. 

The  solutions  being  discussed  were  obtained  without  taking 
into  consideration  contact  of  the  crack  faces  if  this  should 
occur.  In  Section  2,  the  case  was  made  that  there  is  a  range 
of  load  states  where  crack  face  contact  does  not  occur  over 
any  physically-significant  length  scale.  Within  that  range  of 
load  combinations  we  have  explained  that  the  bimaterial  crack 
problem  can  be  analyzed  as  an  open  crack.  Based  on  the  elas¬ 
ticity  solution,  the  range  of  £  for  an  open  crack  (36)  was 
determined.  An  estimate  of  the  open  crack  range  based  on  the 
small-scale  yielding  solulion  is_now  given.  The  displacement 
jump,  Au2.  normalized  by  eoKK/oo,  is  plotted  against  logr,  (f 
s  r/(KK/aj-)),  in  Fig.  6.  Six  load^ states,  £  =  -  1.60,  -  1.43, 
-1.34,  -1.07,  -0.027,  1.28  (KK /<&  =  0.5,0  =  -  90  deg, 
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Fig.  9  Angular  variation  of  stresses  normalized  by  HRR  form  in  elastic 
substrate  (- 180  deg  <fsO  deg)  and  in  elastic-plastic  material  (0  deg 
<  9  <  180  deg)  for  three  load  states.  For  each  load  state,  the  radial 
distances  correspond  to  5  percent  and  0.5  percent  of  the  maximum 
plastic  zone  size. 


-  80.5  deg,  -  75  deg,  -  60  deg,  0  deg,  75  deg)  are  shown,  and 
the  zone  of  crack  face  contact  (if  any)  is  indicated. 

In  Section  2.4  a  practical  definition  for  an  open  crack  (35) 
was  introduced;  however,  values  of  c,  and  c2  were  not  given. 
Appropriate  values  for  c,  and  c2  can  be  determined  from  plots 
of  the  plastic ^one  in  Fig.  4.  The  maximum  extent  of  the  plastic 
zone,  /?p/(KK/ao),  varies  between  0.15  to  0.6,  depending  on 
the  phase  £.  Thus,  a  reasonable  value  for  c,  is  10“3  (this 
corresponds  to  a  distance  which  is  less  than  10“ 2  of  the  plastic 
zone  size),  while  c2  is  assigned  a  value  much  larger  than  unity. 
We  examine  the  displacement  jumps  in  Fig.  6.  The  range  of 
1 t  for  which  the  crack  is  opened  for  distances  f  >  10“ 3  is 
estimated  to  be 

-  1.3s£<  1.3.  (46) 

The  aforementioned  estimate  of  the  open  crack  range  pertains 
to  the  particular  material  combination  being  investigated. 
However,  crack  face  profiles  for  other  material  combinations 
have  shown  that  the  open  crack  range,  when  phrased  in  terms 
of  J,  is  not  strongly  dependent  on  material  combinations.  We 
have  completed  a  detailed  study  of  an  elastic-plastic  material 
bonded  onto  a  rigid  foundation  (e  *  0.1).  Materials  charac¬ 


terized  by  n  =  3,5,  and  10  were  considered.  The  open  crack 
range  differed  from  the  range  given  in  (46)  by  less  than  ten 
percent. 

We  have  previously  stated  that  for  the  present,  our  interest 
is  restricted  to  the  range  of  £  in  which  the  crack  faces  are  not 
in  contact.  Hence,  our  examination  of  detailed  behavior  of 
the  stresses  is  directed  to  load  states  where  the  crack  is  open. 
Stresses  near  the  bond  line  are  especially  of  interest  because 
of  the  likelihood  of  crack  propagation  along  the  interface.  The 
variation  of  stresses  with  radial  distance  from  the  crack  tip  in 
the  elastic-plastic  material  (along  8  =  7.5  deg)  and  in  the  elastic 
material  (along  8  =  -7.5  deg)  is  detailed.  In  Fig.  7  log  (<j0/ 
aa)  is  plotted  against  iogr  for  jhree  states  of  remote  loads,  £ 
=  -0.551,  -0.027,  0.496  (KK/o^T  =  0.5,  <6  =  -30  deg,  0 
deg,  30  deg).  Within  distances  that  can  be  resolved  by  our 
mesh,  10“ 12  <  f  <  103,  the  crack  faces  are  not  in  contact. 
The  stresses  in  material  1  are  shown  in  the  left  column  of  plots; 
two  distinct  slopes  can  be  seen,  a  higher  slope  in  the  elastic 
region  and  a  lower  slope  in  the  plastic  region  with  the  change 
in  slope  occurring  near  the  elastic-plastic  boundary.  Plots  on 
the  right  column  show  the  stresses  in  material  2— similar  trends 
are  observed  and  it  is  apparent  that  the  stress  levels  in  the 
elastic  material  is  limited  by  the  stress  levels  achievable  in  the 
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KK/crJ,  and  KK/oJ,,  respectively. 


elastic-plastic  material.  It  may  be  noted  that  the  behavior  of 
the  stresses  along  radial  distances  with  other  orientations  (e.g., 
6  =  ±67.7  deg,  ±97.5  deg)  are  similar  to  those  observed  in 
Fig.  7.  We  should  point  out  that  the  plotted  stresses  are  the 
actual  values  evaluated  at  the  Gauss  points — no  smoothing 
was  applied  to  the  data  points  in  our  plots  of  the  numerically 
determined  fields. 

To  examine  if  the  actual  fields  are  amenable  to  the  repre¬ 
sentation  given  by  (27),  the  numerically  determined  stresses 
are  normalized  as  follows: 

a  =  - -  a'/a° -  (47) 

1  {J/ao()t0r)t,inJ' 11 

The  plastic  properties  for  the  normalization  pertain  to  material 
1  (recall  that  material  2  is  elastic).  The  variation  of  the  nor¬ 
malized  stresses,  a0,  with  log/1  is  plotted  in  Fig^  8  for  the  three 
load  states,  £  =  -0.551,  -0.027,  0.496  (KK /<&  =  0.5,  <j> 
=  -  30  deg,  0  deg,  30  deg).  For  the  load  states,  £  =  -0.551 
and  -  0.027,  the  normalized  stresses,  a, ,,  well  within  the  plastic 
zone  do  not  seem  to  vary  with  distance,  i.e.,  the  actual  field 
support  the  proposed  form  (27)  and  the  observation  that  hit 


depends  very  weakly  on  distance.  Consider  the  rase  £  =  0.496 
which  displays  the  “strongest”  variation  with  distance  (bottom 
plots  in  Fig.  8).  J'rom  Fig.  4,  the  maximum  length  of  its  plastic 
zone  Rp  is  0.2KK/oo.  Over  the  distances  10  4  <  r/Rp  <  10'  ', 
which  include  all  physically-relevant  length  scales,  the  devia¬ 
tion  of  by  from  its  mean  value  is  less  than  eight  percent.  It 
can  also  be  seen  that  in  every  case  the  stress  levels  in  the  elastic 
substrate  are  set  by  the  stresses  in  the  elastic-plastic  material 
(as  they  must  by  continuity  of  traction). 

Angular  variation  of  6,,  at  fixed  radial  distances  equal  to  5 
and  0.5  percent  of  the  maximum  plastic  zone  size  are  shown 
in  Fig.  9  for  the  three  remote  load  states  discussed  previously. 
In  the  range  -  180  deg  <  6  <  0  deg  the  material  is  elastic, 
while  the  material  is  elastic-plastic  in  the  range  0  deg  <  8  < 
180  deg.  For  the  case,  £  =  -0.551,  the  hoop  stress  reaches 
its  maximum  in  material  1.  On  the  other  hand,  for  £  =  0.496, 
the  hoop  stress  attains  its  maximum  in  material  2.  In  all  cases, 
the  stresses  in  both  materials  are  comparable  in  magnitude  and 
the  normalized  values  are  of  order  unity.  More  importantly, 
we  point  to  the  strong  similarity  of  the  stresses  in  the  elastic- 
plastic  material  to  the  mixed-mode  HRR  fields  for  a  homo- 
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geneous  medium  characterized  by  the  plastic  properties  of  ma¬ 
terial  1.  In  fact,  the  radial  and  angular  stress  distributions  in 
material  1  for  £  =  -  0.027  agreed  closely  with  those  for  Mode 
I  HRR  fields.  Stress  distributions  for  £  =  -0.551  and  £  = 
0.496  are  in  agreement  with  the  mixed-mode  HRR  fields  for 
\F  =  -0.95  and  MP  =  0.78,  respectively  (Symington  et  al., 
1988).  We  defer  the  more  detailed  comparisons  to  Part  HI 
where  the  full  field  solutions  in  the  small  and  finite  strain 
regions  are  presented. 

4.2  Two  Bonded  Elastic-Plastic  Materials.  For  this  analysis 
both  materials  are  elastic-plastic.  The  elastic  properties  of  ma¬ 
terials  1  and  2  are  identical,  £2  =  £,  and  and  v2  =  0.3. 
Thus,  Dundur’s  parameters  (3,  =  02  =  0  and  e  also  vanishes. 
The  plastic  properties  of  material  1  are  characterized  by  yield 
stress  (to,  strain-hardening  exponent  n  of  5,  and  the  material 
constant  a  equal  to  0.1.  Material  2  yields  at  a02.  which  is  four 
times  higher  than  the  yield  stress  of  material  1.  The  strain¬ 


hardening  exponent  of  material  2,  n2,  is  10  and  a2  is  0.1.  The 
remote  elastic  fields  are  given  by  the  standard  K,  and  K„  fields 
and  £  =  <t>. 

Numerical  solutions  were  obtained  for  seven  values  of  the 
elastic  phase  parameter,  <t>  =  -90  deg,  -60  deg,  -30  deg, 
0  deg,  30  deg,  60  deg,  90  deg;  the  values  of  £  are  -  1.57, 
-1.05,  -0.524,  0.0,  0.524,  1.05,  and  1.57.  J  values  were 
extracted  from  numerical  Fields  using  domains  in  which  the 
mean  radii  differed  by  as  much  as  nine  orders  of  magnitude. 
These  values  agreed  to  better  than  two  significant  digits  with 
the  exact  result  (45). 

Since  the  yield  stress  of  material  2  is  four  times  higher  than 
the  yield  stress  of  material  1 ,  the  plastic  zone  in  the  (stronger) 
bottom  material  is  small  compared  to  the  plastic  zone  in  the 
top  material.  To  facilitate  the  comparison  of  plastic  zones  in 
both  materials,  distances  within  material  2  are  normalized  by 
KK/<4  while  distances  within  material  1  are  normalized  by 
KK/<4,.  Plastic  zones  for  six  load  states,  £  =  -  1.57,  -  1.05, 
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elastic-plastic  material  (9  =  -  7.5  deg)  for  three  load  states 


-0.524,0.0,0.524,  1 .05,  (0  =  -90  deg,  -60  deg,  -30  deg, 
0  deg,  30  deg,  60  deg)  are  shown  in  Fig.  10.  (The  computed 
plastic  zone  for  £  =  -  1.57  is  identical  to  that  for  £  =  1.57 
and  is  not  known.)  The  elastic-plastic  boundaries  in  materials 
1  and  2  are  defined  by  the  locii  of  ae  =  am  and  ae  =  a02, 
respectively.  Since  the  stress  levels  in  material  1  are  normalized 
by  <r0|,  while  those  in  material  2  are  normalized  by  a02,  contours 
labeled  by  a,  =  1.0  are  the  elastic-plastic  boundaries.  The 
normalized  plastic  zone  size,  Rpl/(KK/ooi),  for  the  weaker 
material  varies  between  0.15  and  0.6  and  these  values  are  rather 
similar  to  those  found  for  the  homogeneous  medium  subject 
to  combined  K,  and  Kn  loading  (Shih,  1974).  For  the  stronger 
material,  ^/(KK/c^)  varies  between  0.06  and  0.25.  Thus  the 
length  of  the  plastic  zone  in  the  stronger  material  is  about  only 
1/30  of  that  in  the  weaker  material,  i.e.,  the  relative  size  of 
the  plastic  zones  does  not  obey  the  ratio  (crnt/<702)2. 

Hydrostatic  stress  contours  for  the  six  load  states  just  dis¬ 
cussed  are  plotted  in  Fig.  11.  To  permit  a  direct  comparison 
of  the  stresses  in  both  materials,  the  stress  levels  indicated  in 
materials  1  and  2  are  normalized  by  the  same  yield  stress  <j0 ,, 
and  distances  in  both  materials  are  normalized  by  the  same 
length  KK/<toi.  The  stress  contours  for  £  =  -1.57  are  identical 


to  those  for  £  =  1.57;  however,  the  sign  of  the  stress  levels 
have  to  be  changed  in  accordance  with  (19).  It  is  apparent  that 
the  magnitude  of  the  hydrostatic  stress  in  both  materials  vary 
with  the  load  phase  for  moderate  departures  from  £  =  0 
(compare  plots  for  £  =  -0.524,  0,  0.524). 

Because  the  plastic  properties  of  materials  1  and  2  are  dif- 
ferent,  the  crack  tip  fields  do  not  have  a  separable  form  of 
the  HRR  type  (see  discussion  in  Section  2).  For  remote  load 
combinations  in  the  vicinity  of  £  =  <p  =  ±t/2,  (note  that  e 
=  0)  our  numerical  results  show  crack  face  contact  zones  which 
are  comparable  to  the  size  of  the  plastic  zone.  Subsequent 
discussions  are  therefore  restricted  to  those  cases  where  the 
crack  faces  are  clearly  not  in  contact.  _ 

Variation  of  the  normalized  stresses,  5y,  with  log  (r/(KK/ 
ctq))  along  6  =  ±7.5  deg,  for  £  =  -0.524,  0.0,  0.524  (</>  = 
-  30  deg,  0  deg,  30  deg)  is  plotted  in  Fig.  12.  The  stresses  in 
both  materials  are  normalized  by  the  yield  stress  and  the  hard¬ 
ening  exponent  for  material  1  (the  weaker  material).  It  can  be 
seen  that  the  normalized  stresses,  a0,  for  the  two  states,  £  = 
-0.524,  0,  are  only  slightly  dependent  on  radial  distance  in¬ 
dicating  that  a  form  of  the  type  (27)  is  a  suitable  representation 
over  the  radial  distances  considered  in  the  figure.  By  contrast, 
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Fig.  13  Angular  variation  of  normalized  stresses  In  (strong)  elastic- 
plastic  material  (-180  deg  sisf  deg)  and  In  (weak)  elastic-plastic 
material  (0  deg  <ts  180  deg)  for  three  load  states.  For  each  load  state, 
the  radial  distances  correspond  to  five  percent  of  the  maximum  plastic 
zone  size  in  the  weak  and  In  the  strong  elastic-plastic  materials,  re¬ 
spectively. 


stresses  normalized  by  the  yield  stress  and  hardening  exponent 
for  material  2  (the  lower-hardening  material)  display  stronger 
variation  with  rover  the  distances  considered  in  Fig.  12.  How¬ 
ever,  over  vanishingly  small  distances,  10“ 12  <  r  <  10‘4  (not 
shown  in  Fig.  12),  these  (latter)  normalized  stresses  are  nearly 
independent  of  f.  Such  behavior  is  consistent  with  our  obser¬ 
vation  that  as  r  —  0,  the  stresses  are  governed  by  the  strain¬ 
hardening  behavior  of  the  lower-hardening  material,  namely, 
material  2  for  the  material  system  being  considered.  These 
aspects  have  been  discussed  by  Shih  and  Asaro  (1989). 

The  angular  distribution  of  a,j  (normalization  based  on  the 
plastic  properties  of  the  n  =  5  material)  is  plotted  at  two  fixed 
radial  distances  for  each  of  the  load  states,  £  =  -0.524,  0.0, 
0.524  (<t>  =  -  30  deg,  0  deg,  30  deg).  The  first  distance  is 
chosen  to  equal  five  percent  of  the  maximum  extent  of  the 
plastic  zone  in  the  (weaker)  top  material,  and  the  other  distance 
is  chosen  to  equal  five  percent  of  the  maximum  extent  of  the 
plastic  zone  in  the  (stronger)  bottom  material.  These  angular 
distributions  are  shown  in  Fig.  13.  For  the  case  £  =  -0.524, 
the  hoop  stress  reaches  its  maximum  in  material  1 ;  it  is  at  its 
maximum  in  material  2  for  £  =  0.524.  It  is  noteworthy  that 
the  magnitude  of  da  of  both  materials  is  of  order  unity.  The 


similarity  between  the  angular  variation  of  the  stresses  with 
the  mixed-mode  HRR  fields  for  homogeneous  medium  is  dis¬ 
cussed  in  Part  III. 

5  Open  Crack  Range  for  Two  Fracture  Test  Geometries 

For  arbitrary  load  histories,  the  validity  of  the  solutions 
presented  in  this  paper  hinges,  in  part,  on  whether  or  not  the 
crack  faces  make  contact;  if  contact  develops  over  physically- 
significant  size  scales,  the  solutions  are  invalid.  We  have  shown 
that  load  states,  where  contact  between  the  crack  faces  is  ex¬ 
cluded  over  sensible  distances  near  and  away  from  the  crack 
tip,  can  be  expressed  by  a  range  of  £.  We  extend  the  discussion 
to  the  center-cracked-panel  and  confine  attention  to  the  case 
where  the  total  crack  length,  L,  is  small  compared  to  the  width 
of  the  panel.  Using  the  relation  <t>  =  +  9  in  the  linear 

elasticity  prediction  for  an  open  crack  (32),  we  observe  that 
the  right  crack  tip  is  open  if  ^  is  in  the  range, 

-  t/2  +  0s  \ts  t/2.  (48) 

Both  tips  are  open  if 

-ir/2  +  0<)('<xp/2-n.  (49) 
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The  open  crack  range  under  small-scale  yielding  is^ given  by 
(46).  Noting  the  connections  <t>  =  ^  +  ©  and  QQ  =  KK, 
approximating  0  by  2«,  and  using  the  relation  for  *  (13)  in 
(46),  we  obtain 

<(-!)-  «*» 

The  open  crack  range  given  by  (50)  takes  into  consideration 
the  effects  of  small-scale  plasticity.  The  maximum  radial  extent 
of  the  plastic  zone  is 


where  y  varies  between  6.5  (for  *  =  0)  and  1.5  (for  *  =  ±  ir/ 
2).  Using  (51)  in  (50)  leads  immediately  to 

2  +  ln  <i/-<1.3-e^2  +  ln  .  (52) 

Suppose  the  plastic  zone  in  the  center-cracked-pane!  is  con¬ 
fined  to  a  fraction  of  L,  say,  yRp/L  =0.1.  According  to  (52), 
the  crack  tips  are  open  if  -  1 .3  <  4/  <  1 .3.  On  the  other  hand, 
the  open  crack  range  as  determined  by  the  linear  elasticity 
solution  is  given  by  (48).  Suppose  that  rc/L  is  of  order  10" 3 
and  t  =  0.05,  then  fi  (=  -e  In (rc/L))  is  approximately  0.3. 
For  these  values  the  predicted  open  crack  range  by  linear  elas¬ 
ticity  is  -  1.27  <  ^  s  7r/2,  which  is  skewed  about  the  state 
corresponding  to  pure  remote  tension.  By  contrast,  crack  tip 
plasticity  results  in  an  open  crack  range  which  is  nearly  sym¬ 
metric  about  the  pure  tension  state.  Alternatively,  our  results 
indicate  that  states  of  loading  corresponding  to  nearly  remote 
shear  may  produce  contact  over  a  physically-significant  size 
scale  even  with  significant  amounts  of  plastic  yielding  at  the 
crack  tip. 

The  aforementioned  analysis  can  be  extended  to  other  spec¬ 
imen  geometries.  As  an  example,  consider  the  specimen  pro¬ 
posed  by  Charalambides  et  al.  (1988)  which  has  been  analyzed 
by  Suo  and  Hutchinson  (1988).  The  beam-like  specimen  con¬ 
sists  of  a  slender  strip  of  thickness  h  which  is  partially  bonded 
onto  a  thicker  substrate.  As  shown  by  the  latter  authors,  the 
specimen  design  is  such  that  for  a  broad  range  of  remote 
moment  and  axial  loads,  the  loading  mode  at  the  crack  tip 
corresponds  to  <6  =  ir/4  (note  that  the  strip  thickness  h  is  the 
characteristic  crack  dimension  designated  by  L  in  our  for¬ 
mulation).  To  assess  if  the  crack  tip  in  this  specimen  would 
remain  open  under  small-scale  yielding,  we  use  (46)  in  con¬ 
junction  with  (13)  and  (51)  to  get  the  following  condition  for 
an  open  crack: 

*y  Rn  yRn 

-  1.3-f  In  <0<1.3-e  In  .  (53) 

h  h 

For  the  specimen  under  consideration,  <f>  =  rr/4,  so  that  the 
condition  for  the  crack  to  remain  open  is 

-  1.3 - e  In  ^  <ir/4<1.3-e  In  ^  .  (54) 

h  h 

This  condition  suggests  that  for  nearly  all  physically  reasonable 
values  of  interface  fracture  toughness,  the  crack  tip  would  be 
open.  For  example,  if  t  -  0.04  as  in  the  case  presented  in 
Section  4.1,  (54)  would  be  violated  only  if  y Rp/h  ~  4  x  105! 


6  Concluding  Remarks 

(n  Part  I,  full  field  solutions  for  a  complete  boundary  value 
problem  were  presented  from  which  the  near  tip  fields  of  in¬ 
terface  cracks  were  extracted  and  characterized.  The  extent  of 
validity  of  the  small-scale  yielding  solutions  in  terms  of  the 
intensity  of  loading  and  resulting  plastic  zone  size  was  iden¬ 
tified  for  a  center-cracked-panel  geometry.  For  the  most  part 


our  results  were  obtained  numerically,  but  they  led  nonetheless 
to  a  clear  mathematical  structure.  In  particular,  it  was  found 
that  a  plastic  phase  angle,  *  defined  by  (28),  characterizes  the 
family  of  solutions  for  the  small-scale  yielding  fields.  In  the 
present  paper  the  structure  of  the  crack  tip  fields  under  con¬ 
ditions  of  small  scale  yielding  was  explored  in  more  detail  using 
a  boundary  layer  formulation. 

Our  numerical  results  show  that  the  size  of  the  dominant 
plastic  zone  is  controlled  by  the  yield  strength  of  the  weaker 
(lower  yield  strength)  material  and  that  over  length  scales  com¬ 
parable  to  the  dominant  plastic  zone,  the  stress  levels  in  both 
materials  are  set  by  a0,  the  yield  strength  of  the  weaker  material. 
Our  results  also  show  that  the  overall  size  and  shape  of  the 
dominant  plastic  zones  are  similar  to  those  found  for  the  ho¬ 
mogeneous  medium  subject  to  similar  combined  Kt  and  K„ 
loading  (Shih,  1974). 

We  have  observed  that  the  form  of  the  near  tip  fields  is 
governed  by  the  hardening  characteristics  of  the  material  with 
the  least  hardening  capacity.  With  the  stress-strain  behavior 
of  both  materials  being  given  by  piecewise  linear  elastic/power 
law  plastic  relations,  our  numerical  solutions  support  an  ob¬ 
servation  that  as  r  —  0,  the  asymptotic  stress  and  strain  fields 
in  the  lower  hardening  material  approach  those  for  the  same 
material  which  is  bonded  to  a  rigid  substrate.  Indeed,  the  near¬ 
tip  fields  for  the  elastic-plastic  material  bonded  to  an  elastic 
substrate  (see  Figs.  7,  8,  and  9)  are  very  similar  to  the  fields 
presented  in  Part  I  for  the  same  elastic-plastic  material  bonded 
to  a  rigid  substrate.  Furthermore,  the  present  small-scale  yield¬ 
ing  solutions  agreed  remarkably  well  with  the  center-cracked- 
panel,  small-scale  plasticity  solutions  for  the  identical  values 
of*. 

For  a  range  of  £,  the  actual  numerical  fields  appear  to  be 
nearly  separable  in  an  annular  region  well  within  the  plastic 
zone.  These  near-tip  fields  closely  resemble  mixed-mode  HRR 
fields  and  corroborate  a  proposed  representation  for  the  near¬ 
tip  fields  (27).  Our  results  also  show  that  the  stresses  in  the 
higher  hardening  material  are  limited  to  those  levels  that  can 
be  attained  in  the  lower  hardening  material.  These  aspects  and 
the  asymptotic  form  of  near  tip  stress  fields  are  elaborated 
upon  in  another  publication  (Shih  and  Asaro,  1989). 

For  selected  geometries  and  three  material  characterizations, 
(n  =  5, 10,  oo),  solutions  based  on  a  small  strain  J2  flow  theory 
have  recently  been  obtained.  The  plastic  zones  for  the  n  =  5 
and  n  =  10  materials  and  the  hydrostatic  stress  fields  are 
practically  identical  to  those  obtained  for  the  corresponding 
deformation  theory  solutions.  Near  the  interface,  we  did  find 
small  differences  between  the  deformation  and  flow  theory 
solutions.  Away  from  the  interface,  the  deformation  and  flow 
theory  solutions  are  hardly  distinguishable  from  one  another. 

Solutions  which  take  into  account  finite  geometry  changes 
associated  with  blunting  at  the  crack  tip  have  also  been  ob¬ 
tained.  In  this  case,  an  important  length  scale  is  the  crack  tip 
opening  displacement  6,.  Within  a  distance  of  several  5,’s  ahead 
of  the  crack  tip,  the  shear  stresses  near  the  interface  are  com¬ 
parable  to  those  determined  by  the  small  strain  solutions,  but 
are  much  smaller  than  the  values  predicted  by  the  linear  elas¬ 
ticity  fields;  in  the  finite  strain  region,  the  shear  stresses  rise 
slightly  due  to  the  asymmetrical  blunting  of  the  plastically 
deforming  material.  We  find  remarkable  agreement  between 
the  full  field  numerical  solutions  and  the  Mode  I  HRR  stresses 
over  the  distance  36,  <  r  <  66,  for  tension-dominated  load 
states,  -*•/  12  <  *  <  t/12.  For  remote  load  states  which 
include  moderate  amounts  of  shear,  there  is  good  agreement 
between  the  full  field  solutions  and  the  mixed-mode  HRR 
fields,  i.e.,  there  is  a  range  of  load  combinations  in  which  the 
crack  tip  fields  can  be  said  to  be  characterized  by  J  and  *  (or 
M v).  Furthermore,  for  load  states  in  which  o”  >  I  a"  I,  the 
crack  tip  opens  smoothly  and  the  amount  of  opening  scales 
nearly  linearly  with  the  7-integral.  These  aspects  are  discussed 
in  Part  ill  (Shih,  Asaro,  O’Dowd,  1989). 
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Abstract 

The  failure  of  bimaterial  interfaces  with  periodic 
arrays  o)  cracks  or  flaws  is  examined.  An  adhesive 
relation  is  prescribed  for  the  interface  which 
governs  the  normal  and  shear  stress  vs.  normal  and 
shear  displacement  response.  These  relations 
provide  what  amounts  to  a  constitutive  law  for  a 
thin  interface  layer  which  is  envisaged  to  undergo 
a  rupture  process  when  subjected  to  combined 
tensile  and  shear  forces.  On  the  basis  of  these 
models,  we  examine  the  effect  of  geometty  and 
interfacial  properties  on  the  f  ields  of  propagating 
cracks  and  flaws  as  well  as  on  the  ductility,  tough¬ 
ness  and  failure  modes  of  bimaterial  interfaces. 
Both  crystal  plasticity  and  J2  flow  theory  are  used 
to  model  the  plastic  deformation  of  the  materials 
on  either  side  of  the  interface.  In  both  descriptions, 
full  account  is  taken  of  finite  changes  in  geometty, 
large  plastic  strains  and  material  strain  rate  sensi¬ 
tivity.  Within  the  context  of  our  adhesive  models 
the  interface  properties  are  essentially  described  by 
three  parameters,  namely  the  interface  strength 
CTn, a*.  the  interface  normal  separation  be  beyond 
which  all  adhesion  is  lost  and  the  work  <pu  of  pure 
tensile  de-adhesion;  these  three  are  related  through 
relations  of  the  form  <jp„  =  AomMd^,  where  A 
depends  on  the  particular  f  unctional  form  of  the 
adhesive  relations.  Our  results,  performed  for 
metals  with  only  modest  rate  sensitivities,  demon¬ 
strate  the  pivotal  role  that  the  value  of  the  interface 
strength  has  in  establishing  the  mechanisms  of 
interface  failure  and,  thereby,  the  toughness  and 
ductility  of  the  microstructures  containing  the 
interface.  In  addition,  the  effects  of  crack  and  flaw 
size  on  the  stress-strain  fields,  and  failure  modes, 


that  develop  at  flawed  interfaces  are  studied  in  a 
series  of  calculations  involving  flaw  sizes  in  the 
range  15 <1  aj\  <  1500. 


1.  Introduction 

Our  previous  work  on  interface  failure  has 
been  concerned  with  detailed  analyses  of  the 
mechanics  of  interface  cracks,  and  in  particular 
with  the  analysis  of  the  stress  and  deformation 
fields  at  the  tips  of  cracks  lying  on  dissimilar 
material  interfaces.  Both  asymptotic  and  full  field 
analyses  have  been  carried  out,  and  from  these 
the  essential  structure  of  the  crack  tip  fields  has 
been  clarified.  These  analyses  ha,ve  assumed, 
however,  that  the  interfaces  possess  infinite 
strength  and  thus  were  intended  to  provide  a 
comprehensive  picture  of  stress  and  strain  fields 
at,  and  near,  imperfect  interfaces  at  stages  leading 
tip  to  crack  growth.  In  the  present  work,  analyses 
of  interface  debonding  due  to  the  propagation  of 
flaws  and  cracks  ( i.e .  king  flaws)  are  carried  out 
using  an  adhesive  model  for  the  interface  which  is 
also  used  in  the  work  of  Varias  et  al.  |1)  and 
McHugh  et  al.  [2|.  This  adhesive  model  is  moti¬ 
vated  by  the  atomistic  calculations  of  Rose  el  al. 
|3|.  The  adhesive  relations  in  effect  prescribe  a 
separate  constitutive  relation  for  the  interface 
which  can  be  used  to  model  the  effects  of  inter¬ 
face  structure  and  chemistry  on  bonding.  These 
adhesive  laws  are  phenomenological  in  nature 
and  are  meant  to  be  able  to  describe  the  influence 
of  interface  strength,  as  affected  by  structure  and 
chemistry,  on  the  overall  strength,  ductility  and 
toughness  of  microstructures. 
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In  what  follows,  solutions  for  the 
elastic-plastic  fields  and  growth  of  periodic  dis¬ 
tributions  of  interface  cracks  and  flaws  are  pre¬ 
sented.  These  defects  are  taken  to  be  completely 
non-bonded  segments  of  the  interface  whose 
dimensions  are  defined  relative  to  the  parameter 
dc  which  in  the  adhesive  model  measures  the 
separation  across  the  interface  beyond  which 
adhesion  is  completely  lost.  6C  (or  like  parameters 
that  enter  the  adhesive  force  laws)  introduces 
physical  length  scales  into  the  microstructural 
models  that  impart  different  behaviors  to  flaws  of 
different  sizes.  For  cases  where  defect  lengths  are 
large  compared  with  <5C,  i.e.  for  cracks,  the  elas¬ 
tic-plastic  fields  are  compared  with  those  for  sta¬ 
tionary  cracks  as  provided  by  Shih  and  Asaro 
[4-6].  When  the  defect  lengths  are  only  multiples 
of  (10- 100)^,  however,  the  fields  that  develop  at 
the  flaw  tips,  and  the  mechanics  of  debonding, 
can  be  quite  different  from  those  for  cracks. 
Details  for  both  types  of  defect  are  presented  in 
what  follows. 

For  the  specific  cases  presented  here,  the  inter¬ 
face  separates  an  elastic-plastic  material  and  a 
rigid  substrate;  the  crack  tip  field  analyses  of  Shih 
and  Asaro  [4-6]  have  shown,  however,  that  this 
provides  a  very  good  model  for  the  case  of  an 
elastic  substrate  as  well.  Remote  loading  is 
essentially  tensile  and  corresponds  to  remotely 
applied  displacements  normal  to  the  plane  of  the 
interface.  Cases  are  presented  where  fmite-strain- 
rate-dependent  J2  flow  theory  and  a  more  physi¬ 
cally  based  crystal  plasticity  theory  are  used  to 
model  the  behavior  of  the  plastically  deformable 
material.  In  both  cases,  full  account  is  taken  of 
finite  strain  and  geometry  change  effects.  All  the 
analyses  presented  here  are  two  dimensional  and 
correspond  to  states  of  plane  strain. 

Within  our  adhesive  laws,  the  interface  is 
characterized  by  three  parameters:  amax,  the 
maximum  pure  tensile  strength  of  the  interface; 
dc,  defined  above;  <p0,  the  work  expended  in 
completely  separating  the  interface.  These 
parameters  are  directly  related  via  (pQ  =  Aomi>ldc, 
where  the  value  of  the  constant  A  depends  on  the 
particular  functional  form  of  the  adhesive  laws. 
Our  calculations  indicate  that  amax  tends  to 
dominate  the  interface  behavior  and  thereby 
determines  the  mode  of  interface  failure  and  the 
resulting  toughness  or  ductility.  <5C  and  <p0  them¬ 
selves  play  important  roles  but  these  generally 
have  a  quantitative  significance  rather  than  the 
pivotal  qualitative  role  that  the  c/max  level  has. 


Examples  are  shown  of  the  separation  of  flawed 
interfaces,  where  <7max  controls  the  interface 
failure  mode,  which,  in  turn,  strongly  affects 
ductility;  dc  and  <pQ  have  significant  roles  in  sett¬ 
ing  the  values  of  strains  at  which  failure  occurs 
and  thus  the  precise  values  of  ductility.  Other 
examples  involving  long  cracks  are  presented, 
which  demonstrate  the  role  that  amax  plays  in 
establishing  the  value  of  toughness;  these  exam¬ 
ples  include  the  calculation  of  resistance  curves 
for  crack  growth  that  illustrate  the  strong  depend¬ 
ence  of  the  J  integral  vs.  crack  extension  relation 
on  the  value  of  amax.  It  should  be  noted,  however, 
that,  although  our  material  constitutive  relations 
are  strain  rate  dependent,  the  rate  sensitivities  are 
prescribed  to  be  rather  low,  consistent  with  what 
is  the  case  for  f.c.c.  metals  deformed  quasi-stati- 
cally  at  near  room  temperature.  The  significance 
of  this  is  that  for  rate-independent  materials  it  is 
well  known  that  continuum  fields  that  develop  at 
crack  tips  are  self-similar;  this,  in  turn,  means  that 
the  maximum  stress  levels  that  develop  ahead  of 
the  cracks  are  constant  and  are  fixed  in  magni¬ 
tude  by  the  material’s  strain  hardening  properties 
and  yield  strength.  It  is  then  expected  that  inter¬ 
face  failure  will  largely  depend  on  achieving  the 
stress  levels,  within  the  matrix,  necessary  to  over¬ 
come  the  interface’s  strength.  The  :nterface 
strength  is  thus  expected  to  dominate  the  failure 
behavior.  Further  details  on  the  relationship  of 
the  adhesive  properties  to  the  mechanics  of  inter¬ 
face  failure  are  described  in  Appendix  A,  which 
considers  a  simple  example  of  a  uniform  block, 
subjected  to  one-dimensional  tension,  debonding 
from  a  rigid  substrate.  This  analysis  provides  use¬ 
ful  insight  for  interpreting  our  results  on  the 
failure  of  flawed  interfaces  presented  in  Sections 
4  and  5. 

The  plan  of  the  paper  is  as  follows.  In  the  next 
section,  theoretical,  numerical  and  modeling 
issues  are  discussed.  In  Section  3  the  boundary 
value  problem  of  an  imperfect  interface  is 
described.  In  Section  4  the  results  for  growing 
cracks  and  flaws  on  intercrystalline  interfaces,  i.e. 
where  the  deformable  material  is  modeled  using 
crystal  plasticity,  are  described.  In  Section  5  the 
results  for  interface  failure  using  J2  flow  theory 
are  presented;  these  will  include  correlation  of 
the  elastic-plastic  fields  of  growing  cracks  and 
flaws  with  those  for  stationary  cracks  as  found  by 
Shih  and  Asaro  [4-6].  Section  6  contains  discus¬ 
sion  and  conclusions.  Finally,  in  Appendix  A  a 
solution  to  a  simple  boundary  value  problem  of  a 


perfectly  bonded  interface  is  presented.  1'his 
provides  insight  into  the  mechanics  of  the  de- 
adhesion  process. 

Standard  tensor  notation  is  used  throughout. 
Bold  italic  serif  symbols  are  used  to  denote 
vectors  and  bold  italic  sans  serif  symbols  to 
denote  tensors,  the  orders  of  which  will  be  clear 
in  context.  Products  are  indicated  with  dots  and 
products  containing  no  dots  are  dyadic  products. 
Italic  indices  range  from  one  to  the  number  of 
spatial  dimensions  and  repeated  italic  indices  are 
always  summed.  Inverses,  transposes  and  trans¬ 
posed  inverses  are  denoted  with  a  superscript 
-  1 .  T  and  — T  respectively,  and  superposed  dots 
indicate  differentiation  with  respect  to  time  /.  For 
instance. 

A  B  =  .\  ■/)’;  b,b  A  B  =  A  B, 
cil  —  c,c,bb,  cd=cd. 

HA  =  //„•  ,• 1  ..A  b  A  :H=  A., H,k„b,b, 

dll, 

8=  ■■‘b.b  B  c  =  />’ ,c. b 

di 

where  the  base  vectors  b.  are  cartesian  and  inde¬ 
pendent  of  time.  Greek  indices  are  slip  system 
identifiers  and  range  from  one  to  the  number  of 
available  slip  systems. 

2.  Theoretical  and  numerical  considerations 

2.  /.  Material  constitutin’  relations 

For  both  ./-  flow  theory  and  the  crystal  plas¬ 
ticity  theory  used  herein,  the  deformation  gradi¬ 
ent  is  taken  to  be  multiplicably  decomposable 
into  elastic  and  plastic  parts,  via 


where  x  is  the  current  position  of  a  material  par¬ 
ticle  and  X  is  its  position  in  a  reference  state.  Fp  is 
the  deformation  due  solely  to  the  plastic  flow  and 
represents  the  aggregate  effect  of,  for  example, 
dislocation  motion  in  a  ductile  crystalline  mate¬ 
rial,  which  leaves  the  lattice  structure  unchanged. 
F*  is  the  remaining  contribution  to  F  associated 
with  lattice  distortion  and  any  rigid  rotation  of  the 
lattice  and  the  material.  For  both  theories,  the 
kinematics  in  ( 1 )  lead  to  a  velocity  gradient  of  the 
form 

L  =  F  F  ' 

=  F*-[F*y'  +  F*-F?-(FPr'-[F*y'  (2) 


This  leads  to  plastic  and  elastic  parts  of  the  rates 
of  deformation  D  and  spin  W  given  by 

D  =  svni (L  ) 

=  Dp+  D* 

and 

l/V~skew(£  ) 

=  IV>’+  I N* 

The  specific  flow  relations  which  describe  crys¬ 
talline  plasticity  and  Jz  flow  theory  were  devel¬ 
oped.  in  detail,  by  Asaro  [7]  and  Barren  et  at.  [<Sj 
for  the  crystal  theory  used  here,  and  by  Moran  et 
al.  [9J  for  the  J:  flow  theory.  With  these  references 
in  mind,  only  a  brief  statement  of  the  material 
parameters  used  in  the  calculations  presented 
here  are  given  in  Sections  2.1.1  and  2. 1 .2. 

2. 1. 1.  The  crystal  plasticity  model 
In  crystal  plasticity,  plastic  flow  is  assumed  to 
occur  by  slip  along  discrete  sets  of  crystallo¬ 
graphic  slip  systems.  Figure  1  shows  the  slip 
systems  used  in  the  present  study  of  the  de¬ 
adhesion  problems:  a  are  the  time-independent 
cartesian  base  vectors  aligned  with  the  crystal 
lattice  in  the  reference  configuration.  These  slip 
systems  correspond  to  a  two-dimensional  close- 
packed  lattice.  Each  slip  system  is  defined  in  the 
reference  configuration  by  its  unit  crystallo¬ 
id 
I 


Fig.  1.  Slip  systems  of  a  two-dimensional  close-packed 
crystal.  The  lattice  base  vector  a,  is  parallel  to  the  interface 
in  the  reference  configuration. 


graphic  normal  m{a],  and  a  unit  crystallographic 
slip  direction  lying  in  the  slip  plane  *<«)■  With  yM 
being  the  rate  of  shear  on  that  system,  the  plastic 
part  of  the  velocity  gradient  in  (2)  is  given  by 
=  frp.  pp-i 

=  X  '  (3) 

a 


Here  the  sum  is  over  all  available  slip  systems  and 
is  explicitly  indicated. 

In  our  strain-rate-dependent  formulation,  the 
rate  of  shearing  on  each  slip  system  is  given  as  a 
function  of  the  current  state  of  stress  and  material 
state.  The  particular  form  used  here  is  the  adapta¬ 
tion  of  the  Schmid  rule,  described  by  Asaro  [7], 
i.e. 


a.a,  sgn(r„  i  i 


IS.  a,! 


(4) 


where  r(t!i  =s*la  -  e-m*, a)  is  the  current  value  of 
the  resolved  shear  stress.  The  Kirchhoff  stress  is 
given  by  r  =  (det  F)a  where  a  is  the  Cauchy 
stress  and  $*,„,  =  F*-s.;a),  m*la)  =  mla)-{F*)-].  gia. 
is  the  current  value  of  the  slip  system  hardness,  m 
is  the  rate  sensitivity  exponent  and  a(o)  >  0  is  a  ref¬ 
erence  shear  rate.  The  hardness  of  a  particular 
slip  system  g,a  is  obtained  by  the  path-dependent 
integration  of  the  evolution  equation 

l 

g„  =1  >'a  =  J  S  ly,o.,|  dt  (5) 

t>  0  ° 

where  hajt  is  a  matrix  of  (non-negative)  hardening 
moduli,  which  account  for  self-hardening  and 
latent  hardening  among  the  slip  systems,  and  ya  is 
the  accumulated  slip.  The  initial  condition  for  this 
evolution  is  specified  as  gu,(ya  =  0)  =  g„. 

The  strain-hardening  function  haji  is  given  by 


M; 


d}' 


(]otl 


(6a) 


where 


g1. *'., ) =  gu  +  hv ya  +  (g»  -  g„)  tanh  ya 


h0  -  h. 


\  g() 


(6b) 


Here  /i(l  and  /i„  are  the  initial  and  asymptotic  rate 
of  hardening  respectively.  g0  is  the  initial  slip 
system  hardness,  and  g „  can  be  interpreted  as  the 


saturation  hardness  when  hm  =  0.  qa$  is  a  matrix 
of  latent  hardening  ratios,  taken  to  conform  to 
Taylor’s  isotropic  hardening  model,  i.e.  qa$=\ 
for  all  a  and  /3.  The  material  parameters  were 
determined  by  matching  the  relation  (6b)  to 
actual  data  for  a  copper  single  crystal  [10]. 
The  values  are  fi0/g0  =  3.5,  goo/g0=  1-55  and 
h°o  I  go  =  1  ■  We  have  also  assumed  that  a(a}  =  1 0  “ 3 
s“ 1  and  m  =  0.005  which  are  fairly  representative 
values  for  the  range  of  quasi-static  strain  rates 
considered  in  this  paper. 

The  description  of  the  single  crystal’s  constitu¬ 
tive  behavior  is  completed  by  a  specification  of 
its  elastic  response,  which  is  expressed  in  terms 
of  the  Green  strain  of  the  lattice  given  by 
E*-?{(F*)r-F*  - /}  and  the  lattice-based  second 
Piola-Kirchhoff  stress  S*  =  (F*)~l- 
The  elastic  response  may  be  denoted  as  S*  = 
((0O/3£,y*),  where  S*  =  S*a,aJ,'  E*  =  E*alap 
and  <t>  =  <t>(£T)  is  the  Helmholtz  potential  energy 
of  the  lattice  per  unit  reference  volume.  In  rate 
form  the  elasticity  relations  are 

S*  =  K:E* 

K  kilkiajaJaka, 


0-<b 


K"u  0£*,,0£*, 


(7) 


In  this  study,  the  elastic  potential  is  taken  as 
< P  =  U0(£*, ,  +  E*22)2  +  £0(£*n2  +  £*22*") 

+  W£*,2  +  £*2l)2  (8) 

where  the  Lame  constants  have  the  values 
7.0  =  21 72.43g„  and  fi0  =  993.S3g(l. 

Combining  the  kinematical  relations  contained 
in  ( 1 )— (3),  we  get 


E*  =  (FP)~J ■  E-{FP)  ~ 1  —  ly\aAa)  (9) 

a 

where 

A,a)  =  sym{{F*)T-F*-s{a}m^,}  (10) 

and  E  is  the  rate  of  the  total  Green  strain  £  = 
i{FT-F-l). 

Making  use  of  S*  =  R-S  •  (£P)T,  where 
S  =  F~1-t-F~t  is  the  second  Piola-Kirchhoff 
stress  on  the  undeformed  configuration,  and  (3), 
we  get 


S*  =  fr.S-(P)T+2ly,a)B1„l 

a 


(11) 


69 


where 

B,a=sym(s,a)m{a)-S*)  (12) 

Combining  eqns.  (7),  (9)  and  ( 1 1)  we  obtain  the 
complete  constitutive  relation  in  terms  of  E 
and  S . 

S  =  I  :£-Z  (13) 

« 

where 

L  =  ( Furr V:/)~ 1KUm(  F^r '(. F„*)-'a,aiaran 

(14) 

and 

Xtl={Pr'-(K:A{at  +  2Bla))-(RrT  (15) 

2.1.2.  J2  flow  theory 

Although  phenomenological  in  nature,  J:  flow 
theory  can  be  formulated  along  similar  lines  to 
the  physical  crystal  theory  described  above.  For 
isotropic  polvcrystals,  it  is  assumed  that  the 
average  value  of  the  plastic  spin  W*p  among  the 
individual  grains  vanishes,  i.e.  W*p  =  0.  The  flow 
rule  for  D*p  =  sym|Fp-(Fp)~ '}  is  taken  to  be 

D*r=  ep  svm(/?*)  (16a) 

where 

C*  S*  C*  (16b) 

2  a 

f1'  is  the  effective  plastic  strain  rate  given  by  a 
flow  law  with  strain  hardening,  and  S*1  is  the 
deviatoric  part  of  S*  defined  by  S*'-S*~ 
'(S*:C*)(C*)_I.  The  effective  stress  6  is  de¬ 
fined  as 

d:  =  \{S*'  •  C*)\(S*'  ■  C*) 

In  the  present  work,  power  law  material  strain 
rate  sensitivity  and  strain  hardening  were 
assumed  in  the  form 


where  n  is  a  hardening  exponent,  o„  is  a  reference 
stress,  e„(=oJE)  and  et)  are  a  reference  strain 
and  strain  rate  respectively  and  m  is  the  rate- 
sensitivity  exponent.  In  this  study,  we  have  used 
n  —  5  and  m  =  0.005. 

The  elastic  response  is  formulated  in  terms  of  a 
hyperelastic  potential  defined  on  the  intermedr"e 
configuration  \V*(C*).  where  C*  is  the  elastic 
right  Cauchy-Green  tensor,  C*  =  (F*)1  •  F*.  The 
intermediate  configuration  is  thus  obtained  bv 
applying  the  mapping  F p  to  the  undeformed  state. 
Thus  S*  is  given  by 

a  iv* 

S*  =  2  — — -  =  S*(C*)  (19) 

dC*  v  ;  v 


and  the  elasticity  tensor  is  given  by  the  relation 


as*  3:W;* 

^  —  9  —  ^ _ 

“ac*  ac*ac* 


(20) 


The  elastic  potential  used  is  for  a  neo-hookean 
isotropic  elastic  material  extended  into  the  com¬ 
pressible  range  [11]  and  is  given  by 


VV*  =  i/(F*)  +  ~tr(C*)  (21) 

where  F*  -  det  F*. 

The  function  U  represents  the  extension  into 
the  compressible  range  and  is  taken  to  be 

U{F*)  =  Uo(log  F*)2  -  /r„  log  F*  (22) 

Combining ( 19),  (2 1)  and  (22),  we  obtain 

S*-At)  \og(F*)(C*)~]  +  fin\l-(C*)~  '}  (23) 

The  above  relation  reduces  to  Hooke's  law  with 
Lame  constants  A0  and  /x0.  when  the  limit  of  infin¬ 
itesimal  strains  is  considered.  In  this  study 
Afl  =  218.9 (70,  and  /t0=  1  12.8ct0,  corresponding 
to  Young's  modulus  E  of  300afl  and  Poisson's 
ratio  v  of  0.33. 


k  -i 

1  ll/flenl 
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6>H(ip) 

//(ep) 


(17) 


We  assume  that  FI{ip)  evolves  according  to  an 
isotropic  hardening  relation 


1  + 


G/ 


i/" 


(18) 


2.2.  The  adhesive  model 

2  2. 1.  General  considerations 
The  adhesive  model  used  in  the  present  work 
is  motivated  by  the  ab  initio  calculations  of  Rose 
et  al.  [3]  and  Ferrante  and  Smith  [  1 2).  These  have 
led  to  a  universal  adhesive  relation  for  bimetallic 
interfaces.  They  have  found  a  scaled  relation 
between  the  energy  and  interfacial  separation  of 
the  form 

E—  -A£(l  +  6*)  exp(  -<5*)  (24) 
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Fig.  2.  Diagram  of  adhesive  traction  vs.  interface  displace¬ 
ment  in  the  case  of  pure  tensile  straining. 


where  A £  is  the  equilibrium  adhesive  energy, 
<5*  s  (<J  -  6  )/ 1,  6  and  <5eq  are  the  actual  and  equi¬ 
librium  interface  separations  respectively,  and  /  is 
a  scaling  length,  defined  so  that  the  second  deriv¬ 
ative  of  the  scaled  energy-distance  curve  is  unity 
at  equilibrium.  When  (24)  is  differentiated  with 
respect  to  the  interface  separation  a  force- 
separation  (per  unit  area  of  interface)  relation  of 
the  form 


rr  u" 

an  =  £„  7-  exp 

Ocq 


(25) 


results.  £0  =  /?  A £//  is  an  initial  modulus  for  one¬ 
dimensional  tensile  straining  of  the  interface, 
P  =  <5cq//  where  /  is  of  the  order  of  the  electronic 
screening  length,  and  dn  =  6  -  <5eq. 

The  interface  stress  vs.  separation  relation 
is  shown  in  Fig.  2  where  several  parameters 
of  interest  are  illustrated.  For  example,  the 
maximum  stress  attained  from  eqn.  (25)  is 
amax  =  E{){pe)~ '.  Still  another  important  quantity 
obtained  from  (25)  is  the  work  of  interfacial  de¬ 
adhesion  <pn,  obtained  by  calculating  the  work  to 
separate  fully  the  interface,  i.e. 


or. 


Tn.i '  6, 

~p[e 


.4 


(26) 


The  numerology  is  such  that  amax  is  of  the 
order  of  a  tenth  of  the  elastic  modulii;  for  exam¬ 
ple,  if  there  is  no  initial  interface,  £0«(l  +  v)£ 
where  v  is  Poisson’s  ratio  and  £  is  Young’s  modu¬ 
lus.  In  the  more  general  case,  (p0  represents  the 
work  of  creating  two  surfaces  minus  that  which 
was  needed  to  create  the  initial  interface. 

To  explore  the  applicability  of  the  above  to  a 
fracture  theory  further,  it  is  instructive  to  apply  it 
to  one  of  the  many  adhesive  models  proposed  for 
crack  propagation  in  nominally  elastic  materials, 
as  has  been  done  by  Rice  and  Wang  [13].  Such 
models,  e.g.  those  proposed  by  Dugdale  [14], 
Barenblatt  [15],  or  Bilby  et  al.  [16],  consider 
cracks  in  homogeneous  elastic  media  which  open 
smoothly  at  their  tips  under  the  action  of  restrain¬ 
ing  adhesive  forces.  Generally  a  force-displace¬ 
ment  relation  of  the  form  shown  in  Fig.  2  is 
assumed.  The  elastic  crack  tip  singular  field  is 
replaced  by  the  adhesive  zone  whose  length  w  is 
estimated  to  be 
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(27) 


Then,  if  the  above  relations  for  omix  and  <p „ 
are  substituted  in  (27),  we  obtain  w’  =  (e:/2)(£/£u) 
xdcq~3.7(£/£0)c5cq.  If  we  now  apply  this  to  a 
homogeneous  medium  where  £n=(l+v)£~ 
1.33E  (with  v  =  3),  we  obtain  w,a»2.8<5t.q  which  is 
just  a  few  atomic  distances  or  so.  For  cases 
involving  interface  cracks,  £/£0  may  be  larger 
than  unity  as  suggested  by  the  work  of  Rice  and 
Wang  and  the  adhesive  zone  may  be  larger. 

The  specific  intent  here  is  to  explore  the 
implications  of  an  adhesive  theory  such  as  out¬ 
lined  above  to  the  fracture  of  interfaces  when  the 
materials  on  either  side  of  the  interface  undergo 
arbitrary  amounts  of  plastic  deformation.  In  other 
words,  we  seek  a  continuum  analog  to  the  result 
in  (25)  which  may  be  implemented  within  the 
framework  of  a  continuum  elastic-plastic  de¬ 
scription. 

For  this  purpose  we  recall  the  experiments  of 
Argon  and  coworkers  [17,  18]  and  of  Cialone  and 
Asaro  [19]  which  were  concerned  with  estimating 
the  critical  interfacial  stresses  required  to  sepa¬ 
rate  iron  (ferrite)-Fe3C  interfaces.  Conditions  of 
general  plasticity  prevailed  in  the  ductile  iron 
matrices  and  the  scale  of  observation  (by  scan- 


...  '.  jurt'  -Ji> 


ning  electron  microscopy)  was  typically  of  order 
0.05  /on  or  so.  When  displacements  of  the  ferrite 
from  the  carbide  surface  of  some  fraction  of  this 
resolution  was  observed,  debonding  was  said  to 
have  occurred.  Estimates  of  the  average  normal 
tensile  stresses,  i.e.  those  prevailing  over  the  size 
scales  of  fractions  of  a  micrometer,  required  to 
cause  this  were  3o0<  acri[^4a0,  where  o0  is  the 
matrix  yield  strength.  Within  the  context  of  atom¬ 
istic  studies,  dQq  is  in  the  range  (3-4) x  10_,°  m. 
For  -ontinuum  analysis  of  the  kind  considered  in 
this  investigation  an  appropriate  value  for  <5cq  is 
of  the  order  of  0.01  /rm,  i.e.  10-8  m.  To  give  the 
same  work  of  separation,  the  resulting  value  for 
amM  for  a  proposed  continuum  adhesion  model  is 
roughly  10  _  2  of  that  associated  with  the  atomistic 
model.  As  a  fraction  of  £0,  amax  is  3£0/100/Se, 
which  is  approximately  equal  to  £0/ 300.  This 
value  is,  in  fact,  a  reasonable  estimate  of  the  yield 
strength  of  the  ferrite  matrices. 

The  picture  then  in  our  continuum  analog  is 
that  of  a  thin  atomistic  layer,  of  order  100-200  A 
thick  that  undergoes  a  rupture  process  involving 
normal  displacements  of  1 00  A  or  so.  Since 
inelastic  processes  are  expected  to  occur  within 
this  layer,  <pQ  should  now  be  taken  several  times 
larger  than  the  work  required  to  create  free  sur¬ 
faces.  The  cleavage  energy  of  ferrite  has  been 
estimated  to  be  about  three  to  four  times  the 
interface  bonding  energy.  Therefore,  crmax  can  be 
adjusted  upwards  to  3cr0  <  amax  <  4cr0.  This  pheno¬ 
menological  version  of  the  adhesive  law  is  thus 
consistent  with  the  experimental  observations 
noted  above. 

We  adopt  the  above  model  to  describe  con¬ 
tinuum  separation  of  interfaces  and,  in  the  exam¬ 
ple  calculations  to  follow,  the  parameters  omax 
and  <5cq  take  on  values  in  the  ranges  quoted 
above.  Further  theoretical  studies  concerned  with 
the  atomistics  of  interface  separation  will  hope¬ 
fully  provide  further  insight  into  the  details  of 
such  laws  as  will  experiment  aimed  at  comparing 
the  predictions  of  interface  separation  obtained 
from  these  models  with  observation.  Two  pheno¬ 
menological  versions  of  continuum  separation 
relations  are  presented  next. 

2. 2. 2.  Continuum  adhesive  relations 

The  adhesive  Cauchy  tractions  are  governed 
by  relations  of  the  type 


_m,<5.) 


if  <5n  <  6C 
if  6„>dc 


G(<5„,  <5, 


if  <3„<6C 
if  6n  >  dc 


where  on  and  <$n  are  the  normal  traction  and  dis¬ 
placement,  and  a,  and  d,  are  the  tangential  trac¬ 
tion  and  displacement  defined  with  respect  to  the 
mean  surface  as  shown  in  Fig.  3.  serves  as  a 
critical  normal  separation  beyond  which  all 
adhesion  is  lost.  Thus  the  condition  6Jdc--\ 
provides  a  criterion  for  defining  the  current  flaw' 
tip  position  in  the  analyses  of  crack  growth  in  the 
subsequent  sections.  For  the  type  of  interface 
described  above,  dQ  can  be  taken  to  be  about  0.02 
gm.  The  function  F  attains  its  maximum  value, 
designated  by  omax,  at  a  value  of  dn  within  the 
interval  (0,  <5C).  Therefore  the  interface  strength 
omJX  is  the  maximum  normal  traction  that  the 
interface  can  sustain.  When  <5,  =  0,  the  variation 
in  on  v.v.  Sn  is  similar  to  that  shown  in  Fig.  2. 

We  now  use  (24)  to  motivate  a  continuum 
adhesive  relation  for  the  normal  tractions.  To  this 
end  we  restate  the  relation  (24)  in  the  form 

<p=  -A(l  +  J7<3„)  exp(  -  tjdn)  (30) 

The  relation  between  the  normal  stress  and 
normal  displacement  given  by  on  =  At]2dn 
x  exp(  -  t]dn)  is  derived  by  differentiating  <p  of 
(30)  with  respect  to  dn.  The  maximum  stress 
occurs  at  <5n  =  \jr]  and  it  has  the  value  omM  =  Arj/e . 
It  should  be  noted  that  the  work  of  separation  by 

reference  configuration 


mean  surface 


deformed  configuration 


mean  surface 


Fig.  3.  Interface  de-adhesion  model  for  a  deformable 

material  bonded  to  a  rigid  substrate. - ,  position  of 

material  interface  boundary  in  the  reference  configuration. 
The  de-adhesion  process  zone  has  initial  thickness  <50,  which 
is  comparable  with  dc. 
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tension  alone  which  has  been  defined  in  (26)  is 
given  by  <p0  =  eamaJrj.  In  terms  of  the  present 
parameters  an  equivalent  expression  for  the 
normal  stress  is 

o„  = - -  <5n  exp - <5„  (31) 

<Po  \  <Po  I 


this  paper  are  loaded  in  remote  tractions  which 
are  normal  to  the  interface.  Our  detailed  calcu¬ 
lations  show  that  the  amount  of  interface  shearing 
displacement  is  small  compared  with  the  normal 
displacement  so  that  the  actual  work  of  de¬ 
adhesion  is  nearly  equal  to  q>0  which  is  the  work 
of  interfacial  separation  by  tension  alone. 


A  phenomenological  relation  for  the  shear 
stresses  which  is  consistent  with  (29)  is  given  by 


e  or 


a, 


<Po 


(32) 


where  <5C  is  defined  as  the  normal  displacement  at 
which  the  work  of  separation  under  pure  tension 
is  equal  to  95%  of  <p0.  This  is  a  practical  defini¬ 
tion  of  a  characteristic  length  for  total  loss  of 
adhesion  which  is  consistent  with  the  polynomial 
adhesive  relation  to  be  defined  next.  It  should  be 
noted  that  (32)  has  been  scaled  so  that  the  initial 
shear  stiffness  is  identical  with  the  initial  tensile 
stiffness  in  (3  l). 

An  alternative  phenomenological  form  of  a 
continuum  adhesive  relation  can  be  obtained  by 
expressing  the  functions  F  and  G  as  polynomials 
in  dn  and  <5<.  One  version  of  such  a  relation  is 


27  <5 

<$,)  =  —  omax~ 

4  o„ 


G  v  d  n ,  <3 , ) 


27 


<5, 


4 


(33) 

(34) 


As  in  the  previous  relations  both  F  and  G  vanish 
for  values  of  r5n  greater  than  dc.  For  the  form  in 
(33),  F  is  maximum  at  dn  =  d J3  and  6,  =  0,  i.e. 


2.3.  Finite  element  implementation 
2.3.1.  Crystal  plasticity 

The  rate  constitutive  relation  (13)  is  imple¬ 
mented  via  the  one-step  explicit  rate  tangent 
method  of  Pierce  et  al.  [21].  The  time  discrete 
form  of  the  relations  that  result  from  the  algo¬ 
rithm  is 

S'  =  Lan\E'-Xan  (35) 

where  £,an  and  X',an  are  determined  explicitly 
from  the  state  at  time  t. 

To  enforce  equilibrium  at  the  end  of  the  time 
step  t  +A  t,  the  principle  of  virtual  work  is  written 
in  the  reference  configuration,  i.e. 

f  n'*A,:6FdV~  J  t'^'dudS  (36) 

v  s,  •■•s 

Here,  V  is  the  reference  volume,  S;  is  the  bound¬ 
ary  surface  between  the  adjoining  crystals  where 
adhesive  tractions  are  applied  and  5  is  the  regular 
boundary  over  which  tractions  are  imposed. 
t~N  n  is  the  nominal  traction  vector  {N  is  the 
outward  unit  normal  to  crystal  surface), 
n  =  F~ ’  •  r  =  S  FT  is  the  nominal  stress  and  du  is 
an  admissible  variational  displacement  field  that 
is  compatible  with  the  variational  deformation 
gradient  field  6F.  The  explicit  difference  relation 
n,+ti'  -  n'  +  ii'  A/ as  well  as 


n':dF=  S' 

:<5  E  dE =sym{(F')T-  bF) 

(37) 

n':dF  =  S' 

:dE+\S'-{F')7}:dF 

(38) 

is  substituted  into  (36)  to  obtain 


When  A,  =  0  throughout  the  process  of  de¬ 
adhesion.  the  work  of  de-adhesion  is  simply 
given  by 

Equations  (33)  and  (34),  when  interpreted  as  the 
nominal  tractions,  are  those  used  by  Needleman 
[20]  for  void  nucleation. 

The  flaw  and  crack  geometries  considered  in 


|  [<5£:/.,an:£'  +  {S'-(F')T}:(5/:]dK 

V 

xtan-— S'  :c5£dK+—  f  /,+A'-<5«dS 

(39) 

The  adhesive  relations  are  given  in  terms  of 
Cauchy  tractions.  In  terms  of  Cauchy  tractions 
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the  integral  over  St  is 

JVA'-5«dS=  J  o'+&'-6udS 


=  f  a'<3«dS  +  ^- 1 

f  o'-3Hd5]A/ 

J  3  t\J 

5/ 

'5/  ' 

—  /  j  +  1 2 

(40) 

where  o'  are  the  Cauchy  tractions  acting  on  the 
interface  at  time  t.  S'  in  the  current  configuration. 
It  should  be  noticed  that,  differentiating  the  last 
integral,  one  must  take  into  account  the  change  in 
the  interface  with  respect  to  time.  The  /,  term  is 
what  commonly  describes  the  tractions  acting  on 
the  crystal  boundary.  The  term  I2  introduces  new 
coefficients  of  the  nodal  velocities  for  the  inter¬ 
face  nodes,  due  to  the  dependence  of  the  adhe¬ 
sive  traction  rates  and  the  surface  variation  on  the 
nodal  velocities.  It  is  therefore  added  to  the 
global  stiffness  matrix,  which  is  derived  by  the 
usual  procedure.  The  time-discrete  variational 
equation  is  written  in  cartesian  index  form  on  the 
reference  orthonormal  base  vectors  e„  which  in 
this  problem  coincide  with  the  lattice  base  vectors 
a,.  Then  the  usual  finite  element  interpolation 
and  gradient  operators  are  introduced  to  obtain  a 
linear  system  of  equations.  After  applying  the 
appropriate  rate  boundary  conditions  and  solving 
the  resulting  linear  system  of  equations,  the  state 
at  time  t  is  updated  to  time  t  +  A t  in  a  manner  that 
is  consistent  with  the  method  of  Pierce  et  al.  The 
spatial  discretization  is  accomplished  by  using 
quadrilateral  elements  which  are  made  up  of  four 
cross-triangles. 

2.3.2.  J,  flow  theory 

The  starting  point  for  the  implementation 
adopted  here  is  (36).  Using  the  standard  finite 
element  discretization  we  get  a  system  of  non¬ 
linear  equations  for  the  deformation  mapping. 
This  system  of  equations  is  solved  using  a 
Newton-Raphson  iteration.  The  algorithm  is 
explicit  in  the  plastic  flow  direction  R*  and  hard¬ 
ening  moduli  H  which  remain  fixed  throughout 
the  increment  but  fully  implicit  in  the  increment 
of  effective  plastic  strain  Aep.  The  consistent 
tangent  used  in  the  Newton-Raphson  procedure 
follows  from  linearization  of  (36)  with  respect 
to  the  deformation  mapping.  The  linearized 


equation  is 

J  {S>:VuT(Au)}:dFdV+  j  (V/'^'Aw)<5«d.S 
y  s, 

=  /  t'*&,-dudS-~\ n,+&':6FdV  (41) 

S  +  S,  V 

Here  A u  are  the  incremental  displacements,  V0  is 
the  gradient  with  respect  to  undeformed  coor¬ 
dinates  and  V  represents  the  spatial  gradient.  25 
is  used  to  denote  the  gradient  of  n  with  respect  to 
F  at  t  +  At,  i.e. 

(42) 

The  first  term  on  the  left-hand  side  of  (41) 
gives  the  usual  contribution  to  the  stiffness 
matrix.  plays  the  role  of  an  effective 
elastic-plastic  tangent  on  the  reference  configura¬ 
tion.  Full  details  on  the  derivation  of  the  elas¬ 
tic-plastic  tangents  are  given  in  Moran  et  al.  [9]. 
The  second  term  on  the  left-hand  side  takes  into 
account  the  dependence  of  the  interface  tractions 
on  the  deformation  through  (33)  and  (34).  As  dis¬ 
cussed  earlier,  the  adhesive  relation  is  expressed 
in  terms  of  the  Cauchy  traction;  so  the  second 
integral  on  the  left-hand  side  of  (41)  is  re¬ 
written  as 

J  {\r*'-Au)-6udS 

s, 


where  a  is  the  adhesive  Cauchy  traction,  and 
|dx|/|dX|  is  the  ratio  of  the  current  arc  length  to 
the  original  length.  This  additional  stiffness  term, 
due  to  the  interface  tractions,  can  be  obtained  in 
closed  form  and  is  added  to  the  usual  stiffness 
matrix. 

In  this  implementation  the  body  is  discretized 
using  four  node  isoparametric  elements.  Numeri¬ 
cal  difficulties  may  be  encountered  with  this 
element  for  nearly  incompressible  deformations. 
To  overcome  these  numerical  difficulties,  we  use 
an  assumed  strain  method  which  is  essentially  a 
generalization  of  Hughes  [22],  B  method  to  finite 
deformations.  The  deformation  gradient  F  is 
decomposed  into  volumetric  and  deviatoric  parts 

p  z=z  jfy ol ^rdcv  J-  ^dcv  p*o\ 
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Fig.  4.  A  periodic  array  of  flaws  lying  on  an  interface 
between  a  deformable  body  and  a  rigid  substrate.  Because  of 
the  conditions  of  periodicity  and  symmetry,  only  the  shaded 
part  of  the  body  is  required  for  the  analysis. 

where 

pot  =  pi  If  picn  =  p-  1/3  p 

The  volumetric  part  F  of  the  deformation 
gradient  F  is  then  extrapolated  fr<  m  the  reduced 
quadrature  points.  The  resulting  volumetric 
deformation  field  incorporates  fewer  deforma¬ 
tion  modes  and  the  element  performs  well  in  the 
presence  of  near  incompressibility. 

3.  Model  boundary  value  problems 

Figure  4  depicts  the  geometry  of  the  periodic 
array  of  cracks  and  flaws.  Meshes,  spanning  up  to 
six  orders  of  magnitude  in  mesh  size  were  used  to 
resolve  the  elastic-plastic  fields  at  the  tips  of 
flaws  and  cracks  and  along  the  interface.  Sym¬ 
metry,  along  with  the  constraint  of  the  rigid  sub¬ 
strate,  leads  to  the  boundary  conditions,  w,  >  on 
X,  =0  and  X \~w,  corresponding  to  uniaxial 
straining  in  the  X2  direction.  On  the  top  bound¬ 
ary,  i.e.  Xz  =  Ln,  u 2  is  prescribed.  The  crack  or 
flaw  face  is  traction  free  whereas,  on  the  interface, 
tractions  are  prescribed  in  accordance  with  (31) 
and  (32)  or  (33)  and  (34). The  shear  tractions  are 
taken  equal  to  zero  on  the  sides  3^  =  0  and 
A',  =  w  and  on  the  top  boundary  X2  -  T0.  These 
geometry  and  boundary  conditions  lead  to  very 
high  triaxial  stress  states  as  the  deformation 
progresses  into  the  fully  yielded  regime.  Since 
deformation  and  stress  fields  near  notch  tips  are 
strongly  affected  by  triaxial  stress  states,  this 
aspect  of  the  boundary  value  problem  must  be 
kept  in  mind. 

T  hree  types  of  defect  are  considered-,  long 
flaws,  o0/<5c « 1 500,  which  behave  like  cracks; 
short  flaws  where  ajt 5C  =  15;  intermediate-size 
flaws  where  aJdQ  ~  1 50.  Most  of  our  analyses  are 
carried  out  for  these  three  ratios  and  some  addi¬ 
tional  results  are  presented  for  short  flaws  with 


a0/dc  <  50.  In  all  cases  examined,  the  ratio  aQ/b0 
of  the  initial  flaw  length  to  the  initial  ligament  is 
fixed  at  4,  and  the  block  height  is  1.25  bQ. 

To  provide  perspective  for  the  subsequent  dis¬ 
cussions,  we  note  that  for  a  long  flaw  there  exists 
an  annular  zone  near  the  flaw  tip  in  which  the 
stress  and  deformation  state  are  controlled  by  an 
asymptotic  crack  tip  field.  Specifically  the  effects 
of  load  and  crack  geometry  on  the  deformation 
within  this  annular  zone  are  transmitted  through 
the  stress  intensity  factor  K  when  non-linear 
effects  at  the  crack  tip  are  confined  to  a  fraction 
of  its  length.  Under  so-called  contained  or  laTge- 
scale  plasticity,  these  effects  are  transmitted  by 
the  J  integral.  In  this  context,  we  speak  of  an 
autonomous  K -  or  ./-dominated  crack  tip  zone 
where  the  asymptotic  fields  are  a  good  approxi¬ 
mation  to  the  actual  fields  over  distances  that 
extend  approximately  0.05flo-0.1ao.  In  other 
words,  a  flaw  can  be  treated  as  a  crack  if  the  frac¬ 
ture  process  zone  (or  dc  in  this  work)  is 
embedded  within  the  annular  K-  or  ./-dominated 
zone.  Specifically,  flaws  of  the  type  a0/dc«  1500 
turn  out  to  behave  like  cracks.  On  the  contrary, 
flaws  where  a0/(5c«15  cannot  be  treated  as 
cracks  since  dc  is  not  small  compared  with  the 
size  of  the  K-  or  /-dominated  zone.  For  flaws 
characterized  by  flo/<5c  =  150  the  behavior 
depends  upon  the  strength  of  the  interface  and  on 
the  sharpness  of  the  flaw’s  tip  during  its  growth. 

Calculations  were  performed  using  both  the 
exponential  and  the  polynomial  phenomeno¬ 
logical  adhesive  relations  described  above.  Most 
of  what  is  described  in  the  following,  however,  are 
results  obtained  using  the  polynomial  relation  in 
(33)  and  (34). 

Figure  5  is  a  representative  mesh  of  a  crack¬ 
like  flaw.  Two  types  of  notch  geometry  are 
employed:  an  elliptical  notch  tip  (Fig.  5(a))  and  a 
circular  notch  tip  (Fig.  5(b)).  A  representative 
mesh  of  short  flaws  is  shown  in  Fig.  6.  The  bulk 
of  our  results  are  for  flaws  having  elliptical  tips 
although  selected  results  for  flaws  with  rounded 
notch  tips  are  also  discussed. 


4.  Results  using  crystal  plasticity  theory 

4.1.  Long  flam 

The  ratio  a0/h0  of  the  initial  flaw  length  to  its 
opening  for  the  crack  geometry  (i.e.  long  flaw)  is 
500.  The  ratio  of  the  flaw  length  to  the  separation 
length  5C  is  approximately  1500.  Analyses  were 
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Fig.  5.  Finite  element  meshes  used  in  the  study  of  interfaces 
containing  long  or  intermediate-size  flaws  (<z,,/<5c  =  1 500  or 
150):  (a)  notch  with  elliptical  tip;  (b)  notch  with  circular  tip; 
(c)  remote  mesh. 


performed  for  interface  strength  amax  of  5g0,  6 g0 
and  7g(). 

For  all  three  levels  of  interfacial  strength,  the 
process  of  de-adhesion  begins  at  the  notch  tip 
and  continues  as  the  remote  load  increases;  the 
remote  load  is  indicated  by  the  level  of  remote 
tensile  stress  o ".  This  failure  mode  for  all  three 
cases  of  omax  is  shown  in  Figs.  7(a),  7(b)  and  7(c), 
where  the  interface  profiles  are  presented  at 
various  levels  of  remote  nominal  stress  defined  as 
the  total  applied  load  divided  by  the  reference 
cross-sectional  area.  The  interface  displacements 
increase  rapidly  as  the  distance  to  the  crack  tip 
decreases.  From  this  plot,  we  locate  the  current 
CTack  tip,  as  noted  previously,  by  the  position 
where  the  normal  opening  displacement  equals 
dc-  For  the  same  load  states,  the  distribution  of 
jhe  normal  stress  along  the  interface  is  shown  in 
Ftgs.  7(d),  7(e)  ancj  7(f)  The  stresses  increase 


Fig.  6.  Finite  clement  mesh  used  in  the  study  of  interfaces 
with  short  flaws  («„/($,.<  50):  (a)  remote  mesh;  lb)  near-tip 
mesh. 


rapidly  as  the  distance  to  the  crack  tip  decreases 
and  attain  a  maximum  at  a  distance  in  excess  of 
10<5C  ahead  of  the  current  notch  tip.  At  closer 
distances,  the  stresses  decrease  very  rapidly  and 
vanish  in  the  debonded  interface.  As  the  crack 
advances,  the  zone  of  peak  stress  translates  with 
the  crack  tip;  however,  the  high  stress  gradient 
associated  with  the  initial  process  of  crack  growth 
begins  to  weaken.  This  behavior  has  a  ready 
explanation.  A  weaker  stress  gradient  is  a  well- 
known  consequence  of  crack  giowth  [23]  and  the 
process  of  crack  growth  can  only  be  maintained 
by  the  attainment  of  a  peak  stress  of  omax  as 
required  by  the  adhesive  relation.  It  is  apparent 
from  our  results  in  Fig.  7  that  higher  stress  gradi¬ 
ents  develop  along  the  stronger  interface.  The 
effect  of  omax  on  stress  gradient  is  clearly  shown  in 
Fig.  8(a),  where  the  stress  distribution  for  ap¬ 
proximately  the  same  amount  of  crack  growth 
and  normalized  by  the  respective  amax  is  plotted. 
From  this  and  preceding  plots,  it  is  apparent  that 
a  high  normal  stress  extends  over  a  larger  zone  as 
the  strength  of  the  interface  decreases.  The 
hydrostatic  stress  along  the  interface  behaves  in  a 
similar  manner  and  it  is  not  plotted.  The 
deformed  notch  tip  and  the  interface  profile  are 
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rig.  7.  Interface  profiles  and  t22  Kirchhoff  stress  distributions  alone  the  interface  for  an/dc  ~  1 500.  lx,  y)  are  the  cartesian  coordi¬ 
nates  measured  with  respect  to  the  reference  coordinate  system  depicted  in  Fig.  4. 


shown  in  Fig.  8(b);  it  can  be  seen  that  the  inter¬ 
face  displacement  gradient  also  decreases  with 
decreasing  interface  strength. 

The  accumulated  slip  ya  is  also  strongly 
dependent  on  the  interface  strength.  This  can  be 


seen  by  the  contour  of  the  accumulated  slip  at 
y,  =  0.001  shown  in  Fig.  9(a)  for  the  three  differ¬ 
ent  amax  values  and  at  the  same  amount  of  crack 
growth  (Aa  =  20<5c).  The  region  enclosed  by  the 
contour  ( ya  =  0.00 1 )  is  effectively  the  plastic  zone; 
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Fig.  8.  (a)  Distributions  of  r;,  Kirchhoff  stresses,  along  the 
interface,  normalized  by  the  strength  of  the  interface;  (b) 
interface  profiles.  For  each  case  considered  the  amount  of 
crack  growth  is  approximately  20<5C. 


the  plot  demonstrates  the  substantial  increase  in 
the  plastic  zone  with  increase  in  interfacial 
strength.  We  have  observed  that  the  accumulated 
slip  (ahead  of  the  notch  tip)  during  the  process  of 
notch  blunting  and  the  initial  stages  of  crack 
growth  is  in  excess  of  0.1.  After  an  amount  of 
growth  of  about  5<3C  the  accumulated  slip  at  a 
comparable  distance  ahead  of  the  tip  decreases  to 
less  than  a  quarter  of  the  value  at  initiation  of 
growth  and  thereafter  remains  nearly  constant 
with  additional  crack  growth.  The  magnitude  of 
the  plastic  strains  is  strongly  affected  by  amaj(; 
for  example,  ya  at  the  initiation  of  growth  is  about 
0-4  for  omax  =  7g0  while  it  is  less  than  0.2  for 
5g0.  Hydrostatic  stress  contours  of 
Th  =  4.0g0  for  two  interfacial  strengths  om3x  =  6g0 
atd  7g0  are  shown  in  Fig.  9(b)  for  the  same 
^ount  of  crack  growth  as  in  the  previous  plots. 
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Fig.  9.  (a)  Accumulated  slip  contours  of  y,  -  0.00 1  for  three 
values  of  interface  strength;  (b)  hydrostatic  stress  contours  of 
rh  =  4g„  for  two  values  of  interface  strength.  The  amount  of 
crack  growth  is  approximately  206c  in  each  case. 


Hydrostatic  stresses  as  large  as  4g„  are  not 
achieved  in  the  case  crmax  =  5gn.  It  is  apparent 
from  Figs.  9(a)  and  9(b)  that  the  interface  strength 
has  a  major  effect  on  the  size  of  both  the  plastic 
zone  and  the  zone  of  high  hydrostatic  stress  and 
therefore  it  is  expected  to  play  a  prominent  role 
in  process  of  crack  growth. 

The  results  that  have  been  presented  thus  far 
are  based  on  the  polynomial  adhesive  relation, 
described  in  Section  2.2.2.  Calculations  were  also 
performed  using  the  exponential  adhesive  rela¬ 
tions  (31)  and  (32).  In  these  calculations,  the 
interface  strength  and  the  work  of  de-adhesion 
were  prescribed  to  have  the  same  values  as  in  the 
cases  analyzed  using  the  polynomial  relation.  Two 
values  of  interface  strength  omax  =  5g0  and  6g„ 
were  selected.  The  initial  slope  of  the  exponential 
adhesive  relation  is  higher  than  that  of  the  poly¬ 
nomial  relation  so  that  the  adhesive  stress  of  the 
former  reaches  amax  at  a  smaller  interface  dis¬ 
placement.  Since  the  loss  of  adhesion  begins  at 
smaller  interface  displacement,  the  interface 
appears  to  be  more  brittle  when  this  exponential 
relation  governs  the  de-adhesion  process.  This 
behavior  has  been  verified  by  our  calculations 
and  will  be' elaborated  upon  next. 
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The  J  integral  was  evaluated  on  paths  near  the 
growing  crack  tip  and  remote  from  it.  Near  the 
crack  tip  the  value  of  J  was  path  dependent.  On 
contours  remote  from  the  tip  (paths  whose 
characteristic  radius  exceeds  a0 )  the  value  of  J  is 
essentially  independent  of  path.  The  value  of  J  in 
the  remote  region  can  be  viewed  as  the  macro¬ 
scopic  or  global  load  parameter.  With  this  view¬ 
point,  we  employ  J  as  a  load  parameter  for 
correlating  the  process  of  crack  growth.  For  this 
purpose  the  average  of  the  J  values  in  the  remote 
region  is  referred  to  as  the  J  value. 

The  J  resistance  curves  corresponding  to  all 
three  interface  strength  levels  and  the  two  adhe¬ 
sive  relations  are  shown  in  Fig.  10.  The  trend  is 
clear.  In  all  cases  the  onset  of  crack  growth  occurs 
at  J  values  slightly  in  excess  of  <p0  and  the  resis¬ 
tance  curves  rise  rapidly  for  small  amounts  of 
growth.  The  curves  for  amax  =  5g0  level  off  after 
growth  of  about  20<5C.  In  contrast  the  resistance 
curve  for  amax  =  7g0  continues  to  rise  rapidly.  The 
behaviors  of  the  curves  for  amax  =  6 g0  fall  between 
the  two  extremes.  For  interface  strengths  less  than 
about  5 g()  the  resistance  curve  is  controlled  by  <pQ. 
These  curves  demons'rate  the  strong  role  of 
interface  strength  which  is  consistent  with  its  role 
in  establishing  the  stress  distribution  shown  in 
Fig.  8(a)  and  the  plastic  zone  size  shown  in  Fig. 
9(a).  It  can  also  be  seen  that  the  interface 
governed  by  the  exponential  relation  exhibits  a 
faster  rate  of  growth.  This  behavior  is  consistent 
with  the  loss  of  adhesion  being  achieved  at  a 
smaller  interface  displacement  for  the  exponen¬ 
tial  relation  as  has  been  noted  previously.  Never¬ 
theless  we  point  out  that  it  is  omax  that  plays  a 


more  important  role  in  determining  the  crack 
growth  resistance. 

4.2.  Intermediate-size  flaws 

Intermediate-size  flaws  are  defined  by  the 
ratios  a0/h0a  50  and  aa/6c ~  150.  In  this  section 
we  present  the  results  using  elliptical  notch  tips 
for  omax  =  lg0  although  we  have  also  performed  a 
series  of  calculations  for  a  range  of  omax  and  cir¬ 
cular  notch  tip  geometry.  Interface  profiles  are 
shown  in  Fig.  1 1(a)  at  various  levels  of  nominal 
strain  e„om,  defined  as  the  applied  displacement 
on  the  top  boundary  divided  by  the  initial  height 
of  the  crystal.  This  measure  of  the  stretch  of  the 
system  includes  contributions  from  the  bulk 
material  as  well  as  the  adhesive  layer.  Loss  of 
adhesion  begins  at  the  notch  tip  which  advances 
along  the  interface  with  increasing  remote  nomi¬ 
nal  strain.  Figure  11(b)  shows  the  distribution  of 
the  normal  stress  along  the  interface  at  the  same 


Fig.  10.  J  normalized  by  <p0  vs.  crack  extension  for  three 
interface  strengths  based  on  exponential  and  polynomial 
adhesive  relations. 


Fig.  11.  (a)  Interface  profiles  and  (b)  r2;  Kirchhoff  stress 
distributions  along  the  interface  for  intermediate-size  Haws 
(a„/dc  ~  1 50;  <rmlx  =  7g„). 
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Fia.  12.  Nominal  stress  vs.  nominal  strain  curves  for  five 
values  of  crmax  for  short  flaws  (a„/<5c  =  15). 


levels  of  nominal  strain.  The  stresses  reach  a 
maximum  at  a  short  distance  ahead  of  the  current 
tip  and  then  gradually  decrease  as  the  distance 
from  the  flaw  tip  increases.  It  should  be  noted 


that  the  stress  gradients  are  much  smaller  than 
those  of  the  long  flaw  (a0/«5c=  1500). 

4.3.  Short  flaws 

Defects  whose  ratios  aQ/h0  are  about  5  and 
a0/dc  less  than  about  50  are  considered  as  short 
flaws.  For  these,  the  role  of  interface  strength  omax 
in  the  material’s  ductility  and  failure  mode  is 
discussed  first.  Values  of  amax  in  the  range 
6g0-10g0  were  considered.  We  begin  with  the 
case  where  a0/h0  =  5  and  a0/6c~  15.  Figure  12 
shows  the  variation  in  the  nominal  stress  with 
nominal  strain  for  different  strengths  of  the  inter¬ 
face.  In  all  cases  the  nominal  stress  increases 
monotonically  to  a  maximum  and  then  decreases 
when  a  sizable  fraction  of  the  interface  has  failed. 

The  plot  shown  in  Fig.  12  reveals  an  interest¬ 
ing  phenomenon.  There  is  a  substantia]  loss  of 
ductility  when  the  interface  strength  decreases 
from  8go  to  lg0.  This  dramatic  change  in  the 
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ductility  is  due  to  a  change  in  the  failure  mode. 
For  interface  strengths  less  than  8g0  the  initial  loss 
of  adhesion  occurs  at  the  flaw  tip  which  advances 
slowly  with  increasing  remote  strain.  As  the 
remote  strain  increases,  a  zone  of  high  hydro¬ 
static  stress  develops  at  the  midligament  (the  zone 
midway  between  two  neighboring  flaws)  which 
leads  to  a  second  site  of  debonding.  In  this  region 
the  normal  stress  is  nearly  uniform  and  very 
nearly  equal  to  amax  and  the  plastic  strain  is  neg¬ 
ligibly  small.  This  stress-strain  state  is  unstable 
and  indeed  relatively  rapid  loss  of  adhesion 
proceeds  from  this  site  and  advances  towards  the 
notch  tip.  This  process  is  explained  through  the 
results  shown  in  Figs.  13(a)  and  13(b).  Figure 
13(a)  shows  the  interface  profile  for  crmax  =  7 g0  at 
various  levels  of  nominal  strain.  The  interface 
displacement  in  the  vicinity  of  the  midligament  is 
uniformly  large  and  increases  substantially  with 
increasing  nominal  strain.  Figure  13(b)  shows  the 
distribution  of  the  normal  stresses  along  the  inter¬ 
face  at  the  same  levels  of  the  nominal  strain.  The 
stresses  along  the  greater  part  of  the  ligament  are 
uniformly  high  before  debonding  and  at  the  early 
stages  of  debonding.  This  stress  pattern  is  in 
sharp  contrast  with  the  distribution  for  the  long 
and  intermediate  flaw  discussed  in  the  preceding 
sections.  The  rapid  advance  of  the  debonded 
region  towards  the  notch  tip  is  evident  from  the 
loss  of  normal  stress  shown  by  the  curve  for 
cnom  =  0.01.  It  is  note  worthy  that,  in  the  case  of 
the  weaker  interface,  crmaK  <  7gn,  the  total  nominal 
strain  at  complete  debonding  of  the  crystal  from 
the  rigid  substrate  is  not  significantly  affected  by 
interface  strength.  The  total  enom  is  about  equal  to 
(^/Lq.  This  is  not  surprising  since  the  crystal 
develops  small  plastic  strains  throughout  the 
process  of  debonding.  This  latter  aspect  is  dis¬ 
cussed  in  some  detail  in  connection  with  the  per¬ 
fect  interface  problem  in  Appendix  A. 

For  interfacial  strengths  greater  than  7gc  there 
is  only  one  site  of  debonding.  This  occurs  only  at 
the  tip  of  the  flaw  and  advances  towards  the  mid- 
ligament  in  a  gradual  manner.  This  process  is 
responsible  for  the  gradual  decrease  of  the 
nominal  stress  for  o,mx  =  8gn,  9g0  and  10g„  as 
shown  in  T  ig.  12.  The  process  of  debonding  for 
om  x  =  8g0  is  shown  in  detail  in  Fig.  13(c)  where 
the  interface  profiles  at  three  levels  of  nominal 
strain  are  plotted.  The  distribution  of  the  normal 
stresses  along  the  interface,  for  the  same  levels  of 
nominal  strain,  are  shown  in  Fig.  13(d).  The 
normal  stresses  exhibit  a  weak  maximum  near  the 


current  flaw  tip  which  is  in  contrast  with  the  stress 
pattern  in  Fig.  13(b).  It  may  be  noted  that  the 
plastic  strain  of  the  crystal  is  much  larger  for 
interfacial  strengths  in  excess  of  7g0  and  this 
results  in  higher  ductility. 

We  have  investigated  the  effects  caused  by 
varying  omM  with  <5C  fixed.  Changing  omax  under 
constant  <5C  results  in  a  change  of  the  de-adhesion 
energy  <p0.  The  effects  of  varying  crmax  and  dc  but 
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Fig.  14.  Nominal  sires'  v.s.  nominal  strain  for  short  flaws:  (a) 
<p„  =  0.22a,, g„;  (b)  </>,,  =  0.29ang(1;  (c)  the  interfaces  with 
parameters  a„/6c=13.7,  o,ma,  =  St:<i  and  n„/6c  *»  1 5.4, 
a™«  =  9gn  have  the  same  energy  of  de-adhesion 
(<?>„  =  0.3 3n08n)-  These  interfaces  have  strengths  that  differ 
11.8%  with  respect  to  their  average  value.  Thp  interfaces 
with  the  same  strength  have  energies  of  de-adhesion  that 
differ  1 1 .8%  with  respect  to  their  average  value. 
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holding  <Pq  constant  are  discussed  ne  .  Figure 
14(a)  show  the  effects  of  varying  amax  and  (5C  so  as 
to  keep  (p0  at  the  fixed  value  of  0.22 a0g0.  The 
values  crmax  =  6g0  and  7 g0  were  chosen  so  that  the 
same  failure  mode  (debonding  at  the  notch  tip 
and  at  the  midligament)  develops  in  both  cases. 
For  these  cases,  the  ductility  as  measured  by  the 
work  or  by  the  nominal  strain  to  total  separation 
are  nearly  equal.  In  Fig.  14(b)  we  compare  two 
cases  where  aniax  =  7 gQ  and  8g0.  The  lower  inter¬ 
face  strength  produces  debonding  at  two  sites 
while  the  higher  strength  causes  debonding  at  the 
notch  tip  only.  The  values  of  dc  are  such  as  to 
maintain  <p0  =  0.29a0gQ.  Very  different  ductilities 
between  the  two  cases  are  seen  although  the  work 
<p0  for  interfacial  separation  is  the  same.  Three 
cases,  which  produce  de-adhesion  only  at  the 
notch  tip,  are  shown  in  Fig.  14(c).  Two  cases 
correspond  to  omax  =  Sgn,  a0/dc  =  13.7  and 
ama,  =  9g0,  ^o/<5c  —  1 5.4;  these  have  the  same 
energy  of  de-adhesion,  i.e.  <p0  =  O.33n0g0.  The 
interface  strengths  differ  by  11.8%.  The  third 
case  with  omax  =  8gn,  fl0/t5c  ®  15.4  has  a  smaller 
value  of  cpQ  which  is  1 1.8%  smaller  than  that  of 
the  two  preceding  cases.  It  can  be  seen  that  duc¬ 
tility  is  more  greatly  affected  by  changing  omax 
(keeping  (p0  constant)  than  by  changing  <p0  (keep¬ 
ing  omax  constant). 

We  may  conclude  from  this  series  of  calcula¬ 
tions  that  the  ductility  is  strongly  affected  by  the 
interface  strength  when  de-adhesion  is  accom¬ 
panied  by  large-scale  plastic  deformation.  In  the 
context  of  the  present  model,  large-scale  plas¬ 
ticity  can  be  said  to  occur  when  the  nominal 
plastic  strain  of  the  matrix  (or  crystal)  is  compar¬ 
able  with  dJL$.  In  contrast,  if  little  plastic  defor¬ 
mation  develops  during  the  process  of 
de-adhesion,  then  the  energy  of  interfacial  sepa¬ 
ration  appears  to  dominate  the  overall  response 
of  the  material. 

The  third  series  of  calculations  was  aimed  at 
sorting  size  effects  (geometrically  similar 
notches),  relative  to  the  fracture  process  zone,  on 
ductility  and  failure  mode.  With  the  aspect  ratio 
fl0//tn  of  the  flaw  fixed  at  5,  we  carried  out  calcu¬ 
lations  for  ermax=7g0  and  a0/6c  approximately 
equal  to  15,  23  and  46.  Plots  of  nominal  stress  vs. 
nominal  strain  for  the  three  a0/bc  ratios  are 
shown  in  Fig.  15.  It  is  apparent  that  the  ductility 
decreases  as  the  fracture  process  zone  size  de¬ 
creases  (relative  to  the  defect  size),  since  the 
fracture  process  zone  size  is  scaled  by  the  work  of 
de-adhesion  (for  fixed  omM).  For  the  two  ratios 


a0/dca  15  and  23,  debonding  initiates  at  the 
notch  tip  and  the  midligament,  as  indicated  by  the 
results  in  Fig.  16(a)  for  tf0/<5c  ~  23.  for  this  failure 
mode  the  amount  of  plastic  flow  that  develops 
with  interfacial  de-adhesion  is  small  and  this 


Fig.  15.  Nominal  stress  vs.  nominal  strain  for  short  flaws. 
The  interfaces  have  the  same  strength  and  different  <;„/<),.  and 
h.,/<5c  ratios ( oma(  =  7 g„;  ujh„  =  5). 


Fig.  16.  Interface  profiles  for  short  flaws  (o„,ax  =  7g0; 
ajh „=5):  (a)  aa/6c  =»  23;  (b)  a„/6t »  46. 
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thereby  leads  to  the  rapid  drop  in  the  nominal 
stress.  In  contrast,  debonding  initiates  at  the  notch 
tip  of  the  larger  flaw  (a0jdca 46)  and  the  failure 
process  proceeds  from  this  site  as  indicated  by 
the  results  in  Fig.  16(b).  Thus,  while  the  maxi¬ 
mum  load  is  reached  at  a  smaller  nominal  strain, 
the  drop  in  the  nominal  stress  may  be  more 
gradual  than  those  for  the  two  smaller  flaws  since 
larger  amounts  of  plastic  flow  are  expected  to 
develop  as  the  debonding  process  of  the  larger 
flaw  continues. 

It  is  clear,  based  on  the  nominal  stress-strain 
curves  presented  in  this  section,  that  the  initial 


stiffness  of  the  system  increases  with  increase  in 
the  interface  strength  and  in  the  ra'io  a0/<5c. 
A  quantification  of  the  above  effect  is  presented 
in  Appendix  A,  where  an  expression  for  the 
initial  stiffness  of  a  crystal  perfectly  bonded  on 
a  rigid  substrate  is  given. 

5.  Results  using  J2  flow  theory 

5.1.  Long  flaws 

As  in  the  crystal  plasticity  calculations  the 
gei  -netry  of  the  long  flaw  is  described  by  the 
ratios  <20//20  =  500  and  a0/(5c~1500.  We  present 
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Fig.  1 7.  (a)-(cj  Crack-opening  profiles  and  (d)-(f)  t21  Kirchhoff  stress  distributions  along  the  interface,  for  three  interface 
strengths.  In  each  case,  a„/dt  *  1 500. 
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results  for  three  different  interface  strengths 
omlx  =  3.0a0,  3.5 au  and  4.0(7,,.  These  values 
roughly  correspond  to  omM  =  6.0g0,  7.0g0  and 
S.Og0  in  terms  of  the  hardness  parameter  used  in 
crystal  plasticity. 

The  deformation  produces  an  essentially  ten¬ 
sile  opening  mode  at  the  crack  tip.  Figures  17(a), 
1 7(b)  and  1 7(c)  show  the  crack  profiles  at  various 
stages  of  remote  loading  for  the  cases  where 
amax  ~  3.0O(,,  3.5cr0  and  4.0a0.  Such  profiles  are 
similar  for  all  three  cases  of  interface  strength  and 
indicate  that  the  cracks  continue  to  open,  and  to 
debond,  in  a  primarily  tensile  mode.  As  was 


noted  previously,  a  normal  displacement  <$n  =  dc 
corresponds  to  a  complete  loss  of  adhesion  and 
this  is  our  criterion  for  defining  the  current  crack 
tip  position.  De-adhesion  starts  at  a  lower  remote 
stress  for  the  lower  interface  strength  amax  =  3a0 
and  proceeds  at  a  considerably  faster  rate.  The 
crack-opening  process  was  found  to  be  smooth  in 
all  three  cases  examined  and  continued  to  occur 
only  at  the  crack  tip  as  the  remote  displacement 
increases.  It  may  be  noted  that  the  amount  of 
crack  tip  blunting  increases  with  increasing  inter¬ 
face  strength.  Figures  17(d),  17(e)  and  17(f)  show 
the  normal  traction  along  the  bond  line  for  the 
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Fig.  18.  (aHc)  Contours  of  effective  plastic  strain  fi*  =  0.001 ;  (d)-(f )  hydrostatic  stress  contours  of  t„  =  2  a„  for  different  values  of 
interface  strength.  The  amount  of  crack  growth  is  approximately  the  same  for  all  three  values  of  am!lx:  curve  ( 1)  A  .j  =  0;  curve  (2) 
A  a"  10<5c;curve(3)Afl»45<5c. 
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remote  stress  levels  corresponding  to  those  in 
Figs.  17(a),  17(b)  and  17(c).  It  is  seen  that  the 
stress  fields  essentially  translate  with  the  moving 
crack  tip  and  that  the  peak  stress  is  maintained  at 
amlx  (as  it  must  for  crack  growth  to  continue).  For 
these  interface  strengths  the  stress  gradients 
increase  slightly  as  umrix  increases. 

Figure  18  shows  effective  plastic  strain  con¬ 
tours  at  ep  =  0.001  for  three  different  levels  of 
remote  stress.  This  contour  is  effectively  the  elas¬ 
tic-plastic  boundary.  The  amount  of  crack 
growth  is  roughly  the  same  for  each  interface 
strength  considered.  In  all  cases,  yielding  is  con¬ 
fined  to  a  fraction  of  the  ligament  although  the 
deformation  is  clearly  outside  the  range  of  small- 
scale  yielding.  It  can  be  seen  that  for  the  same 
amount  of  crack  growth  the  plastic  zone  for  the 
strongest  interface  is  nearly  three  times  larger 
than  that  for  the  weakest  interface.  The  develop¬ 
ment  of  a  mode-I-like  crack  tip  field  is  evident 
and  the  leading  edge  of  the  plastic  zone  essen¬ 
tially  translates  with  the  moving  crack  tip. 

Figure  18  also  shows  contours  of  hydrostatic 
stress  for  the  same  amount  of  crack  growth  for 
the  three  interface  strengths.  The  zone  enclosed  by 
the  hydrostatic  tension  contour  given  by  rh  =  2a0 
translates  with  the  crack  tip.  For  the  weak 
strength  interfaces  the  contour  size  increases 
slightly.  In  contrast,  the  contour  size  increases 
substantially  for  the  stronger  interface  strengths. 
The  translation  and  expansion  of  the  contour  are 
a  result  of  the  stress  redistributing  from  a  high 
gradient  stationary  crack  tip  field  to  a  field  of 
weaker  stress  gradient. 

The  shape  of  the  plastic  zones  for  a  growing 
crack  suggests  that,  at  least  for  these  modest  levels 
of  crack  growth,  a  mode-I-like  crack  tip  field  is 
maintained  which  scales  approximately  with  the 
level  of  the  J  integral.  To  explore  this,  the  value  of 
J  was  computed  at  each  numerically  obtained 
crack  growth  solution  using  a  finite  deformation 
implementation  of  the  domain  integral  method 
[24],  For  each  solution,  J  values  were  extracted 
from  annular  domains  whose  mean  radii  range 
from  dL  to  KlhT.  J  values  from  domains  com¬ 
parable  in  size  with  dc  and  the  crack  tip  opening 
showed  strong  path  dependence  as  they  must, 
lleyond  these  near-tip  domains,  J  is  essentially 
path  independent,  the  variation  in  values  being 
less  than  1%.  In  subsequent  discussions  the  ./ 
value  referred  to  is  the  average  of  these  nearly 
path-independent  values,  and  this  can  serve  as  a 
load  parameter. 


Figure  19  shows  the  radial  variation  in  Kirch- 
hoff  stress  r22/o0  along  the  bond  line  for  the 
levels  of  remote  loading  indicated.  Also  shown  on 
these  plots  for  omM  =  3.0a0,  3.5 a0  and  4.0a0  are 
the  Hutchinson-Rice-Rosengren  (HRR)  fields  of 
mixity  0.97,  0.98  and  1.0  respectively.  A  mixity  of 
1.0  corresponds  to  pure  mode  I  and  departure 
from  unity  indicates  deviation  from  pure  mode  I 
conditions.  Thus  the  above  mixities  indicate 
that  the  near-tip  fields  become  more  mode  I 
like  as  the  interface  strength  increases.  This 
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Fig.  19.  Comparison  of  near-tip  fields  with  mixed-mode 
HRR  fields  for  three  interface  strengths,  (a)  am>,=  3a„. 
4/0  =  0.97;  i  -  3.5  ,Vf’  =  0.98;  (c)  oma>  =  4a,„ 
A/p  =1.0. 


behavior  is  expected  since  the  amount  of  plastic 
yielding  increases  substantially  with  interface 
strength  and  this  promotes  mode-I-like  condi¬ 
tions.  The  values  for  the  mixities  are  suggested  by 
the  small-scale  yielding  analysis  discussed  by  Shih 
and  Asaro  [5,  6]  and  Shih  et  al.  [24],  In  plotting 
the  HRR  fields  for  a  particular  extension  we  have 
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Fig.  20.  J  resistance  curves  for  a  long  flaw  for  three  interface 
strengths. 


used  the  J  value  calctAited  from  the  numerical 
fields  for  that  particular  crack  extension  and  we 
have  repositioned  the  HRR  stresses  along  the  axis 
of  the  interface  line  by  a  distance  corresponding 
to  that  amount  of  crack  extension.  We  can  see  the 
good  agreement  between  the  growing-crack  full 
field  solution  and  the  HRR  stresses  over  dis¬ 
tances  large  compared  with  the  fracture  process 
zone  size  as  scaled  by  t3c.  This  agreement  suggests 
that  the  remote  J  value  which  is  essentially  path 
independent  is  an  appropriate  load  charac¬ 
terization  parameter  for  the  near-tip  fields,  at 
least  for  modest  amounts  of  crack  extension  and 
for  the  interface  strengths  considered. 

The  crack  grows  stably  under  rising  remote 
displacement  for  the  levels  of  interface  strength 
considered  in  this  section.  To  examine  the  effect 
of  ar  ,  (the  adhesive  strength)  and  crack  tip  plas¬ 
ticity  on  fracture  toughness,  we  have  plotted  the 
computed  J  values  vs.  crack  growth  in  Fig.  20  for 
the  three  levels  of  amix  considered.  Since  J  is 
normalized  by  the  respective  <plt,  differences  in 


(x-a0)/<5c  (x-a#)/<5e 


t  ig.  2).  (a),  (b)  Flaw-opening  profiles  and  (c),  (d)  normal  Kirchhoff  stress  distributions  along  the  interface,  for  two  interface 
strengths.  In  each  case,  a„/6<  *>  1 50. 
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the  resistance  curves  are  primarily  due  to  plas¬ 
ticity.  For  example  for  crmax  -  3.0o0  the  value  /,  of 
/  at  the  onset  of  growth,  is  close  to  <p0 .  The  com¬ 
puted  J  values  rise  moderately  to  a  steady  state 
value  of  about  2.0  (p0  as  the  crack  grows.  This 
behavior  is  indeed  suggested  by  the  plastic  zone 
size  shown  in  Fig.  18(a).  For  the  omax  =  3.5o0  case, 
Jt  is  about  twice  <p0;  however,  the  J  required  for 
further  crack  growth  rises  more  rapidly.  For 
the  strongest  case  considered  (ermaj  =  4.0a0), 
/, »  2.5  ( p0 ,  which  is  larger  than  those  for  the  other 
two  cases.  Moreover  the  J  required  for  further 
crack  growth  rises  more  rapidly  than  that  for 
the  omax  =  3.5  ct0  case.  The  behavior  of  the  resis¬ 
tance  curve  for  the  strongest  interface  is  also 
suggested  by  the  large  increase  in  plastic  zone  size 
shown  in  Fig.  18(c).  The  behavior  of  Jx  and  the  JR 
curve  are  thus  strongly  affected  by  crack  tip  plas¬ 
ticity,  which  depends  on  amax. 

The  calculations  which  led  to  the  above  results 
used  the  elliptical  notch  tip  mesh  shown  in  Fig. 
5(a).  Similar  calculations  were  carried  out  using 


the  circular  notch  tip  mesh  shown  in  Fig.  5(b). 
The  differences  in  the  stress  and  deformation 
fields  and  J  resistance  curves  between  the  two 
series  of  calculations  were  not  significant,  sug¬ 
gesting  that  notch  acuity  is  not  an  important 
geometric  factor  for  long  flaws. 

5.2.  Intermediate-size  flaws 

Calculations  were  carried  out  for  an  intermedi¬ 
ate-size  defect,  where  a0/hu~  50  and  aQ/dca  150 
with  interface  strengths  of  amax  =  3o0  and  4o0. 
Figure  21  shows  the  flaw  tip  opening  profiles, 
along  with  the  normal  stresses  developed  ahead 
of  the  flaw  tip,  at  three  remote  levels  of  tensile 
stress.  It  is  clear  from  these  figures  that  the  flaws 
open  smoothly  and  quickly  develop  sharp  pro¬ 
files  much  like  those  characteristic  of  long  flaws 
as  discussed  in  the  previous  section.  The  stress 
fields  take  on  forms  that  resemble  those  of  the 
long  flaws  described  above,  with  the  stress  con¬ 
centration  refocusing  ahead  of  the  current  flaw 
tip.  However,  it  may  be  noted  that  the  stress 
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Fig.  22.  (;i),  (b)  Contours  of  effective  plastic  strain,  cp  =  0.001,  for  different  values  of  interface  strength  and  (c),  (d)  hydrostatic 
stress  contours  of  rx  =  2 au  for  different  values  of  interface  strength  at  the  same  amount  of  flaw  growth:  curve  ( 1 )  An  =  0;  curve 
(2)  A  a  a  1 34c; curve  (3)  A  a  *>  306c. 
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gradients  are  weaker  than  those  for  the  long  flaw 
{a0/dc  =  1500).  Figure  22  shows  contours  of 
effective  plastic  strain  and  hydrostatic  stress  at 
three  stages  of  flaw  growth.  We  see  that  for  the 
amax  =  3cr0  case  we  have  contained  yielding 
although  for  the  omai  =  4 o0  case  the  ligament  has 
fully  yielded  before  any  significant  flaw  growth 
takes  place.  Again  we  see  that  the  hydrostatic 
stress  contours  translate  with  the  flaw  tip  and 
increase  in  size  as  the  flaw  advances. 

The  7  integral  was  again  found  to  be  path  inde¬ 
pendent  to  within  1%  for  remote  paths  and  again 
we  address  the  question  of  7  dominance  by  com¬ 
parison  with  HRR  fields.  The  comparison  is 
shown  in  Fig.  23.  We  have  used  mode  1  HRR 
fields  for  the  comparisons  since  the  size  of  the 
plastic  zone  is  comparable  with  or  larger  than 
flaw  length.  For  the  CTmax  =  3a(l  case  the  flaw 
grows  under  contained  yielding  and  the  agree¬ 
ment  between  the  numerical  fields  and  the  HRR 
fields  seems  to  improve  after  moderate  amounts 


Fig.  23.  Comparison  of  near  tip  Fields  with  mode  I  HRR 
fields  for  two  interface  strengths:  (a)  <7m„  =  3a„;  (b) 


of  flaw  growth.  This  result  is  in  contrast  with  the 
case  for  the  stronger  interface.  Initially  the  agree¬ 
ment  between  HRR  fields  and  numerical  fields  is 
good  and  the  agreement  deteriorates  as  the  flaw 
advances.  We  observe  that  the  large  discrepancy 
between  the  flaw  tip  fields  for  the  stronger  inter¬ 
face  and  the  HRR  fields  develops  as  the  ligament 
becomes  fully  yielded. 

The  7  resistance  curves  for  both  interface 
strengths  are  shown  in  Fig.  24.  These  are  similar 
in  form  to  the  case  of  a  long  flaw.  A  relatively 
flat  resistance  curve  for  the  weak  interface,  with  7, 
close  to  (p(t  because  of  the  small  amount  of  plas¬ 
ticity  and  a  rising  resistance  curve  with  7,  ~  2.5 </>„ 
f°r  am»  =  4fl0. 

Similar  calculations  for  an  intermediate-size 
flaw  were  carried  out  using  the  circular  notch  tip 
mesh  in  Fig.  5(b).  For  the  weak  interface,  no  sig¬ 
nificant  differences  between  the  stress  and  defor¬ 
mation  fields  were  seen.  However,  for  the 
stronger  interface  ( om3x  =  4a0)  debonding  initiates 
slightly  ahead  of  the  notch  tip  and  propagates 
towards  the  midligament.  This  suggests  that 
notch  acuity  may  be  an  important  geometric 
factor  for  intermediate-size  flaws.  Notch  acuity 
also  has  a  strong  effect  on  the  behavior  of  short 
flaws  and  this  has  been  studied  by  Varias  et  al. 
(25). 

5.3.  Short  flaws 

Figure  25  shows  the  opening  profiles  and 
normal  stresses  along  the  bond  line  for  a  short 
flaw,  where  a0/h0=5,  an/(\  ~  15  and  omax  =  3o() 
and  4o0.  It  can  be  seen  that  the  stress  fields 
do  not  resemble  the  fields  ahead  of  a  sharp 
crack.  There  is  very  little  stress  concentration  at 
the  tip  of  the  flaw  and  the  normal  stress  is  nearly 
uniform  along  the  ligament  and  its  magnitude  is 


Fig.  24.  J  resistance  curves  for  an  intermediate-size  (taw  for 
two  interface  strengths. 
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Fie.  25.  :a  .  b  ■;  Interface  profiles  and  (c).  (d)  normal  Kirehhoff  stress  distributions  along  the  interface,  for  two  interface  strengths. 
In  each  case.  <r|p/rV  =  I  5. 


close  to  amax.  Furthermore  the  J  integral  is 
strongly  path  dependent.  This  behavior  is  true  for 
both  interface  strengths  considered.  It  can  be 
seen  from  the  interface  profiles  that,  while 
debonding  initiates  at  the  flaw  tip,  further  de¬ 
adhesion  is  not  confined  to  the  tip.  When  maxi¬ 
mum  load  is  reached,  the  entire  interface  has 
undergone  significant  de-adhesion  and  the 
remote  load  then  drops  abruptly.  We  conclude 
that  the  short  flaw  does  not  behave  like  a  crack. 

6.  Discussion  and  conclusions 

The  adhesive  relation  of  the  exponential  and 
polynomial  type  is  completely  specified  by  the 
de-adhesion  energy  (p0  and  the  interface  strength 
omax.  The  adhesive  relation  has  a  characteristic 
length  <5C  which  is  determined  by  <p0  and  crmax.  For 
example,  dc  is  given  by  16<pn/9<7max  for  the  poly¬ 
nomial  relation.  Thus,  <5C  and  <7max  also  specify  the 
adhesive  relation.  To  varying  degrees,  all  three 


parameters  affect  the  ductility  and  the  failure 
modes  of  the  system.  For  example,  in  the  case  of 
the  homogeneous  block  debonding  from  a  rigid 
substrate,  as  discussed  in  Appendix  A,  the  level 
of  matrix  plastic  strains  that  are  achieved  during 
debonding  compared  with  dJL^  determines  how 
the  adhesive  parameters  control  ductility.  When 
the  nominal  strain  enom  of  the  system  at  separa¬ 
tion  is  comparable  with  dJL^,  ductility  is  con¬ 
trolled  by  cp().  For  stronger  interfaces,  where  cm„„ 
becomes  large  compared  w'ith  d^/L^,  om.M  plays 
the  dominant  role.  This  example  has  direct  rele¬ 
vance  to  the  more  complex  cases  involving  the 
separation  of  interfaces  containing  periodic 
arrays  of  short  flaws  with  a0fh0  approximately  5 
as  discussed  in  Section  4.3.  In  these  cases,  both 
the  process  of  failure  ( i.e .  the  failure  mode)  and 
the  resulting  toughness  as  measured  by  the  total 
work  expended  on  the  flawed  body  can  be 
affected.  We  have  shown  that,  when  omax  is  large 
relative  to  the  material's  yield  strength  or  when 
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the  flaw  size  is  large  compared  with  <5C,  the  pro¬ 
cess  of  de-adhesion  initiates  and  proceeds  from  a 
single  site  which  is  the  tip  of  the  flaw.  On  the 
contrary,  when  the  relative  flaw  size  a0/dc  is 
smaller  than  40,  and  omiX  is  small,  de-adhesion 
can  occur  at  two  sites:  at  the  flaw  tip  and  at  the 
midligament.  For  even  blunter  short  flaws,  de¬ 
adhesion  may  initiate  only  at  the  midligament 
[25], 

In  the  case  of  long  flaws  the  plastic  strains  well 
within  the  flaw  tip  plastic  zone  are  large  and  the 
failure  behavior  is  then  controlled  by  the  plastic 
zone  size  relative  to  <5C.  For  example,  the  fracture 
toughness  and  the  growth  resistance  are  con¬ 
trolled  primarily  by  %  when  the  plastic  zone  size 
during  the  process  of  flaw  growth  is  comparable 
with  dc.  However,  when  the  interface  is  suffi¬ 
ciently  strong  that  the  plastic  zone  that  develops 
at  the  onset  of  flaw  growth  is  large  compared 
with  6C,  it  is  the  interface  strength  amM  that 
controls  the  initiation  toughness  and  flaw 
growth  resistance  curve.  Such  behavior  under 
large-scale  plasticity  can  be  clearly  seen  in  Fig.  10 
for  the  case  of  crystal  plasticity  and  in  Fig.  20 
where  the  material  is  modeled  by  J2  flow  theory. 

We  observe  that  the  long  flaw'  (a0/dc  «  1500) 
behaves  like  a  crack  for  the  range  of  interface 
strengths  and  load  levels  considered  for  both 
material  characterizations.  For  certain  interface 
strengths  and  load  levels,  the  fields  of  the  inter¬ 
mediate-size  flaw  (a0/< 5C  ~  150)  have  some  crack¬ 
like  features.  In  the  case  of  the  short  flaw  («„/bc 
<50),  crack-like  behavior. was  not  observed  for 
any  interface  strength  or  load  level.  We  reached 
these  conclusions  through  direct  comparisons  of 
the  full  field  solutions  for  stationary  and  growing 
flaw  with  singular  fields  of  the  HRR  type.  Specifi¬ 
cally,  the  stationary  interface  long  flaw  fields 
in  the  annular  zone,  which  encompass  the  finite 
strain  zone  and  the  fracture  process  zone,  are 
well  approximated  by  mixed-mode  HRR  fields. 
Furthermore,  the  annular  zone  of  high  stresses 
essentially  translates  with  the  moving  long  flaw  tip 
and  the  fields  remain  in  agreement  with  the  HRR 
fields  for  limited  amounts  of  growth.  Such 
observations  allow  us  to  conclude  that  the  near- 
lip  state  is  characterized  by  J  and  the  mode 
mbeity.  Recent  results  of  Needleman  [26]  aTe  con¬ 
sistent  with  our  results  described  above.  For 
example  in  his  studies  of  de-adhesion  along  an 
imperfect  interface  he  observes  the  transition 
from  a  more  or  less  uniform  separation  profile  to 
crack-like  profiles  as  the  ratio  of  defect  length  to 


the  characteristic  dimension  6C  increases.  For  the 
longest  defect  considered,  i.e.  a0/dc  =  7500,  the 
shape  of  his  plastic  zone  contour  is  consistent 
with  the  tension-dominated  plastic  zone  contour 
given  by  Shih  and  Asaro  [4], 

We  must  point  out  that  the  boundary  value 
problem  dealt  with  in  this  paper  gives  rise  to  very 
high  triaxial  stress  states  over  the  entire  ligament 
as  the  deformation  progresses  into  the  fully 
yielded  state.  Such  stress  states  may  not  be  typical 
of  load  states  encountered  in  actual  systems. 
Investigations  for  a  plastically  deformable  mate¬ 
rial  on  an  elastic  substrate  using  boundary 
conditions  w'hich  lead  to  less  severe  triaxiality  are 
in  progress. 
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Appendix  A 

A 1.  Failure  of  a  perfect  interface 

Insight  into  the  debonding  process  can  be 
gained  by  examining  a  deformable  crystal  which 
is  perfectly  bonded  on  a  rigid  substrate.  We  con¬ 
sider  a  rectangular  block  which  is  perfectly 
bonded  at  its  bottom  edge  to  a  rigid  substrate  and 
subject  to  the  boundary  conditions  shown  in  Fig. 
A 1 .  The  block  deforms  under  plane  strain  condi¬ 
tions.  As  indicated  in  the  figure,  the  crystal  has 
two  slip  systems,  which  are  symmetric  with 
respect  to  the  direction  of  the  applied  load. 
Because  of  the  boundary  conditions  imposed,  a 
slip  system  with  direction  along  the  X j  axis  (see 


t  j  —  0  ,du2>0 


Fig.  Al.  Boundary  value  problem  of  a  perfect  interface,  u, 
and  t,  are  the  components  of  the  displacements  and  tractions 
respectively. 


f 


Fig.  1)  cannot  operate.  To  make  the  analysis 
tractable,  we  assume  a  rate-independent  plastic 
response  where  slip  is  possible  when  the  resolved 
shear  stress  r{o)  reaches  a  critical  value,  which  is 
equal  to  the  hardness  g(ya)  of  the  crystal.  The 
results  that  we  present  are  not  particular  to  the 
details  of  the  constitutive  response  and  they  are 
representative  of  the  rate-dependent  material 
response  of  the  type  considered  in  the  main 
sections  of  this  paper.  The  elastic  response  of  the 
crystal  is  assumed  to  be  isotropic. 

The  boundary  conditions  that  are  applied  to 
the  deformable  block  are  similar  to  those  applied 
to  the  model  boundary  value  problem  shown  in 
Fig.  4.  The  displacement  of  the  top  boundary'  is 
monotonically  increasing  while  the  left  and  the 
right  boundaries  are  held  fixed.  These  three 
boundaries  are  free  of  tangential  tractions.  The 
bottom  boundary  is  held  to  the  rigid  substrate  by 
adhesive  tractions.  Since  the  displacement  along 
the  A',  direction  is  zero,  only  normal  adhesive 
tractions  act  across  the  interface,  which  obey  the 
general  relation  o„  =  f(d),  where  6  is  the  normal 
displacement  of  the  bottom  boundary.  The  ana¬ 
lytic  solution  for  infinitesimal  strains  is  discussed. 

The  initial  stiffness  of  the  system  is  given  by 
the  relation 


da 

,[ii 

fl-v-2-T  1  1 

d  (u/U) 

[t/“0 

l  1-v  /  '  w'( 0) 

(Al) 


where  E  and  v  are  Young's  modulus  and  Pois¬ 
son's  ratio  respectively,  and  /'( 0)  is  the  slope  of 
the  adhesive  relation  at  zero  displacement.  If  the 
product  L0/'( 0)  of  the  height  of  the  block  and  the 
initial  stiffness  is  large  compared  with  £,  the 
system  responds  as  if  the  adhesion  is  infinitely 
strong.  Since  /'( 0)  scales  with  the  ratio  omM/dc, 
an  increase  in  the  interface  strength  or  decrease 
in  the  ratio  dc/L0  results  in  a  stiffer  initial 
response  of  the  system.  This  behavior  is  also 
observed  in  the  short-flaw  problem. 

Owing  to  the  symmetry  conditions,  the 
resolved  shear  stress  of  both  slip  systems  has  the 
same  absolute  value  and  so  does  plastic  slip  along 
the  two  slip  directions.  Initial  yielding  occurs 
when  r(a|  reaches  the  initial  hardness  g(,  of  the 
crystal,  or  equivalently  when  the  stress  o[-o22) 
reaches  the  value  a0,  i.e. 

2(1  ~v)  go 
°  (l-.2v)  sin20 


t 
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The  total  strain  increment  is 


of  the  system  is  given  by 


d£*-*de/  +  de/ 

=  de,/  +  X  dy(o)J(/ni(o>s/o)  +  m}a)siM) 
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where 
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d£22  = — : - +  dy  sin  26 

To 

dfn  =  0  =  d£e,,  -  dy  sin  28 
d£3,  -  d£e33  =  dfp33  =  0 
dy,  =  ~dy2  =  dy 

It  should  be  noted  that  the  shearing  strains  are 
zero.  The  resolved  shear  stress  on  the  first  slip 
system  is  given  by 

T  =  g(  7a)  =  <?( 2 y)  =  K -  CT,  I )  sin  2  <9 

A  similar  relation  holds  for  r2.  To  determine  the 
stresses,  the  above  relations  must  be  supple¬ 
mented  by  the  elastic  stress-strain  relations. 

Combining  the  above  relations  and  integrating 
the  resulting  expression,  we  get  the  relation  for 
the  stress  a  as  a  function  of  the  slip  of  the  first 
system  y,  i.e. 


where  u  is  the  displacement  of  the  top  boundary 
of  the  crystal  and  <5m  is  the  normal  interface  dis¬ 
placement  when  the  normal  adhesive  traction  is 
at  its  maximum  value  of  amax.  The  above  relations 
show  that  the  nominal  strain  at  the  peak  stress 
decreases  when  amM  or  <3 m/L0  decreases. 

After  reaching  the  peak  stress  the  crystal 
unloads  elastically  as  the  process  ot  de-adhesion 
proceeds.  The  crystal  is  completely  debonded 
from  the  substrate  when  its  bottom  boundary 
displacement  d  reaches  a  value  equal  to  <5C.  The 
residual  stress  upon  complete  debonding  is  given 
by  the  difference  between  the  elastic  and  the  elas¬ 
tic-plastic  solutions  for  the  system  when  the 
applied  stress  is  of  magnitude  omM.  The  total 
engineering  strain  at  complete  debonding  is 
determined  by  the  permanent  plastic  strain  and 
the  elastic  strain  associated  with  the  residual 
stress,  i.e. 
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ym.,x  sin  20  (A5) 


The  crystal  stress  o  reaches  a  peak  value  of 
(which  is  the  strength  of  the  interface).  At  this 
point  the  plastic  slip  y  also  reaches  its  maximum 
value  designated  by  ym3X.  Thereafter  its  value 
remains  constant.  Under  the  assumption  that  ymax 
is  comparable  with  the  elastic  strains,  we  obtain 
the  result 


1  — v  — 2 v2  om3X  2(1  -  v2)  gu 
\  sin  2 6  E  (sin  20)2  E  j 

(  (sin  2df  E j 


(A3) 


The  above  relation  is  obtained  by  assuming  that 
ovV,/  is  given  by  eqn.  (ob)  and  taking  the  asymp¬ 
totic  expansion  for  small  ya.  A  numerical  solution 
of  (A2)  has  shown  that  (A3)  is  rather  accurate. 
This  is  not  surprising  because  of  the  character  of 
the  highly  constrained  deformation  caused  by  the 
Prescribed  boundary  conditions. 

At  the  peak  stress  omtx,  the  nominal  strain  enom 


The  first  term  on  the  right-hand  side  of  (A5)  is 
the  contribution  of  the  adhesive  relation.  The 
second  term  is  due  to  the  plastic  deformation  of 
the  crystal  along  the  X2  direction  and  the  third 
term  is  the  elastic  strain  due  to  the  a, ,  and  a-,, 
residual  stresses.  The  last  two  terms  vanish  when 
the  interface  strength  is  smaller  than  the  con¬ 
strained  yield  strength  o0  of  the  crystal.  Thus, 
when  the  plastic  deformation  of  the  crystal  is 
small,  the  engineering  strain  »/L0  of  the  system  at 
complete  debonding  is  nearly  given  by  <5 C/L„. 
This  result  was  observed  in  the  short-defect 
problem  where  the  interface  strength  is  smaller 
than  8g0. 

At  complete  debonding,  the  work  done  by  the 
applied  load  is 

W-J/Wdd  +  Lo  J  g(ya)dya 
0  0 
£7 

+  2(i-2)(^  sin  26 f  (A6) 
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where 

w  2r 

7^=  j  g(ya)dy*  “2g0ymax  (A7) 

^0  0 

The  above  relation  for  the  plastic  work  is  valid 
for  plastic  slip  which  is  comparable  with  the 
elastic  strain. 

Equation  (A6)  shows  that  the  work  done  by 
the  applied  load  is  the  sum  of  the  work  of  de¬ 
adhesion,  the  work  of  plastic  deformation  of  the 
crystal  and  the  elastic  energy  stored  in  the  crystal 
(due  to  the  residual  stresses).  The  last  two  terms 
vanish  when  the  crystal  deforms  only  elastically 
so  that  the  work  of  the  applied  load  is  totally 
spent  in  the  de-adhesion  process.  Therefore  the 
ductility  of  the  system  is  controlled  by  the  work  of 
de-adhesion,  when  the  plastic  strain  of  he  crystal 
is  small.  In  contrast,  the  ductility  depends  on  both 
the  work  of  de-adhesion  and  the  work  of  plastic 
deformation,  when  the  plastic  strain  of  the  crystal 
is  comparable  with  <3 JL0.  The  plastic  work  at 
system  failure  depends  primarily  on  interface 
strength.  In  fact  the  interface  strength  is  the  only 
parameter  of  the  adhesive  relation  governing  the 
plastic  work  to  failure  in  this  problem.  In  other 
words,  the  ductility  depends  on  interface 
strength.  This  simple  problem  is  similar  in  some 
respects  with  the  short-defect  problems  discussed 
in  Section  4.3  and  provides  insight  into  those 
problems.  It  also  explains  the  observations  w'ith 
regard  to  the  mechanical  response  of  a  system 
containing  a  planar  interface  populated  by  a 
periodic  array  of  cavities  [AlJ. 

The  solution  to  the  present  boundary  value 
problem  assumes  an  even  simpler  form  when  J2 
flow  theory  is  used.  In  this  theory,  the  material 
hardens  according  to  d  =  H(i p).  Initial  yielding 
occurs  when  d=  o,„  i.e.  when  the  applied  stress  is 
given  by 
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After  yielding  has  occurred,  the  strain  increments 
are  given  by  the  relations 


drs,  = 


di/  -  dd 


■  dc%,  +  dep 


dcM  =  dc,,  =  0  =  dccu  -\dep 
d£en  =  dce3, 


It  should  be  noted  that  all  shear  strains  vanish.  We 
also  have  the  relation 

d=  a  -  an  =  H(ep) 

Combining  the  above  relations  and  integrating 
the  resulting  expression,  we  obtain 

o  =  — - — {i£eP  +  (l- v)H(eP)}  (A8) 

1  -  2v 


following  the  procedure  described  in  the  crystal 
plasticity  case,  similar  relations  for  engineering 
strain  at  peak  stress,  total  engineering  strain  and 
ductility  are  derived: 
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(All) 


Finally  we  note  that  the  perfect  interface 
problem  does  not  have  always  an  equilibrium 
solution  during  unloading.  The  existence  of  an 
equilibrium  solution  depends  on  the  interface  and 
the  material  properties.  During  unloading,  the 
bottom  boundary  displacement  increment  is 
larger  than  that  prescribed  on  the  top  boundary, 
i.e. 


0<dn<dd 


The  above  inequalities,  combined  with  the 
stress-strain  relation  during  unloading,  result  in 
the  following  necessary  condition  for  the  exis¬ 
tence  of  equilibrium  solution: 


„  df  1-v  1 

d<)  1  -  v-  2v2  L,, 


(A12) 


The  left  inequality  is  automatically  satisfied  in  the 
unloading  part  of  the  adhesive  relations.  The 
right  inequality  gives  a  lower  bound  for  the  tan¬ 
gent  of  the  adhesive  relation  during  unloading. 


The  issue  of  existence  of  an  equilibrium  solution 
also  arises  in  problems  of  debonding  of  imperfect 
interfaces,  which  were  discussed  in  the  main 
sections  of  this  paper.  It  is  clear  that  the  existence 
of  an  equilibrium  solution  can  be  given  in  a 
manner  similar  to  (A12).  However,  an  explicit 
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relation  cannot  be  found  owing  to  the  complexity 
of  the  boundary  value  problem. 

Reference  for  Appendix  A 

At  A.  G.  Varias,  R.  J.  Asaro  and  C.  F.  Shih,  to  bo  published. 
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